Source author record

Felix Otto

Felix Otto appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

37works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

37 published item(s)

preprint2022arXiv

A tree-free approach to regularity structures: The regular case for quasi-linear equations

We give a motivation and gentle introduction into the regularity structure and model introduced by Otto, Sauer, Smith and Weber, which fall into the framework of Hairer, but have a greedier index set than the one given by trees. We do this here for a simple quasi-linear parabolic equation and assume that the driving noise is so regular that no renormalization is needed. We introduce the abstract model space $\mathsf{T}$ and its grading, the pre-model $\mathbfΠ$, the centered model $Π_x$, the structure group $\mathsf{G}$, and the re-centering transformations $Γ_{xy}$. Using integration and reconstruction, we establish the desired estimates on $Π_x$ and $Γ_{xy}$, which here are deterministic since we deal with the regular case.

preprint2022arXiv

Bias in the representative volume element method: periodize the ensemble instead of its realizations

We study the Representative Volume Element (RVE) method, which is a method to approximately infer the effective behavior $a_{\text{hom}}$ of a stationary random medium. The latter is described by a coefficient field $a(x)$ generated from a given ensemble $\langle\cdot\rangle$ and the corresponding linear elliptic operator $-\nabla\cdot a\nabla$. In line with the theory of homogenization, the method proceeds by computing $d = 3$ correctors (d denoting the space dimension).To be numerically tractable, this computation has to be done on a finite domain: the so-called "representative" volume element, i. e. a large box with, say, periodic boundary conditions. The main message of this article is: Periodize the ensemble instead of its realizations. By this we mean that it is better to sample from a suitably periodized ensemble than to periodically extend the restriction of a realization $a(x)$ from the whole-space ensemble $\langle\cdot\rangle$. We make this point by investigating the bias (or systematic error), i. e. the difference between $a_{\text{hom}}$ and the expected value of the RVE method, in terms of its scaling w. r. t. the lateral size $L$ of the box. In case of periodizing $a(x)$, we heuristically argue that this error is generically $O(L^{-1})$. In case of a suitable periodization of $\langle\cdot\rangle$, we rigorously show that it is $O(L^{-d})$. In fact, we give a characterization of the leading-order error term for both strategies, and argue that even in the isotropic case it is generically non-degenerate. We carry out the rigorous analysis in the convenient setting of ensembles $\langle\cdot\rangle$ of Gaussian type with integrable covariance, which allow for a straightforward periodization and which make the Price theorem and the Malliavin calculus available for optimal stochastic estimates of correctors.

preprint2022arXiv

Thermodynamically consistent and positivity-preserving discretization of the thin-film equation with thermal noise

In micro-fluidics not only does capillarity dominate but also thermal fluctuations become important. On the level of the lubrication approximation, this leads to a quasi-linear fourth-order parabolic equation for the film height $h$ driven by space-time white noise. The gradient flow structure of its deterministic counterpart, the thin-film equation, which encodes the balance between driving capillary and limiting viscous forces, provides the guidance for the thermodynamically consistent introduction of fluctuations. We follow this route on the level of a spatial discretization of the gradient flow structure. Starting from an energetically conformal finite-element (FE) discretization, we point out that the numerical mobility function introduced by Grün and Rumpf can be interpreted as a discretization of the metric tensor in the sense of a mixed FE method with lumping. While this discretization was devised in order to preserve the so-called entropy estimate, we use this to show that the resulting high-dimensional stochastic differential equation (SDE) preserves pathwise and pointwise strict positivity, at least in case of the physically relevant mobility function arising from the no-slip boundary condition. As a consequence, this discretization gives rise to a consistent invariant measure, namely a discretization of the Brownian excursion (up to the volume constraint), and thus features an entropic repulsion. The price to pay over more naive discretizations is that when writing the SDE in Itô's form, which is the basis for the Euler-Mayurama time discretization, a correction term appears. To conclude, we perform various numerical experiments to compare the behavior of our discretization to that of the more naive finite difference discretization of the equation.

preprint2021arXiv

Sharp boundary $\varepsilon$-regularity of optimal transport maps

In this paper we develop a boundary $\varepsilon$-regularity theory for optimal transport maps between bounded open sets with $C^{1,α}$-boundary. Our main result asserts sharp $C^{1,α}$-regularity of transport maps at the boundary in form of a linear estimate under certain assumptions: The main quantitative assumptions are that the local nondimensionalized transport cost is small and that the boundaries are locally almost flat in $C^{1,α}$. Our method is completely variational and builds on the recently developed interior regularity theory.

preprint2020arXiv

