-
Boundary Layer Estimates in Stochastic Homogenization
Authors:
Peter Bella,
Julian Fischer,
Marc Josien,
Claudia Raithel
Abstract:
We prove quantitative decay estimates for the boundary layer corrector in stochastic homogenization in the case of a half-space boundary. Our estimates are of optimal order and show that the gradient of the boundary layer corrector features nearly fluctuation-order decay; its expected value decays even one order faster. As a corollary, we deduce estimates on the accuracy of the representative volu…
▽ More
We prove quantitative decay estimates for the boundary layer corrector in stochastic homogenization in the case of a half-space boundary. Our estimates are of optimal order and show that the gradient of the boundary layer corrector features nearly fluctuation-order decay; its expected value decays even one order faster. As a corollary, we deduce estimates on the accuracy of the representative volume element method for the computation of effective coefficients: our understanding of the decay of boundary layers enables us to improve the order of convergence of the RVE method for $d\geq 3$.
△ Less
Submitted 19 March, 2024;
originally announced March 2024.
-
Shock-driven nucleation and self-organization of dislocations in the dynamical Peierls model
Authors:
Yves-Patrick Pellegrini,
Marc Josien
Abstract:
Dynamic nucleation of dislocations caused by a stress front ('shock') of amplitude $σ_{\rm a}$ moving with speed $V$ is investigated by solving numerically the Dynamic Peierls Equation with an efficient method. Speed $V$ and amplitude $σ_{\rm a}$ are considered as independent variables, with $V$ possibly exceeding the longitudinal wavespeed $c_{\rm L}$. Various reactions between dislocations take…
▽ More
Dynamic nucleation of dislocations caused by a stress front ('shock') of amplitude $σ_{\rm a}$ moving with speed $V$ is investigated by solving numerically the Dynamic Peierls Equation with an efficient method. Speed $V$ and amplitude $σ_{\rm a}$ are considered as independent variables, with $V$ possibly exceeding the longitudinal wavespeed $c_{\rm L}$. Various reactions between dislocations take place such as scattering, dislocation-pair nucleation, annihilation, and crossing. Pairs of edge dislocation are always nucleated with speed $v\gtrsim c_{\rm L}$ (and likewise for screws with $c_{\rm L}$ replaced by $c_{\rm S}$, the shear wavespeed). The plastic wave exhibits self-organization, forming distinct `bulk' and `front' zones. Nucleations occur either within the bulk or at the zone interface, depending on the value of $V$. The front zone accumulates dislocations that are expelled from the bulk or from the interface. In each zone, dislocation speeds and densities are measured as functions of simulation parameters. The densities exhibit a scaling behavior with stress, given by $((σ_a/σ_{\rm th})^2-1)^β$, where $σ_{\rm th}$ represents the nucleation threshold and $0<β<1$.
△ Less
Submitted 8 August, 2023;
originally announced August 2023.
-
Bias in the representative volume element method: periodize the ensemble instead of its realizations
Authors:
Nicolas Clozeau,
Marc Josien,
Felix Otto,
Qiang Xu
Abstract:
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 metho…
▽ More
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.
△ Less
Submitted 27 May, 2022;
originally announced May 2022.
-
Stochastic homogenization and geometric singularities : a study on corners
Authors:
Marc Josien,
Claudia Raithel,
Mathias Schäffner
Abstract:
In this contribution we are interested in the quantitative homogenization properties of linear elliptic equations with homogeneous Dirichlet boundary data in polygonal domains with corners. To begin our study of this situation, we consider the setting of an angular sector in 2 dimensions : Unlike in the whole-space, on such a sector there exist non-smooth harmonic functions (these depend on the an…
▽ More
In this contribution we are interested in the quantitative homogenization properties of linear elliptic equations with homogeneous Dirichlet boundary data in polygonal domains with corners. To begin our study of this situation, we consider the setting of an angular sector in 2 dimensions : Unlike in the whole-space, on such a sector there exist non-smooth harmonic functions (these depend on the angle of the sector). Here, we construct extended homogenization correctors corresponding to these harmonic functions and prove growth estimates for these which are quasi-optimal, namely optimal up to a logarithmic loss. Our construction of the corner correctors relies on a large-scale regularity theory for a-harmonic functions in the sector, which we also prove and which, as a by-product, yields a Liouville principle. We also propose a nonstandard 2-scale expansion, which is adapted to the sectoral domain and incorporates the corner correctors. Our final result is a quasi-optimal error estimate for this adapted 2-scale expansion.
△ Less
Submitted 24 January, 2022;
originally announced January 2022.
-
Representative volume elements for matrix-inclusion composites -- a computational study on periodizing the ensemble
Authors:
Matti Schneider,
Marc Josien,
Felix Otto
Abstract:
We investigate volume-element sampling strategies for the stochastic homogenization of particle-reinforced composites and show, via computational experiments, that an improper treatment of particles intersecting the boundary of the computational cell may affect the accuracy of the computed effective properties. Motivated by recent results on a superior convergence rate of the systematic error for…
▽ More
We investigate volume-element sampling strategies for the stochastic homogenization of particle-reinforced composites and show, via computational experiments, that an improper treatment of particles intersecting the boundary of the computational cell may affect the accuracy of the computed effective properties. Motivated by recent results on a superior convergence rate of the systematic error for periodized ensembles compared to taking snapshots of ensembles, we conduct computational experiments for microstructures with circular, spherical and cylindrical inclusions and monitor the systematic errors in the effective thermal conductivity for snapshots of ensembles compared to working with microstructures sampled from periodized ensembles. We observe that the standard deviation of the apparent properties computed on microstructures sampled from the periodized ensembles decays at the scaling expected from the central limit theorem. In contrast, the standard deviation for the snapshot ensembles shows an inferior decay rate at high filler content. The latter effect is caused by additional long-range correlations that necessarily appear in particle-reinforced composites at high, industrially relevant, volume fractions. Periodized ensembles, however, appear to be less affected by these correlations. Our findings provide guidelines for working with digital volume images of material microstructures and the design of representative volume elements for computational homogenization.
△ Less
Submitted 13 March, 2021;
originally announced March 2021.
-
The annealed Calderon-Zygmund estimate as convenient tool in quantitative stochastic homogenization
Authors:
Marc Josien,
Felix Otto
Abstract:
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 optim…
▽ More
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.
△ Less
Submitted 18 May, 2020;
originally announced May 2020.
-
Quantitative homogenization for the case of an interface between two heterogeneous media
Authors:
Marc Josien,
Claudia Raithel
Abstract:
In this article we are interested in quantitative homogenization results for linear elliptic equations in the non-stationary situation of a straight interface between two heterogenous media. This extends the previous work [Josien, 2019] to a substantially more general setting, in which the surrounding heterogeneous media may be periodic or random stationary and ergodic. Our main result is a quanti…
▽ More
In this article we are interested in quantitative homogenization results for linear elliptic equations in the non-stationary situation of a straight interface between two heterogenous media. This extends the previous work [Josien, 2019] to a substantially more general setting, in which the surrounding heterogeneous media may be periodic or random stationary and ergodic. Our main result is a quantification of the sublinearity of a homogenization corrector adapted to the interface, which we construct using an improved version of the method developed in [Fischer and Raithel, 2017]. This quantification is optimal up to a logarithmic loss and allows to derive almost-optimal convergence rates.
△ Less
Submitted 2 December, 2019;
originally announced December 2019.
-
Local precised approximation in multiscale problems with local defects
Authors:
Xavier Blanc,
Marc Josien,
Claude Le Bris
Abstract:
We proceed here with our systematic study, initiated in [3], of multiscale problems with defects, within the context of homogenization theory. The case under consideration here is that of a diffusion equation with a diffusion coefficient of the form of a periodic function perturbed by an $L^r (R^d ) , 1 < r < +$\infty$…
▽ More
We proceed here with our systematic study, initiated in [3], of multiscale problems with defects, within the context of homogenization theory. The case under consideration here is that of a diffusion equation with a diffusion coefficient of the form of a periodic function perturbed by an $L^r (R^d ) , 1 < r < +$\infty$$ , function modeling a localized defect. We outline the proof of the following approximation result: the corrector function, the existence of which has been established in [3,4], allows to approximate the solution of the original multiscale equation with essentially the same accuracy as in the purely periodic case. The rates of convergence may however vary, and are made precise, depending upon the $L^r$ integrability of the defect. The generalization to an abstract setting is mentioned. Our proof exactly follows, step by step, the pattern of the original proof of Avellaneda and Lin in [1] in the periodic case, extended in the works of Kenig and collaborators [13], and borrows a lot from it. The details of the results announced in this Note are given in our forthcoming publications [2,12].
△ Less
Submitted 18 January, 2019;
originally announced January 2019.
-
Precised approximations in elliptic homogenization beyond the periodic setting
Authors:
Xavier Blanc,
Marc Josien,
Claude Le Bris
Abstract:
We consider homogenization problems for linear elliptic equations in divergence form. The coecients are assumed to be a local perturbation of some periodic background. We prove $W^{1,p}$ and Lipschitz convergence of the two-scale expansion, with explicit rates. For this purpose, we use a corrector adapted to this particular setting, and dened in [10, 11], and apply the same strategy of proof as Av…
▽ More
We consider homogenization problems for linear elliptic equations in divergence form. The coecients are assumed to be a local perturbation of some periodic background. We prove $W^{1,p}$ and Lipschitz convergence of the two-scale expansion, with explicit rates. For this purpose, we use a corrector adapted to this particular setting, and dened in [10, 11], and apply the same strategy of proof as Avellaneda and Lin in [1]. We also propose an abstract setting generalizing our particular assumptions for which the same estimates hold.
△ Less
Submitted 18 December, 2018;
originally announced December 2018.
-
Decomposition and pointwise estimates of periodic Green functions of some elliptic equations with periodic oscillatory coefficients
Authors:
Marc Josien
Abstract:
This article is about the $\mathbb{Z}^d$-periodic Green function $G_n(x,y)$ of the multiscale elliptic operator $Lu=-{\rm div}\left( A(n\cdot) \cdot \nabla u \right)$, where $A(x)$ is a $\mathbb{Z}^d$-periodic, coercive, and Hölder continuous matrix, and $n$ is a large integer. We prove here pointwise estimates on $G_n(x,y)$, $\nabla_x G_n(x,y)$, $\nabla_y G_n(x,y)$ and…
▽ More
This article is about the $\mathbb{Z}^d$-periodic Green function $G_n(x,y)$ of the multiscale elliptic operator $Lu=-{\rm div}\left( A(n\cdot) \cdot \nabla u \right)$, where $A(x)$ is a $\mathbb{Z}^d$-periodic, coercive, and Hölder continuous matrix, and $n$ is a large integer. We prove here pointwise estimates on $G_n(x,y)$, $\nabla_x G_n(x,y)$, $\nabla_y G_n(x,y)$ and $\nabla_x \nabla_y G_n(x,y)$ in dimensions $d \geq 2$. Moreover, we derive an explicit decomposition of this Green function, which is of independent interest. These results also apply for systems.
△ Less
Submitted 24 July, 2018;
originally announced July 2018.
-
Mathematical properties of the Weertman equation
Authors:
Marc Josien
Abstract:
We derive here some mathematical properties of the Weertman equation and show it is the limit of an evolution equation. The Weertman equation is a semilinear integrodifferential equation involving a fractional Laplacian. In addition to this purely theoretical interest, the results proven here give a solid ground to a numerical approach that we have implemented (see arXiv:1704.04489).
We derive here some mathematical properties of the Weertman equation and show it is the limit of an evolution equation. The Weertman equation is a semilinear integrodifferential equation involving a fractional Laplacian. In addition to this purely theoretical interest, the results proven here give a solid ground to a numerical approach that we have implemented (see arXiv:1704.04489).
△ Less
Submitted 18 September, 2017;
originally announced September 2017.
-
Fourier-based numerical approximation of the Weertman equation for moving dislocations
Authors:
Marc Josien,
Yves-Patrick Pellegrini,
Frédéric Legoll,
Claude Le Bris
Abstract:
This work discusses the numerical approximation of a nonlinear reaction-advection-diffusion equation, which is a dimensionless form of the Weertman equation. This equation models steadily-moving dislocations in materials science. It reduces to the celebrated Peierls-Nabarro equation when its advection term is set to zero. The approach rests on considering a time-dependent formulation, which admits…
▽ More
This work discusses the numerical approximation of a nonlinear reaction-advection-diffusion equation, which is a dimensionless form of the Weertman equation. This equation models steadily-moving dislocations in materials science. It reduces to the celebrated Peierls-Nabarro equation when its advection term is set to zero. The approach rests on considering a time-dependent formulation, which admits the equation under study as its long-time limit. Introducing a Preconditioned Collocation Scheme based on Fourier transforms, the iterative numerical method presented solves the time-dependent problem, delivering at convergence the desired numerical solution to the Weertman equation. Although it rests on an explicit time-evolution scheme, the method allows for large time steps, and captures the solution in a robust manner. Numerical results illustrate the efficiency of the approach for several types of nonlinearities.
△ Less
Submitted 14 April, 2017;
originally announced April 2017.
-
From the Newton equation to the wave equation : the case of shock waves
Authors:
Xavier Blanc,
Marc Josien
Abstract:
We study the macroscopic limit of a chain of atoms governed by the Newton equation. It is known from the work of Blanc, Le Bris, Lions, that this limit is the solution of a nonlinear wave equation, as long as this solution remains smooth. We show, numerically and mathematically that, if the distances between particles remain bounded, it is not the case any more when there are shocks -at least for…
▽ More
We study the macroscopic limit of a chain of atoms governed by the Newton equation. It is known from the work of Blanc, Le Bris, Lions, that this limit is the solution of a nonlinear wave equation, as long as this solution remains smooth. We show, numerically and mathematically that, if the distances between particles remain bounded, it is not the case any more when there are shocks -at least for a convex nearest-neighbour interaction potential with convex derivative.
△ Less
Submitted 3 May, 2016;
originally announced May 2016.
-
New bounds for the inhomogenous Burgers and the Kuramoto-Sivashinsky equations
Authors:
Michael Goldman,
Marc Josien,
Felix Otto
Abstract:
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…
▽ More
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].
△ Less
Submitted 20 March, 2015;
originally announced March 2015.