On the existence of the Green function for elliptic systems in divergence form

We study the existence of the Green function for an elliptic system in divergence form $-\nabla\cdot a\nabla$ in $\mathbb{R}^d$, with $d>2$. The tensor field $a=a(x)$ is only assumed to be bounded and $λ$-coercive. For almost every point $y \in \mathbb{R}^d$, the existence of a Green's function $G(a; \cdot, y)$ centered in $y$ has been proven in [J. Conlon, A. Giunti and F.Otto, "Green's function for elliptic systems: Delmotte-Deuschel bounds", 2017]. In this paper, we show that the set of points $y \in \mathbb{R}^d$ for which $G(a; \cdot, y)$ does not exist has zero $p$-capacity, for an exponent $p >2$ depending only on the dimension $d$ and the ellipticity ratio of $a$.

preprint2020arXiv

The annealed Calderon-Zygmund estimate as convenient tool in quantitative stochastic homogenization

This article is about the quantitative homogenization theory of linear elliptic equations in divergence form with random coefficients. We derive gradient estimates on the homogenization error, i.e. on the difference between the actual solution and the two-scale expansion of the homogenized solution, both in terms of strong norms (oscillation) and weak norms (fluctuation). These estimates are optimal in terms of scaling in the ratio between the microscopic and the macroscopic scale. The purpose of this article is to highlight the usage of the recently introduced annealed Calderon-Zygmund (CZ) estimates in obtaining the above, previously known, error estimates. Moreover, the article provides a novel proof of these annealed CZ estimate that completely avoids quenched regularity theory, but rather relies on functional analysis. It is based on the observation that even on the level of operator norms, the Helmholtz projection is close to the one for the homogenized coefficient (for which annealed CZ estimates are easily obtained). In this article, we strive for simple proofs, and thus restrict ourselves to ensembles of coefficient fields that are local transformations of Gaussian random fields with integrable correlations and Hölder continuous realizations. As in earlier work, we use the natural objects from the general theory of homogenization, like the (potential and flux) correctors and the homogenization commutator. Both oscillation and fluctuation estimates rely on a sensitivity calculus, i.e. on estimating how sensitively the quantity of interest does depend on an infinitesimal change in the coefficient field, which is fed into the Spectral Gap inequality. In this article, the annealed CZ estimate is the only form in which elliptic regularity theory enters.

preprint2020arXiv

Uniform energy distribution in a pattern-forming system of surface charges

We consider a variational model for a charge density $u\in\{-1,1\}$ on a (hyper)plane, with a short-range attraction coming from the interfacial energy and a long-range repulsion coming from the electrostatic energy. This competition leads to pattern formation. We prove that the interfacial energy density is (asymptotically) equidistributed at scales large compared to the scale of the pattern. We follow the strategy laid out in [G. Alberti, R. Choksi, F. Otto, Uniform energy distribution for an isoperimetric problem with long-range interactions, J. A.M.S.]. The challenge comes from the reduced screening capabilities of surface charges compared to the volume charges considered in that paper.

preprint2019arXiv

Quantitative estimates in stochastic homogenization for correlated coefficient fields

This paper is about the homogenization of linear elliptic operators in divergence form with stationary random coefficients that have only slowly decaying correlations. It deduces optimal estimates of the homogenization error from optimal growth estimates of the (extended) corrector. In line with the heuristics, there are transitions at dimension $d=2$, and for a correlation-decay exponent $β=2$; we capture the correct power of logarithms coming from these two sources of criticality. The decay of correlations is sharply encoded in terms of a multiscale logarithmic Sobolev inequality (LSI) for the ensemble under consideration --- the results would fail if correlation decay were encoded in terms of an $α$-mixing condition. Among other ensembles popular in modelling of random media, this class includes coefficient fields that are local transformations of stationary Gaussian fields. The optimal growth of the corrector $ϕ$ is derived from bounding the size of spatial averages $F=\int g\cdot\nablaϕ$ of its gradient. This in turn is done by a (deterministic) sensitivity estimate of $F$, that is, by estimating the functional derivative $\frac{\partial F}{\partial a}$ of $F$ w.~r.~t.~the coefficient field $a$. Appealing to the LSI in form of concentration of measure yields a stochastic estimate on $F$. The sensitivity argument relies on a large-scale Schauder theory for the heterogeneous elliptic operator $-\nabla\cdot a\nabla$. The treatment allows for non-symmetric $a$ and for systems like linear elasticity.

preprint2016arXiv

A Liouville theorem for elliptic systems with degenerate ergodic coefficients

We study the behavior of second-order degenerate elliptic systems in divergence form with random coefficients which are stationary and ergodic. Assuming moment bounds like Chiarini and Deuschel [Arxiv preprint 1410.4483, 2014] on the coefficient field $a$ and its inverse, we prove an intrinsic large-scale $C^{1,α}$-regularity estimate for $a$-harmonic functions and obtain a first-order Liouville theorem for subquadratic $a$-harmonic functions.

preprint2016arXiv

Convergence of the thresholding scheme for multi-phase mean-curvature flow

We consider the thresholding scheme, a time discretization for mean curvature flow introduced by Merriman, Bence and Osher. We prove a convergence result in the multi-phase case. The result establishes convergence towards a weak formulation of mean curvature flow in the BV-framework of sets of finite perimeter. The proof is based on the interpretation of the thresholding scheme as a minimizing movement scheme by Esedoglu et. al.. This interpretation means that the thresholding scheme preserves the structure of (multi-phase) mean curvature flow as a gradient flow w. r. t. the total interfacial energy. More precisely, the thresholding scheme is a minimizing movement scheme for an energy functional that $Γ$-converges to the total interfacial energy. In this sense, our proof is similar to the convergence results of Almgren, Taylor and Wang and Luckhaus and Sturzenhecker, which establish convergence of a more academic minimizing movement scheme. Like the one of Luckhaus and Sturzenhecker, ours is a conditional convergence result, which means that we have to assume that the time-integrated energy of the approximation converges to the time-integrated energy of the limit. This is a natural assumption, which however is not ensured by the compactness coming from the basic estimates.

preprint2016arXiv

Correlation structure of the corrector in stochastic homogenization

Recently, the quantification of errors in the stochastic homogenization of divergence-form operators has witnessed important progress. Our aim now is to go beyond error bounds, and give precise descriptions of the effect of the randomness, in the large-scale limit. This paper is a first step in this direction. Our main result is to identify the correlation structure of the corrector, in dimension $3$ and higher. This correlation structure is similar to, but different from that of a Gaussian free field.

preprint2016arXiv

Green's function for elliptic systems: existence and Delmotte-Deuschel bounds

We prove that for an open domain $D \subset \mathbb{R}^d $ with $d \geq 2 $ , for every (measurable) uniformly elliptic tensor field $a$ and for almost every point $y \in D$ , there exists a unique Green's function centred in $ y $ associated to the vectorial operator $ -\nabla \cdot a\nabla $ in D. In particular, when $d > 2$ this result also implies the existence of the fundamental solution for elliptic systems, i.e. the Green function for $ -\nabla \cdot a\nabla $ in $ \mathbb{R}^d $. Moreover, introducing an ensemble $\langle\cdot \rangle$ over the set of uniformly elliptic tensor fields, under the assumption of stationarity we infer for the fundamental solution $G$ some pointwise bounds for $\langle |G(\cdot; x,y)|\rangle$, $\langle|\nabla_x G(\cdot; x,y)|\rangle$ and $\langle |\nabla_x\nabla_y G(\cdot; x,y)|\rangle$. These estimates scale optimally in space and provide a generalization to systems of the bounds obtained by Delmotte and Deuschel for the scalar case.

preprint2016arXiv

Stochastic homogenization of linear elliptic equations: Higher-order error estimates in weak norms via second-order correctors

We are concerned with the homogenization of second-order linear elliptic equations with random coefficient fields. For symmetric coefficient fields with only short-range correlations, quantified through a logarithmic Sobolev inequality for the ensemble, we prove that when measured in weak spatial norms, the solution to the homogenized equation provides a higher-order approximation of the solution to the equation with oscillating coefficients. In the case of nonsymmetric coefficient fields, we provide a higher-order approximation (in weak spatial norms) of the solution to the equation with oscillating coefficients in terms of solutions to constant-coefficient equations. In both settings, we also provide optimal error estimates for the two-scale expansion truncated at second order. Our results rely on novel estimates on the second-order homogenization corrector, which we establish via sensitivity estimates for the second-order corrector and a large-scale $L^p$ theory for elliptic equations with random coefficients. Our results also cover the case of elliptic systems.

preprint2016arXiv

The corrector in stochastic homogenization: optimal rates, stochastic integrability, and fluctuations

We consider uniformly elliptic coefficient fields that are randomly distributed according to a stationary ensemble of a finite range of dependence. We show that the gradient and flux $(\nablaϕ,a(\nabla ϕ+e))$ of the corrector $ϕ$, when spatially averaged over a scale $R\gg 1$ decay like the CLT scaling $R^{-\frac{d}{2}}$. We establish this optimal rate on the level of sub-Gaussian bounds in terms of the stochastic integrability, and also establish a suboptimal rate on the level of optimal Gaussian bounds in terms of the stochastic integrability. The proof unravels and exploits the self-averaging property of the associated semi-group, which provides a natural and convenient disintegration of scales, and culminates in a propagator estimate with strong stochastic integrability. As an application, we characterize the fluctuations of the homogenization commutator, and prove sharp bounds on the spatial growth of the corrector, a quantitative two-scale expansion, and several other estimates of interest in homogenization.

preprint2015arXiv

A higher-order large-scale regularity theory for random elliptic operators

We develop a large-scale regularity theory of higher order for divergence-form elliptic equations with heterogeneous coefficient fields $a$ in the context of stochastic homogenization. The large-scale regularity of $a$-harmonic functions is encoded by Liouville principles: The space of $a$-harmonic functions that grow at most like a polynomial of degree $k$ has the same dimension as in the constant-coefficient case. This result can be seen as the qualitative side of a large-scale $C^{k,α}$-regularity theory, which in the present work is developed in the form of a corresponding $C^{k,α}$-"excess decay" estimate: For a given $a$-harmonic function $u$ on a ball $B_R$, its energy distance on some ball $B_r$ to the above space of $a$-harmonic functions that grow at most like a polynomial of degree $k$ has the natural decay in the radius $r$ above some minimal radius $r_0$. Though motivated by stochastic homogenization, the contribution of this paper is of purely deterministic nature: We work under the assumption that for the given realization $a$ of the coefficient field, the couple $(ϕ,σ)$ of scalar and vector potentials of the harmonic coordinates, where $ϕ$ is the usual corrector, grows sublinearly in a mildly quantified way. We then construct "$k$th-order correctors" and thereby the space of $a$-harmonic functions that grow at most like a polynomial of degree $k$, establish the above excess decay and then the corresponding Liouville principle.

preprint2015arXiv

An isoperimetric problem with Coulomb repulsion and attraction to a background nucleus

We study an isoperimetric problem the energy of which contains the perimeter of a set, Coulomb repulsion of the set with itself, and attraction of the set to a background nucleus as a point charge with charge $Z$. For the variational problem with constrained volume $V$, our main result is that the minimizer does not exist if $V - Z$ is larger than a constant multiple of $\max(Z^{2/3}, 1)$. The main technical ingredients of our proof are a uniform density lemma and electrostatic screening arguments.

preprint2015arXiv

Anchored Nash inequalities and heat kernel bounds for static and dynamic degenerate environments

We introduce anchored versions of the Nash inequality. They allow to control the $L^2$ norm of a function by Dirichlet forms that are not uniformly elliptic. We then use them to provide heat kernel upper bounds for diffusions in degenerate static and dynamic random environments. As an example, we apply our results to the case of a random walk with degenerate jump rates that depend on an underlying exclusion process at equilibrium.

preprint2015arXiv

Branched Microstructures in the Ginzburg-Landau Model of Type-I Superconductors

We consider the Ginzburg-Landau energy for a type-I superconductor in the shape of an infinite three-dimensional slab, with two-dimensional periodicity, with an applied magnetic field which is uniform and perpendicular to the slab. We determine the optimal scaling law of the minimal energy in terms of the parameters of the problem, when the applied magnetic field is sufficiently small and the sample sufficiently thick. This optimal scaling law is proven via ansatz-free lower bounds and an explicit branching construction which refines further and further as one approaches the surface of the sample. Two different regimes appear, with different scaling exponents. In the first regime, the branching leads to an almost uniform magnetic field pattern on the boundary; in the second one the inhomogeneity survives up to the boundary.

preprint2015arXiv

Hölder regularity for a non-linear parabolic equation driven by space-time white noise

We consider the non-linear equation $T^{-1} u+\partial_tu-\partial_x^2π(u)=ξ$ driven by space-time white noise $ξ$, which is uniformly parabolic because we assume that $π'$ is bounded away from zero and infinity. Under the further assumption of Lipschitz continuity of $π'$ we show that the stationary solution is - as for the linear case - almost surely Hölder continuous with exponent $α$ for any $α<\frac{1}{2}$ w. r. t. the parabolic metric. More precisely, we show that the corresponding local Hölder norm has stretched exponential moments. On the stochastic side, we use a combination of martingale arguments to get second moment estimates with concentration of measure arguments to upgrade to Gaussian moments. On the deterministic side, we first perform a Campanato iteration based on the De Giorgi-Nash Theorem as well as finite and infinitesimal versions of the $H^{-1}$-contraction principle, which yields Gaussian moments for a weaker Hölder norm. In a second step this estimate is improved to the optimal Hölder exponent at the expense of weakening the integrability to stretched exponential.

preprint2015arXiv

Interpolation inequalities in pattern formation

We prove some interpolation inequalities which arise in the analysis of pattern formation in physics. They are the strong version of some already known estimates in weak form that are used to give a lower bound of the energy in many contexts (coarsening and branching in micromagnetics and superconductors). The main ingredient in the proof of our inequalities is a geometric construction which was first used by Choksi, Conti, Kohn, and one of the authors in the study of branching in superconductors.

preprint2015arXiv

New bounds for the inhomogenous Burgers and the Kuramoto-Sivashinsky equations

We give a substantially simplified proof of near-optimal estimate on the Kuramoto-Sivashinsky equation from [F. Otto, "Optimal bounds on the Kuramoto-Sivashinsky equation", JFA 2009], at the same time slightly improving the result. The result in the above cited paper relied on two ingredients: a regularity estimate for capillary Burgers and an a novel priori estimate for the inhomogeneous inviscid Burgers equation, which works out that in many ways the conservative transport nonlinearity acts as a coercive term. It is the proof of the second ingredient that we substantially simplify by proving a modified Kármán-Howarth-Monin identity for solutions of the inhomogeneous inviscid Burgers equation. This gives a new interpretation of the results obtained in [F. Golse, B. Perthame "Optimal regularizing effect for scalar conservation laws", Rev. Mat. Iber., 2013].

preprint2015arXiv

Quantitative stochastic homogenization: local control of homogenization error through corrector

This note addresses the homogenization error for linear elliptic equations in divergence-form with random stationary coefficients. The homogenization error is measured by comparing the quenched Green's function to the Green's function belonging to the homogenized coefficients, more precisely, by the (relative) spatial decay rate of the difference of their second mixed derivatives. The contribution of this note is purely deterministic: It uses the expanded notion of corrector, namely the couple of scalar and vector potentials $(ϕ,σ)$, and shows that the rate of sublinear growth of $(ϕ,σ)$ at the points of interest translates one-to-one into the decay rate.

preprint2015arXiv

Sublinear growth of the corrector in stochastic homogenization: Optimal stochastic estimates for slowly decaying correlations

We establish sublinear growth of correctors in the context of stochastic homogenization of linear elliptic PDEs. In case of weak decorrelation and "essentially Gaussian" coefficient fields, we obtain optimal (stretched exponential) stochastic moments for the minimal radius above which the corrector is sublinear. Our estimates also capture the quantitative sublinearity of the corrector (caused by the quantitative decorrelation on larger scales) correctly. The result is based on estimates on the Malliavin derivative for certain functionals which are basically averages of the gradient of the corrector, on concentration of measure, and on a mean value property for $a$-harmonic functions.

preprint2014arXiv

Annealed estimates on the Green function

We consider a random, uniformly elliptic coefficient field $a(x)$ on the $d$-dimensional integer lattice $\mathbb{Z}^d$. We are interested in the spatial decay of the quenched elliptic Green function $G(a;x,y)$. Next to stationarity, we assume that the spatial correlation of the coefficient field decays sufficiently fast to the effect that a logarithmic Sobolev inequality holds for the ensemble $\langle\cdot\rangle$. We prove that all stochastic moments of the first and second mixed derivatives of the Green function, that is, $\langle|\nabla_x G(x,y)|^p\rangle$ and $\langle|\nabla_x\nabla_y G(x,y)|^p\rangle$, have the same decay rates in $|x-y|\gg 1$ as for the constant coefficient Green function, respectively. This result relies on and substantially extends the one by Delmotte and Deuschel \cite{DeuschelDelmotte}, which optimally controls second moments for the first derivatives and first moments of the second mixed derivatives of $G$, that is, $\langle|\nabla_x G(x,y)|^2\rangle$ and $\langle|\nabla_x\nabla_y G(x,y)|\rangle$. As an application, we are able to obtain optimal estimates on the random part of the homogenization error even for large ellipticity contrast.

preprint2014arXiv

Corrector estimates for elliptic systems with random periodic coefficients

We consider an elliptic system of equations on the torus $\left[ -\frac{L}{2}, \frac{L}{2} \right)^d$ with random coefficients $A$, that are assumed to be coercive and stationary. Using two different approaches we obtain moment bounds on the gradient of the corrector, independent of the domain size $L$. In the first approach we use Green function representation. For that we require $A$ to be locally Hölder continuous and distribution of $A$ to satisfy Logarithmic Sobolev inequality. The second method works for non-smooth (possibly discontinuous) coefficients, and it requires that statistics of $A$ satisfies Spectral Gap estimate.

preprint2014arXiv

Moment bounds for the corrector in stochastic homogenization of a percolation model

We study the corrector equation in stochastic homogenization for a simplified Bernoulli percolation model on $\mathbb{Z}^d$, $d>2$. The model is obtained from the classical $\{0,1\}$-Bernoulli bond percolation by conditioning all bonds parallel to the first coordinate direction to be open. As a main result we prove (in fact for a slightly more general model) that stationary correctors exist and that all finite moments of the corrector are bounded. This extends a previous result in [GO1], where uniformly elliptic conductances are treated, to the degenerate case. With regard to the associated random conductance model, we obtain as a side result that the corrector not only grows sublinearly, but slower than any polynomial rate. Our argument combines a quantification of ergodicity by means of a Spectral Gap on Glauber dynamics with regularity estimates on the gradient of the elliptic Green's function.

preprint2014arXiv

On annealed elliptic Green function estimates

We consider a random, uniformly elliptic coefficient field $a$ on the lattice $\mathbb{Z}^d$. The distribution $\langle \cdot \rangle$ of the coefficient field is assumed to be stationary. Delmotte and Deuschel showed that the gradient and second mixed derivative of the parabolic Green function $G(t,x,y)$ satisfy optimal annealed estimates which are $L^2$ resp. $L^1$ in probability, i.e. they obtained bounds on $\langle |\nabla_x G(t,x,y)|^2 \rangle^{\frac{1}{2}}$ and $\langle |\nabla_x \nabla_y G(t,x,y)| \rangle$, see T. Delmotte and J.-D. Deuschel: On estimating the derivatives of symmetric diffusions in stationary random environments, with applications to the $\nablaϕ$ interface model, Probab. Theory Relat. Fields 133 (2005), 358--390. In particular, the elliptic Green function $G(x,y)$ satisfies optimal annealed bounds. In a recent work, the authors extended these elliptic bounds to higher moments, i.e. $L^p$ in probability for all $p<\infty$, see D. Marahrens and F. Otto: {Annealed estimates on the Green function}, arXiv:1304.4408 (2013). In this note, we present a new argument that relies purely on elliptic theory to derive the elliptic estimates (see Proposition 1.2 below) for $\langle |\nabla_x G(x,y)|^2 \rangle^{\frac{1}{2}}$ and $\langle |\nabla_x \nabla_y G(x,y)| \rangle$.

preprint2014arXiv

Quantitative estimates on the periodic approximation of the corrector in stochastic homogenization

In the present contribution we establish quantitative results on the periodic approximation of the corrector equation for the stochastic homogenization of linear elliptic equations in divergence form, when the diffusion coefficients satisfy a spectral gap estimate in probability, and for $d>2$. The main difference with respect to the first part of [Gloria-Otto, arXiv:1409.0801] is that we avoid here the use of Green's functions and more directly rely on the De Giorgi-Nash-Moser theory.

preprint2014arXiv

Quantitative results on the corrector equation in stochastic homogenization

We derive optimal estimates in stochastic homogenization of linear elliptic equations in divergence form in dimensions $d\ge 2$. In previous works we studied the model problem of a discrete elliptic equation on $\mathbb{Z}^d$. Under the assumption that a spectral gap estimate holds in probability, we proved that there exists a stationary corrector field in dimensions $d>2$ and that the energy density of that corrector behaves as if it had finite range of correlation in terms of the variance of spatial averages - the latter decays at the rate of the central limit theorem. In this article we extend these results, and several other estimates, to the case of a continuum linear elliptic equation whose (not necessarily symmetric) coefficient field satisfies a continuum version of the spectral gap estimate. In particular, our results cover the example of Poisson random inclusions.

preprint2014arXiv

Upper bounds on Nusselt number at finite Prandtl number

We study Rayleigh Bénard convection based on the Boussinesq approximation. We are interested in upper bounds on the Nusselt number $\mathrm{Nu}$, the upwards heat transport, in terms of the Rayleigh number $\mathrm{Ra}$, that characterizes the relative strength of the driving mechanism and the Prandtl number $\mathrm{Pr}$, that characterizes the strength of the inertial effects. We show that, up to logarithmic corrections, the upper bound $\mathrm{Nu}\lesssim \mathrm{Ra}^{\frac{1}{3}}$ of Constantin and Doering in 1999 persists as long as $\mathrm{Pr}\gtrsim \mathrm{Ra}^{\frac{1}{3}}$ and then crosses over to $\mathrm{Nu}\lesssim\mathrm{Pr}^{-\frac{1}{2}}\mathrm{Ra}^{\frac{1}{2}}$. This result improves the one of Wang by going beyond the perturbative regime $\mathrm{Pr} \gg \mathrm{Ra}$. The proof uses a new way to estimate the transport nonlinearity in the Navier-Stokes equations capitalizing on the no-slip boundary condition. It relies on a new Calderón-Zygmund estimate for the non-stationary Stokes equations in $L^1$ with a borderline Muckenhoupt weight.

preprint2013arXiv

A reduced model for domain walls in soft ferromagnetic films at the cross-over from symmetric to asymmetric wall types

We study the Landau-Lifshitz model for the energy of multi-scale transition layers -- called "domain walls" -- in soft ferromagnetic films. Domain walls separate domains of constant magnetization vectors $m^\pm \in \mathbb{S}^2$ that differ by an angle $2α$. Assuming translation invariance tangential to the wall, our main result is the rigorous derivation of a reduced model for the energy of the optimal transition layer, which in a certain parameter regime confirms the experimental, numerical and physical predictions: The minimal energy splits into a contribution from an asymmetric, divergence-free core which performs a partial rotation in $\mathbb{S}^2$ by an angle $2θ$, and a contribution from two symmetric, logarithmically decaying tails, each of which completes the rotation from angle $θ$ to $α$ in $\mathbb{S}^1$. The angle $θ$ is chosen such that the total energy is minimal. The contribution from the symmetric tails is known explicitly, while the contribution from the asymmetric core is analyzed in [7]. Our reduced model is the starting point for the analysis of a bifurcation phenomenon from symmetric to asymmetric domain walls. Moreover, it allows for capturing asymmetric domain walls including their extended tails (which were previously inaccessible to brute-force numerical simulation).

preprint2013arXiv

Invariant measure of the stochastic Allen-Cahn equation: the regime of small noise and large system size

We study the invariant measure of the one-dimensional stochastic Allen-Cahn equation for a small noise strength and a large but finite system. We endow the system with inhomogeneous Dirichlet boundary conditions that enforce at least one transition from -1 to 1. (Our methods can be applied to other boundary conditions as well.) We are interested in the competition between the energy that should be minimized due to the small noise strength and the entropy that is induced by the large system size. Our methods handle system sizes that are exponential with respect to the inverse noise strength, up to the critical exponential size predicted by the heuristics. We capture the competition between energy and entropy through upper and lower bounds on the probability of extra transitions between -1 and 1. These bounds are sharp on the exponential scale and imply in particular that the probability of having one and only one transition from -1 to +1 is exponentially close to one. In addition, we show that the position of the transition layer is uniformly distributed over the system on scales larger than the logarithm of the inverse noise strength. Our arguments rely on local large deviation bounds, the strong Markov property, the symmetry of the potential, and measure-preserving reflections.

preprint2013arXiv

Uniform logarithmic Sobolev inequalities for conservative spin systems with super-quadratic single-site potential

We consider a noninteracting unbounded spin system with conservation of the mean spin. We derive a uniform logarithmic Sobolev inequality (LSI) provided the single-site potential is a bounded perturbation of a strictly convex function. The scaling of the LSI constant is optimal in the system size. The argument adapts the two-scale approach of Grunewald, Villani, Westdickenberg and the second author from the quadratic to the general case. Using an asymmetric Brascamp-Lieb-type inequality for covariances, we reduce the task of deriving a uniform LSI to the convexification of the coarse-grained Hamiltonian, which follows from a general local Cramér theorem.

preprint2012arXiv

An optimal error estimate in stochastic homogenization of discrete elliptic equations

This paper is the companion article to [Ann. Probab. 39 (2011) 779--856]. We consider a discrete elliptic equation on the $d$-dimensional lattice $\mathbb{Z}^d$ with random coefficients $A$ of the simplest type: They are identically distributed and independent from edge to edge. On scales large w.r.t. the lattice spacing (i.e., unity), the solution operator is known to behave like the solution operator of a (continuous) elliptic equation with constant deterministic coefficients. This symmetric "homogenized" matrix $A_{\mathrm{hom}}=a_{\mathrm{hom}}\mathrm{Id}$ is characterized by $ξ\cdot A_{\mathrm{hom}}ξ=<(ξ+\nablaϕ)\cdot A(ξ+\nablaϕ)>$ for any direction $ξ\in\mathbb{R}^d$, where the random field $ϕ$ (the "corrector") is the unique solution of $-\nabla^*\cdot A(ξ+\nablaϕ)=0$ in $\mathbb{Z}^d$ such that $ϕ(0)=0$, $\nablaϕ$ is stationary and $<\nablaϕ>=0$, $<\cdot>$ denoting the ensemble average (or expectation).

preprint2011arXiv

An optimal variance estimate in stochastic homogenization of discrete elliptic equations

We consider a discrete elliptic equation on the $d$-dimensional lattice $\mathbb{Z}^d$ with random coefficients $A$ of the simplest type: they are identically distributed and independent from edge to edge. On scales large w.r.t. the lattice spacing (i.e., unity), the solution operator is known to behave like the solution operator of a (continuous) elliptic equation with constant deterministic coefficients. This symmetric ``homogenized'' matrix $A_{\mathrm {hom}}=a_{\mathrm {hom}}\operatorname {Id}$ is characterized by $ξ\cdot A_{\mathrm {hom}}ξ=\langle(ξ+\nablaϕ)\cdot A(ξ+\nablaϕ)\rangle$ for any direction $ξ\in\mathbb {R}^d$, where the random field $ϕ$ (the ``corrector'') is the unique solution of $-\nabla^*\cdot A(ξ+\nablaϕ)=0$ such that $ϕ(0)=0$, $\nablaϕ$ is stationary and $\langle\nablaϕ\rangle=0$, $\langle\cdot\rangle$ denoting the ensemble average (or expectation). It is known (``by ergodicity'') that the above ensemble average of the energy density $\mathcal {E}=(ξ+\nablaϕ)\cdot A(ξ+\nablaϕ)$, which is a stationary random field, can be recovered by a system average. We quantify this by proving that the variance of a spatial average of $\mathcal {E}$ on length scales $L$ satisfies the optimal estimate, that is, $\operatorname {var}[\sum \mathcal {E}η_L]\lesssim L^{-d}$, where the averaging function [i.e., $\sumη_L=1$, $\operatorname {supp}(η_L)\subset\{|x|\le L\}$] has to be smooth in the sense that $|\nablaη_L|\lesssim L^{-1-d}$. In two space dimensions (i.e., $d=2$), there is a logarithmic correction. This estimate is optimal since it shows that smooth averages of the energy density $\mathcal {E}$ decay in $L$ as if $\mathcal {E}$ would be independent from edge to edge (which it is not for $d>1$). This result is of practical significance, since it allows to estimate the dominant error when numerically computing $a_{\mathrm {hom}}$.

preprint2011arXiv

The Formation and Coarsening of the Concertina Pattern

The concertina is a magnetization pattern in elongated thin-film elements of a soft material. It is a ubiquitous domain pattern that occurs in the process of magnetization reversal in direction of the long axis of the small element. Van den Berg argued that this pattern grows out of the flux closure domains as the external field is reduced. Based on experimental observations and theory, we argue that in sufficiently elongated thin-film elements, the concertina pattern rather bifurcates from an oscillatory buckling mode. Using a reduced model derived by asymptotic analysis and investigated by numerical simulation, we quantitatively predict the average period of the concertina pattern and qualitatively predict its hysteresis. In particular, we argue that the experimentally observed coarsening of the concertina pattern is due to secondary bifurcations related to an Eckhaus instability. We also link the concertina pattern to the magnetization ripple and discuss the effect of a weak (crystalline or induced) anisotropy.

preprint2010arXiv

A compactness result for Landau state in thin-film micromagnetics

We deal with a nonconvex and nonlocal variational problem coming from thin-film micromagnetics. It consists in a free-energy functional depending on two small parameters $\eps$ and $η$ and defined over $S^2-$vector fields $m$ that are tangent at the boundary of a two-dimensional domain $Ω$. We are interested in the behavior of minimizers as $\eps, η\to 0$. The minimizers tend to be in-plane away from a region of length scale $\eps$ (generically, an interior vortex ball or two boundary vortex balls) and of vanishing divergence, so that $S^1-$transition layers of length scale $η$ (Néel walls) are enforced by the boundary condition. We first prove an upper bound for the minimal energy that corresponds to the cost of a vortex and the configuration of Néel walls associated to the viscosity solution, so-called Landau state. Our main result concerns the compactness of vector fields $m_{\eps, η}$ of energies close to the Landau state in the regime where a vortex is energetically more expensive than a Néel wall. Our method uses techniques developed for the Ginzburg-Landau type problems for the concentration of energy on vortex balls, together with an approximation argument of $S^2-$vector fields by $S^1-$vector fields away from the vortex balls.