Processing math: 60%
Research article Special Issues

Homogenisation of high-contrast brittle materials

  • This paper is an overview on some recent results concerning the variational analysis of static fracture in the so-called high-contrast brittle composite materials. The paper is divided into two main parts. The first part is devoted to establish a compactness result for a general class of free-discontinuity functionals with degenerate (or high-contrast) integrands. The second part is focussed on some specific examples which show that the degeneracy of the integrands may lead to non-standard limit effects, which are specific to this high-contrast setting.

    Citation: Caterina Ida Zeppieri. Homogenisation of high-contrast brittle materials[J]. Mathematics in Engineering, 2020, 2(1): 174-202. doi: 10.3934/mine.2020009

    Related Papers:

    [1] Manuel Friedrich . Griffith energies as small strain limit of nonlinear models for nonsimple brittle materials. Mathematics in Engineering, 2020, 2(1): 75-100. doi: 10.3934/mine.2020005
    [2] Stefano Almi, Giuliano Lazzaroni, Ilaria Lucardesi . Crack growth by vanishing viscosity in planar elasticity. Mathematics in Engineering, 2020, 2(1): 141-173. doi: 10.3934/mine.2020008
    [3] Giacomo Canevari, Arghir Zarnescu . Polydispersity and surface energy strength in nematic colloids. Mathematics in Engineering, 2020, 2(2): 290-312. doi: 10.3934/mine.2020015
    [4] Zhiwen Zhao . Exact solutions for the insulated and perfect conductivity problems with concentric balls. Mathematics in Engineering, 2023, 5(3): 1-11. doi: 10.3934/mine.2023060
    [5] Paola F. Antonietti, Chiara Facciolà, Marco Verani . Unified analysis of discontinuous Galerkin approximations of flows in fractured porous media on polygonal and polyhedral grids. Mathematics in Engineering, 2020, 2(2): 340-385. doi: 10.3934/mine.2020017
    [6] Raimondo Penta, Ariel Ramírez-Torres, José Merodio, Reinaldo Rodríguez-Ramos . Effective governing equations for heterogenous porous media subject to inhomogeneous body forces. Mathematics in Engineering, 2021, 3(4): 1-17. doi: 10.3934/mine.2021033
    [7] Patrick Dondl, Martin Jesenko, Martin Kružík, Jan Valdman . Linearization and computation for large-strain visco-elasticity. Mathematics in Engineering, 2023, 5(2): 1-15. doi: 10.3934/mine.2023030
    [8] Federico Cluni, Vittorio Gusella, Dimitri Mugnai, Edoardo Proietti Lippi, Patrizia Pucci . A mixed operator approach to peridynamics. Mathematics in Engineering, 2023, 5(5): 1-22. doi: 10.3934/mine.2023082
    [9] Chiara Gavioli, Pavel Krejčí . Deformable porous media with degenerate hysteresis in gravity field. Mathematics in Engineering, 2025, 7(1): 35-60. doi: 10.3934/mine.2025003
    [10] Luca Formaggia, Alessio Fumagalli, Anna Scotti . A multi-layer reactive transport model for fractured porous media. Mathematics in Engineering, 2022, 4(1): 1-32. doi: 10.3934/mine.2022008
  • This paper is an overview on some recent results concerning the variational analysis of static fracture in the so-called high-contrast brittle composite materials. The paper is divided into two main parts. The first part is devoted to establish a compactness result for a general class of free-discontinuity functionals with degenerate (or high-contrast) integrands. The second part is focussed on some specific examples which show that the degeneracy of the integrands may lead to non-standard limit effects, which are specific to this high-contrast setting.


    In this note we analyse the large-scale behaviour of high-contrast composite materials which can undergo fracture. In a variational setting, the microscopic behaviour of high-contrast composites is typically described by means of scale-dependent energy functionals with "degenerate" integrands. For brittle materials the scale-dependent energies are of the general form

    Fε(u)=Ωfε(x,u)dx+Sugε(x,νu)dHn1, (1.1)

    where ε>0 describes both the composite-microstructure and the degeneracy of the mechanical properties of the material (cf. (1.3)). In (1.1) the variable u:ΩRnR belongs to SBV(Ω), the space of special functions of bounded variation in Ω. In this simplified scalar setting, u represents an anti-plane displacement and Ω is the cross-section of an infinite cylindrical body. Being u an SBV-function, discontinuities are allowed and the discontinuity set of u, denoted by Su, models the cracks in the material. The deformation gradient Du can be decomposed into the sum of a bulk part udx and a surface part (u+u)νuHn1Su, where u denotes the approximate gradient of u, u+ and u the traces of u on both sides of Su, and νu denotes the (generalised) normal to Su. The volume term in Fε represents the elastic energy stored in the unfractured part of the material, whereas the surface term in Fε accounts for the presence of cracks. According to the Griffith criterion, in brittle materials, already for the smallest crack-amplitude, there is no interaction between the two lips of the crack, so that the corresponding fracture energy does not depend on [u]=u+u.

    For finite-contrast brittle materials, the limit behaviour of energies of type (1.1) is by-now well-understood and the corresponding theory provides a rigorous micro-to-macro upscaling for brittle fracture. In fact, if fε and gε satisfy (mild regularity assumptions and) standard growth and coercivity conditions of type

    c1|ξ|pfε(x,ξ)c2(1+|ξ|p)andc3gε(x,ν)c4, (1.2)

    for every ε>0, x,ξRn, νSn1, for some p>1, and 0<c1c2<+, 0<c3c4<+, then in [18] Giacomini and Ponsiglione showed, among other, that the limit behaviour of Fε is captured by a scale-independent free-discontinuity functional of the same type as Fε; i.e.,

    F0(u)=Ωf0(x,u)dx+Sug0(x,νu)dHn1,

    with f0 and g0 also satisfying (1.2). Under these assumptions, Giacomini and Ponsiglione also showed that volume and surface energy decouple in the limit, so that the energy density f0 is not affected by the presence of the surface term in Fε, whereas the surface energy density g0 is not affected by the volume term in Fε. In a recent work, Cagnetti, Dal Maso, Scardia and Zeppieri [12] generalised the asymptotic analysis carried out in [18] and devised (nearly optimal) sufficient conditions which ensure a macroscopic bulk-surface energy decoupling for a wide class of finite-contrast vectorial free-discontinuity functionals which may also depend on [u]. The class of periodic free-discontinuity functionals originally analysed by Braides, Defranceschi and Vitali [9] satisfy the sufficient conditions provided in [12]. Moreover, random free-discontinuity functionals with stationary finite-contrast integrands can be also seen as a special instance of those treated in [12], as shown by Cagnetti, Dal Maso, Scardia and Zeppieri in [13]. Therefore, a volume-surface interaction can be ruled out for a large class of finite-contrast free-discontinuity functionals. In this setting, in particular, microscopic brittle energies always converge to macroscopic brittle energies. However, the general theory established in [9,12,13,18] is not well-suited for studying the large-scale behaviour of those brittle composites whose different constituents have very different mechanical properties from one another. Indeed, in this case the integrands fε and gε in (1.1) may exhibit a so-called high-contrast behaviour and satisfy (1.2) only in a subset Ωε of Ω.

    In the last decade there has been an ever increasing interest in the study of high-contrast free-discontinuity functionals and in the derivation of their effective properties. In particular, the case where (at least) one of the conditions in (1.2) is violated in "many small" periodically distributed regions inside Ω has been considered (see, e.g., [3,4,5,6,11,14,16,17,19,20,21]). Depending on the type of degeneracy of fε and gε, nonstandard limit effects have been also observed. These nonstandard effects are typical of the high-contrast setting and arise from a nontrivial volume-surface limit interaction, which cannot be excluded in this degenerate setting. In fact, in the two companion papers [4,16], Barchiesi, Dal Maso and Zeppieri show that when only gε is degenerate, already for very simple free-discontinuity functionals of Mumford-Shah type, a bulk-surface interaction cannot be ruled out. Namely, a volume-surface coupling can be observed when homogenising a material made of "many" purely brittle inclusions periodically distributed in a connected unbreakable structure, whose fracture-resistance is assumed to be infinite. This coupling produces a homogeneous material whose overall behaviour is of ductile (or cohesive) type; in other words, the homogenised surface energy explicitly depends on [u]. A similar phenomenon is also observed by Barchiesi, Lazzaroni and Zeppieri [6] who show that a ductile behaviour can be seen as the macroscopic effect of a nontrivial volume-surface interaction in the homogenisation of two purely brittle materials with a high-contrast bulk energy. Moreover, in the recent work [19] Pellet, Scardia and Zeppieri prove, instead, that nonstandard constitutive laws may arise when homogenising two purely brittle materials with a high-contrast surface energy. The functionals analysed in [6] and [19] are both of type

    Fαε,βεε(u)=Ωaε(xε)|u|2dx+Subε(xε)dHn1 (1.3)

    where the elastic modulus aε and the fracture resistance (or fracture thoughness) bε are Q-periodic functions and in the unit periodicity cell Q:=(1/2,1/2)n are defined as

    aε(y)={αεify¯Qr1ifyQ¯Qrbε(y)={βεify¯Qr1ifyQ¯Qr

    with αε,βε[0,1], r(0,1), and Qr:=(r/2,r/2)n. Since αε,βε are not bounded away from zero, the functions aε and bε can be degenerate. In their turn, the integrands fε(y,ξ)=aε(y)|ξ|2 and gε(y,ν)=bε(y) in (1.3) will not satisfy, in general, the coercivity conditions in (1.2).

    The limit case αε=βε=0 corresponds to the case of periodically perforated brittle materials studied by Cagnetti and Scardia [14] and by Focardi, Gelli, and Ponsiglione [17] (see also Barchiesi and Focardi [5] for more general free-discontinuity functionals). In spite of the strong degeneracy of the coefficients aε and bε, which in this case are equal to zero in a "large" portion of Ω, in this case it can be proven that the functionals F0,0ε exhibit a limit behaviour which is qualitatively similar to that of free-discontinuity functionals with coercive integrands. Namely, in this case bulk and surface terms do not interact in the limit.

    The aim of this note is to show that, contrary to the coercive case, where general homogenisation results can be proven to describe the limit behaviour of a large class of free-discontinuity functionals, in the non-coercive setting, already for special functionals of type (1.3), a unified homogenisation theory cannot be established. In fact, the limit behaviour of Fαε,βεε is highly sensitive both to the choice of the parameters αε and βε and to their vanishing rate compared to the period of the microstructure ε.

    This note is divided into two main parts and organised as follows. In first part we will deal with sequences of general free-discontinuity functionals of type (1.1) whose coefficients fε and gε are "weakly coercive" or "degenerate"; i.e., they satisfy the lower bounds in (1.2) only in a set Ωε which is obtained removing from Ω many small periodically distributed connected regions. We will use the localisation method of Γ-convergence [8,15] to prove that these kind of functionals are (pre)compact. That is, up to subsequences, they always Γ-converge to a free-discontinuity functional of type

    Ωf(x,u)dx+Sug(x,[u],νu)dHn1.

    Moreover, the limit integrands f and g are non-degenerate and satisfy coercivity conditions of type (1.2) for some positive constants ˆc1,ˆc3 which are strictly smaller than c1,c3, respectively. In this part of the analysis a pivotal role is played by an extension result for SBV-functions defined in periodically perforated domains, proved by Cagnetti and Scardia [14] (see also the later variant in [5]).

    In the second part of this note we will specialise the general theory to some prototypical and yet relevant model cases. Namely, we will briefly review the case of perforated (or porous) brittle materials studied by Cagnetti and Scardia [14] and by Focardi, Gelli and Ponsiglione [17] (see also [5]), the case of high-contrast brittle materials with soft inclusions treated by Barchiesi, Lazzaroni and Zeppieri in [6], and eventually the case of high-contrast brittle materials with weak inclusions analysed by Pellet, Scardia and Zeppieri in [19]. In particular we will show that the choice of the integrands fε and gε in (1.1) strongly affects the form of the Γ-limit which can give rise to macroscopic models accounting for damage as well as to models accounting for cohesive fracture.

    In this part we will use the localisation method of Γ-convergence [8,15] to prove a convergence result for a general class of free-discontinuity functionals of brittle type, with degenerate coefficients.

    In the choice of the convergence to compute the Γ-limit, a crucial role will be played by an extension result for SBV-functions defined in periodically perforated domains due to Cagnetti and Scardia [14,Theorem 1.3] and by a later variant due to Barchiesi and Focardi [5,Theorem 1].

    We list below a few notation which will be used throughout the paper.

    ΩRn denotes an open and bounded set with Lipschitz boundary. The set A(Ω) denotes the collection of all open subsets of Ω;

    Q denotes the open unit cube of Rn centred at the the origin, whereas for xRn and r>0 we set Qr(x):=rQ+x;

    ● for νSn1 we denote with Qν the open unit cube of Rn centred at the the origin, with one face orthogonal to ν and for xRn and r>0 we set Qνr(x):=rQν+x;

    ● for xRn and νSn1 we denote by Πν(x) the hyperplane through x and perpendicular to ν; i.e., Πν(x):={yRn:(yx)ν=0}. If x=0 we simply write Πν;

    ● For uL1(Ω) and m>0 the function um denotes the truncated function of u at level m; i.e., um:=(um)(m);

    ● For ξRn we denote by uξ the linear function with gradient equal to ξ; i.e., uξ(x):=ξx, for every xRn;

    ● For xRn, tR, and νSn1 we denote with uν,tx the piecewise constant function taking values 0,t and jumping across the hyperplane Πν(x); i.e.,

    uν,tx(y):={tif(yx)ν0,0if(yx)ν<0.

    The functional setting we are going to consider in this note is that of SBV, the space of special functions of bounded variation. We recall here only the definition of the spaces which are relevant for our analysis and we refer the reader to [2] for a comprehensive treatment on the subject. We set

    SBV(Ω):={uBV(Ω): Du=uLn+(u+u)νudHn1Su}.

    Here Su denotes the approximate discontinuity set of u, νu is the generalised normal to Su, u+ and u are the traces of u on both sides of Su. In this paper we work with the following vector subspace of SBV(Ω)

    SBVp(Ω):={uSBV(Ω): uLp(Ω) and Hn1(Su)<+},

    where p>1. We consider also the larger space of generalised special functions of bounded variation in Ω,

    GSBV(Ω):={uL1(Ω): umSBV(Ω) for all mN},

    as well as

    GSBVp(Ω):={uGSBV(Ω): uLp(Ω) and Hn1(Su)<+}.

    We consider also

    SBVpc(Ω):={uSBV(Ω):u=0Lna.e.Hn1(Su)<+};

    it is known (see [2,Theorem 4.23]) that every u in SBVpc(Ω)L(Ω) is piecewise constant in the sense of [2,Definition 4.21], namely there exists a Caccioppoli partition (Ei) of Ω such that u is constant Ln-a.e. in each set Ei. Moreover, we set

    P(Ω):={uSBVpc(Ω):u(x){0,1} Lna.e.inΩ}.

    For u,wL1(U), in what follows, by "u=w near U" we mean that there exists a neighbourhood V of U in Rn such that u=w Ln-a.e. in VU. Let fk:Rn×Rn[0,+) be Carathéodory functions such that

    (H1) there exist p>1 and 0<c1c2<+ such that for every (x,ξ)Rn×Rn and for every kN

    c1|ξ|pfk(x,ξ)c2(1+|ξ|p); (2.1)

    (H2) fk(x,0)=0 for every xRn and for every kN.

    Let moreover gk:Rn×Sn1(0,+) be Borel functions such that

    (H3) there exist 0<c3c4<+ such that for every (x,ν)Rn×Sn1 and every kN

    c3gk(x,ν)c4; (2.2)

    (H4) gk(x,ν)=gk(x,ν), for every (x,ν)Rn×Sn1 and every kN.

    Let ΩRn be open bounded and with Lipschitz boundary and let KQ be compact and such that QK has a Lipschitz boundary. We define

    E:=RniZn(K+i);

    the set E is open, connected, Q-periodic, and has a Lipschitz boundary. Let εk be a sequence of positive numbers such that εk0 as k+ and denote by Ωk the εkQ-periodic set defined as Ωk:=ΩεkE (see Figure 1).

    Figure 1.  Schematic of a high-contrast composite material.

    Let moreover αk,βk[0,1] and consider the sequence of functionals Fk:L1(Ω)[0,+] defined as

    Fk(u):={Ωkfk(x,u)dx+αkΩΩkfk(x,u)dx+SuΩkgk(x,νu)dHn1+βkSu(ΩΩk)gk(x,νu)dHn1ifuSBVp(Ω),+otherwise in L1(Ω). (2.3)

    We observe that thanks to assumption (H2) the functionals Fk decrease by truncation, wherease they do not satisfy the standard coercivity conditions required, e.g., in [9,12,18] since the coefficients αk,βk are not bounded away from zero.

    Due to the possible degeneracy of the coefficients αk and βk, the functionals Fk are not, in general, equi-coercive with respect to the strong L1(Ω)-convergence. Similarly as in [5,11,17,19], in what follows we give a notion of convergence on L1(Ω) which is weaker that the L1(Ω)-convergence and ensures the equi-coercivity of the functionals Fk. This will be done by appealing to [5,Theorem 1]. For the readers' convenience we recall here a slightly simplified version of this result which is useful for our purposes.

    Theorem 2.1 (cf. Theorem 1 in [5]). Let (uk)SBVp(Ωk) be such that

    supk(Ωk|uk|pdx+Ωk|uk|pdx+Hn1(SukΩk))<+. (2.4)

    Then, there exist (˜uk)SBVp(Ω), with ˜uk=uk a.e. in Ωk, and a function uGSBVp(Ω)Lp(Ω) such that (up to subsequences) ˜uku in L1(Ω).

    If moreover supkukL(Ωk)<+ then uSBVp(Ω)L(Ω) and ˜uku in Lp(Ω).

    Let (uk)L1(Ω) be a sequence satisfying

    supkukLp(Ωk)<+andsupkFk(uk)<+.

    Then, clearly (uk)SBVp(Ωk); moreover in view of (H1) and (H3) the sequence (uk) satisfies the uniform bound (2.4). Therefore invoking Theorem 2.1 immediately yields the existence of a function uGSBVp(Ω)Lp(Ω) and a sequence (˜uk)SBVp(Ω) with ˜uk=uk a.e. in Ωk, such that (up to subsequences not relabelled) ˜uku in L1(Ω).

    This observation motivates the choice of the following notion of convergence on L1(Ω).

    Definition 2.2 (Convergence). Let (uk) be a sequence in L1(Ω). We say that (uk) converges to a function uL1(Ω), and we write uku, if there exists a sequence (˜uk)L1(Ω) such that ˜uk=uk a.e. in Ωk, and ˜uk converges to u in L1(Ω).

    Remark 2.3 (Uniqueness of the limit). We observe that since C(K):=Ln(QK)>0, then the limit in the sense of Definition 2.2 is well-defined. Indeed, assume that uku1 and uku2. Then by definition there exist (˜u1,k),(˜u2,k)L1(Ω) such that ˜u1,k=˜u2,k=uk in Ωk and ˜u1,ku1 and ˜u2,ku2 in L1(Ω). Therefore

    0=limk+Ωk|˜u1,k˜u2,k|dx=limk+Ω|˜u1,k˜u2,k|χΩkdx=C(K)Ω|u1u2|dx,

    where the last inequality follows by the Riemann-Lebesgue Theorem applied to the εkQ-periodic function χΩk. Then, since C(K)>0 we necessarily have u1=u2 a.e. in Ω.

    We notice moreover that the convergence uku readily implies

    limk+ukuL1(Ωk)=0.

    Remark 2.4 (Convergence of truncated functions). Let (uk)L1(Ω) be such that uku for some uL1(Ω). Let mN and denote by (umk) the sequence of truncated functions of uk at level m, then umkum where um denotes the truncated function of u at level m. Indeed, set vk:=(˜uk)m, then vk=umk a.e. in Ωk, moreover since ˜uku in L1(Ω) then (˜uk)mum in L1(Ω), and actually in any Lp(Ω).

    In what follows we study the Γ-convergence of the functionals Fk with respect to the convergence as in Definition 2.2. To this end we give the following sequential notion of Γ-convergence.

    Definition 2.5 (Sequential Γ-convergence). Let Fk,F:L1(Ω)[0,+]; we say that the functionals Fε -converge to F with respect to the convergence as in Definition 2.2 if for every uL1(Ω) the two following conditions are satisfied:

    (i) (Ansatz-free lower bound) For every (uk)L1(Ω) with uku we have

    F(u)lim infk+Fk(uk);

    (ii) (Existence of a recovery sequence) There exists (ˉuk)L1(Ω) with ˉuku such that

    F(u)lim supk+Fk(ˉuk).

    Remark 2.6. It is standard to show that F is lower semicontinuous with respect to the convergence as in Definition 2.2 and hence with respect to the strong L1(Ω)-convergence.

    For every uL1(Ω) we consider the functionals

    Γ-lim infk+Fk(u):=inf{lim infk+Fk(uk):uku} (2.5)

    and

    Γ-lim supk+Fk(u):=inf{lim supk+Fk(uk):uku}. (2.6)

    It is easy to show that the infima in (2.5) and (2.6) are actually attained.

    In what follows we also use the compact notation

    F(u):=Γ-lim infk+Fk(u) and F"(u):=Γ-lim supk+Fk(u). (2.7)

    It is immediate to see that Definition 2.5 is equivalent to F=F"=F in L1(Ω).

    Remark 2.7 (The case αk,βk=0). In the case of porous brittle materials [5,14,17], which corresponds to the parameter choice αk,βk=0, the Γ-convergence of the functionals Fk can be equivalently studied with respect to the strong L1(Ω)-convergence. Indeed, in this case a sequence (uk) with equibounded energy can be replaced by the L1(Ω)-converging sequence (˜uk) given by Theorem 2.1, without changing the energy.

    The following proposition shows that the domain of the Γ-limit of Fk (if it exists) is GSBVp(Ω).

    Proposition 2.8 (Domain of the Γ-limit). Let F and F" be as in (2.7); then

    domF=domF"=GSBVp(Ω).

    Proof. We first show that GSBVp(Ω)domF". By the growth conditions (2.1) and (2.2) we have Fk(u)G(u) where

    G(u):={c2Ω(1+|u|p)dx+c3Hn1(SuΩ)in GSBVp(Ω)+otherwise in L1(Ω) (2.8)

    The functional G is lower semicontinuous with respect to the strong L1(Ω)-convergence, hence we have

    inf{lim supk+Fk(uk):ukuinL1(Ω)}G(u).

    Since the convergence in Definition 2.2 is weaker than the L1(Ω)-convergence we then have F"≤G, and thus the desired inclusion.

    We now prove that domFGSBVp(Ω). To this end, let udomF then there exists (uk)L1(Ω) with uku such that lim infkFk(uk)=F(u)<+. Then, up to subsequences (not relabelled) we have supkF(uk)<+, thus in particular (uk)SBVp(Ω).

    Let mN and let umk be the truncated function of uk at level m; then (umk)SBVp(Ω)L(Ω). Since the functionals Fk decrease by truncation, for every fixed mN it also holds supkFk(umk)<+. Therefore, for mN fixed we can appeal to Theorem 2.1 to deduce the existence of a sequence (vk)L1(Ω) such that vk=umk a.e. in Ωk and of a function vSBVp(Ω) such that up to subsequence (not relabelled) vkv in L1(Ω). Since uku we have

    0=limk+Ωk|vkumk|dx=limk+Ω|vk(˜uk)m|χΩkdx=C(K)Ω|vum|dx,

    therefore v=um a.e. in Ω. Eventually, the arbitrariness of mN yields uGSBVp(Ω).

    In this section we show that, up to subsequences, the functionals Fk Γ-converge to a free-discontinuity functional of the form

    F(u)=Ωf(x,u)dx+Sug(x,u+u,νu)dHn1.

    for some f and g. Moreover, we show that, despite the degeneracy of the coefficients αk,βk, the limit integrands f and g satisfy standard coercivity conditions similar to (2.1) and (2.2), respectively.

    If not otherwise specified, in what follows the Γ-convergence of the functionals Fk is always understood in the sense of Definition 2.5.

    To prove the existence of a Γ-convergent subsequence of Fk we make use of the so-called localisation method [8,15] which we adapt to the sequential notion of Γ-convergence as in Definition 2.5.

    We start by localising the functionals Fk; that is we consider Fk:L1(Ω)×A(Ω)[0,+] defined as

    Fk(u,U):={Ukfk(x,u)dx+αkUUkfk(x,u)dx+SuUkgk(x,νu)dHn1+βkSu(UUk)gk(x,νu)dHn1ifuSBVp(U),+otherwise in L1(Ω), (2.9)

    where Uk:=UεkE.

    We also define the localised versions of (2.5) and (2.6); i.e., for every UA(Ω) we consider the functionals defined as

    F(,U):=Γ-lim infk+Fk(,U),F"(,U):=Γ-lim supk+Fk(,U). (2.10)

    Remark 2.9 (Properties of F,F"). It is easy to show that F and F" are lower semicontinuous with respect to the convergence in Definition 2.2, local, and that they decrease by truncation. Moreover, as set functions they are both increasing, whereas F is also superadditive.

    Remark 2.10 (On assumption (H2)). If we drop assumption (H2) the functionals Fk will not decrease by truncation, but rather satisfy

    Fk(um,U)Fk(u,U)+c2Ln(Uk{|u|m})+αkc2Ln(UUk{|u|m}). (2.11)

    If αk is infinitesimal, the inequality in (2.11) implies

    F(um,U)F(u,U)+c2muL1(Ω), (2.12)

    (and analogously for F"). In fact, by definition of Γ-liminf there exists a sequence (uk)L1(Ω) such that uku and F(u,U)=lim infk+Fk(uk,U). Then if umk is the truncated function of uk at level m, by (2.11) we get

    Fk(umk,U)Fk(uk,U)+c2Ln(Uk{|uk|m})+αkc2Ln(UUk{|uk|m})Fk(uk,U)+c2Ln(Uk{|˜uk|m})+αkc2Ln(Ω)Fk(uk,U)+c2Ln(U{|˜uk|m})+αkc2Ln(Ω)

    where ˜uk is as in Definition 2.2 and thus ˜uku in L1(Ω). Therefore, taking the liminf as k+ gives

    lim infk+Fk(umk,U)F(u,U)+c2Ln(U{|u|m}),

    hence (2.12) follows by the definition of Γ-liminf, taking into account that umkum, and by the Chebyshev inequality. Therefore, (2.12) ensures that F "almost" decreases by truncation up to an error which becomes small for m large. Inequality (2.12) is then enough to to carry out the Γ-convergence analysis below (cf. [12]). Hence, if αk is infinitesimal assumption (H2) can be dropped.

    However, if the sequence αk is uniformly bounded from below, we have no control on the term αkc2Ln(UUk{|uk|m}), therefore from (2.11) we cannot infer (2.12). Since with we want to study the Γ convergence of Fk for any choices of αk[0,1], assumption (H2) is actually necessary.

    In general the set functions F(u,) and F"(u,) are not inner regular. Then we consider their inner regular envelopes defined as:

    F(u,U):=sup{F(u,V):V⊂⊂U,VA(Ω)}.

    and

    F"(u,U):=sup{F"(u,V):V⊂⊂U,VA(Ω)}.

    Remark 2.11 (Properties of F,F"). The functionals F and F" are lower semicontinuous with respect to the convergence in Definition 2.2 [15,Remark 15.10], local [15,Remark 15.25], and it is immediate to check that they decrease by truncation. Furthermore, as set functions, they are both increasing and F is superadditive [15,Remark 15.10].

    The following compactness result is the analogue of [15,Theorem 16.9], when the sequential notion of Γ-convergence in Definition 2.5 is considered. We omit its proof since it is standard.

    Proposition 2.12 (Compactness by Γ-convergence). Let Fk be the localised functionals as in (2.9). Then there exists a subsequence (Fkj)(Fk) such that the corresponding functionals F and F" defined in (2.10) satisfy F=F".

    We now set

    F:=F=F". (2.13)

    In what follows we show that actually F coincides with the Γ-limit of the subsequence (Fkj). To this end we start noticing that by monotonicity we always have F"=FFF". Therefore, if we show that F"=F"; i.e., that F" is inner regular, we immediately get F=F"=F and therefore that Fkj(,U) Γ convergesto {\mathcal F}(\cdot, U) forevery U\in A(\Omega) $, as desired.

    A crucial preliminary result needed to prove the inner-regularity of F" is the so-called fundamental estimate, which has to hold uniformly in k. Since the Γ-limit is computed with respect to the convergence in Definition 2.2, the fundamental estimate we need is non-standard. Namely, we have to prove that the error in the fundamental estimate tends to zero when uku. This is achieved by first showing that the error goes like ukuLp(Ωk) and then by resorting to a truncation argument.

    We notice that an analogous estimate for degenerate functionals defined in Sobolev spaces can be found in [10,Proposition 3.3]. Whereas in the SBV-setting, for functionals of Mumford-Shah type with degenerate surface energy it can be found in the recent [19,Lemma 4.4].

    Following [10] we start showing how to construct suitable cut-off functions which are constant in iZnεk(K+i). To this end let δ>0 be small enough so that the set Kδ:={xRn:dist(x,K)<δ} satisfies Kδ⊂⊂Q. Let ψC0(Q) be a cut-off function between K and Kδ (that is 0ψ1, ψ1 on K, and sptψKδ) such that |ψ|2η.

    For kN and iZn, we define the operator Rki:W1,loc(Rn)W1,loc(Rn) as

    Rki(ϕ)(x):=(1ψ(xεki))ϕ(x)+ψ(xεki)εkKδ+εkiϕ(y)dy.

    By definition we have that

    Rki(ϕ)(x)=ϕ(x)if xεkKδ+εki,

    while Rki is constant in εkK+εki, namely we have

    Rki(ϕ)(x)=εkKδ+εkiϕ(y)dyif xεkK+εki.

    Finally, we consider the operator Rk:W1,loc(Rn)W1,loc(Rn) defined as

    Rk(ϕ)(x):={Rki(ϕ)(x)ifxεkKδ+εki,iZn,ϕ(x)otherwise.

    Let URn be open and bounded and let ϕW1,(U) then Rk(ϕ) is uniformly bounded in k. More precisely, we have

    Rk(ϕ)L(U;Rn)(2δd+1)ϕL(U;Rn), (2.14)

    where d denotes the diameter of Kδ. In fact,

    Rk(ϕ)L(U;Rn)2εkδsupiϕεKδ+εiϕ(y)dyL(εkKδ+εki;Rn)+ϕL(U;Rn)

    and

    ϕεkKδ+εkiϕ(y)dyL(εkKδ+εki;Rn)εkdϕL(U;Rn).

    In the next proposition we make use of the operator Rk to construct cut-off functions whose gradient vanishes in RnεkE; these cut-off functions are then used to prove the desired fundamental estimate.

    Proposition 2.13 (Fundamental estimate). For every η>0, and for every U, U", VA(Ω), with U⊂⊂U", there exist two constants M(η)>0 and kηN satisfying the following property: For every k>kη, for every uL1(Ω) with uSBVp(U"), and for every vL1(Ω) with vSBVp(V), there exists a function φC0(Ω) with φ=1 in a neighbourhood of U, sptφU" and 0φ1 such that

    Fk(φu+(1φ)v,UV)(1+η)(Fk(u,U")+Fk(v,V))+M(η)||uv||Lp(SεkE) (2.15)

    where S:=(U"U)V.

    Proof. Let U,U",VA(Ω) be as in the statement. Let η>0 be fixed and choose NN in a way such that

    1Nmax{c2Ln((U"U)V),3p1max{c2c1,c4c3}}<η. (2.16)

    Let moreover UA(Ω) be such that U⊂⊂U⊂⊂U" and consider the open sets

    U⊂⊂U1⊂⊂⊂⊂U3N⊂⊂U"

    where

    Ul:={x:dist(x,U)<dist(U,U)3Nl},for everyl=1,,3N.

    We notice that by definition of Ul we have that

    dist(Ul,Ul+1)=13N,for everyl=1,,3N1. (2.17)

    For every j=0,,N1 let ϕj be a cut-off function between U3j+1 and U3j+2 with |ϕj|<4N.

    Let kηN be such that

    εkd<13Nfor everyk>kη, (2.18)

    where d:=diam(Kδ)<2.

    If iZn is such that (εkK+εki)Ul for every k>kη, then thanks to (2.17)–(2.18) we can deduce that (εkK+εki)(RnUl+1)=. Therefore the functions φj:=Rk(ϕj) are cut-off functions between the sets U3j and U3(j+1), for every j=0,,N1 (where we have set U0:=U).

    Now let uSBVp(U") and vSBVp(V); for every j=0,,N1 fixed we have

    Fk(φju+(1φj)v,UV)=Fk(u,(UV)¯U3j)+Fk(v,VU3(j+1))+Fk(φju+(1φj)v,V(U3(j+1)¯U3j))Fk(u,U")+Fk(v,V)+Fk(φju+(1φj)v,V(U3(j+1)¯U3j)).

    We set

    wj:=φju+(1φj)v,Sj:=V(U3(j+1)¯U3j)

    and estimate the term Fk(wj,Sj). We clearly have

    Fk(wj,Sj)=Fk(wj,SjεkE)+Fk(wj,SjεkE). (2.19)

    By construction φj=0 in RnεkE, therefore appealing to (2.1) and (2.2) we deduce

    Fk(wj,SjεkE)=αkSjεkEfk(x,φju+(1φj)v)dx+βk(SjεkE)Swjgk(x,νwj)dHn1c2αk(Ln(SjεkE)+SjεkE|u|pdx+SjεkE|v|pdx)+c4βk(Hn1((SjεkE)Su)+Hn1((SjεkE)Sv))c2αkLn(SjεkE)+c2c1(αkSjεkEfk(x,u)dx+αkSjεkEfk(x,v)dx)+c4c3(βk(SjεkE)Sugk(x,νu)dHn1+βk(SjεkE)Svgk(x,νv)dHn1)c2αkLn(SjεkE)+max{c2c1,c4c3}(Fk(u,SjεkE)+Fk(v,SjεkE)). (2.20)

    Moreover, again invoking (2.1) and (2.2), in εkE we have

    Fk(wj,SjεkE)=SjεkEfk(x,wj)dx+SjεkESwjgk(x,νwj)dHn1c2(Ln(SjεkE)+SjεkE|φj(uv)+φju+(1φj)v|pdx)+c4(Hn1((SjεkE)Su)+Hn1((SjεkE)Sv))c2Ln(SjεkE)+c23p1φjpL(U;Rn)SjεkE|uv|pdx+3p1max{c2c1,c4c3}(Fk(u,SjεkE)+Fk(v,SjεkE)). (2.21)

    Since |ϕj|4N, combining the definition of φj with (2.14) gives

    φjL(U;Rn)(2δd+1)4N. (2.22)

    In view of (2.19), by gathering (2.20)–(2.22) we then obtain for every j=0,,N1

    Fk(wj,Sj)c2Ln(Sj)+3p1max{c2c1,c4c3}(Fk(u,Sj)+Fk(v,Sj))+M(η)SjεkE|uv|pdx,

    where

    M(η):=c23p1(2δd+1)p(4N)p.

    Therefore there exists j{0,,N1} such that

    Fk(wj,Sj)1NN1i=0Fk(wj,Sj)c2NLn((U"U)V)+3p1Nmax{c2c1,c4c3}(Fk(u,U")+Fk(v,V))+M(η)(U"U)VεkE|uv|pdx.

    Finally the thesis follows from (2.16) by choosing φj as a cut-off function and setting S:=(U"U)V.

    Thanks to the fundamental estimate Proposition 2.13 we are now able to prove the following abstract Γ-convergence result for the sequence of localised functionals Fkj(,U).

    Theorem 2.14 (Abstract Γ-convergence and properties of the Γ-limit). Let F be as in (2.13), then:

    a. (locality and lower semicontinuity) for every UA(Ω), the functional F(,U) is local and lower semicontinuous with respect to the L1(Ω)-convergence;

    b. (measure property) for every uGSBVp(Ω), the set function F(u,) is the restriction to A(Ω) of a Radon measure on Ω;

    c. (Γ-convergence) for every UA(Ω) it holds F(,U)=F(,U)=F"(,U)on GSBVp(Ω)

    d. (translational invariance in u) for every uL1(Ω) and UA(Ω) there holds F(u+s,U)=F(u,U) for every sR.

    Proof. Since the L1(Ω)-convergence implies the convergence in the sense of Definition 2.2, property 1 immediately follows from Remark 2.11. In view of Remark 2.11, property 2 follows by the De Giorgi and Letta criterion (see, e.g., [15,Theorem 14.23]) once we show that for every uGSBVp(Ω) the set function F(u,) is subadditive. In its turn, the subadditivity of F(u,) follows from Proposition 2.13. Since in our setting this proof is not entirely standard, we discuss it in detail for the readers' convenience.

    We start observing that on GSBVp(Ω) the functional F satisfies the following limsup-type inequality: For every uGSBVp(Ω) and for every U,UA(Ω) with U⊂⊂U, there exists a sequence (uj)GSBVp(U)L1(Ω) with uju such that

    lim supj+Fkj(uj,U)F(u,U)

    (see, e.g., [15,Proposition 16.4 and Remark 16.5] also recalling that the infimum in the definition of F" is actually attained).

    Now let U,VA(Ω) and let uGSBVp(Ω)L(Ω). Fix any U⊂⊂U, V⊂⊂V, U,VA(Ω). Choose an open set U" such that U⊂⊂U"⊂⊂U and two sequences (uj)GSBVp(U")L1(Ω) and (vj)GSBVp(V)L1(Ω), with uju and vju such that

    lim supj+Fkj(uj,U")F(u,U),lim supj+Fkj(vj,V)F(u,V). (2.23)

    Since the functionals Fk decrease by truncation, we can additionally assume that ujL(Ω),vjL(Ω)uL(Ω); clearly,

    limj+ujuLp(Ωkj)=limj+vjuLp(Ωkj)=0. (2.24)

    Let η>0 be fixed and arbitrary. The fundamental estimate Proposition 2.13 provides us with constants M(η)>0 and jηN and with a sequence (φj) of cut-off functions between U and U" such that

    Fkj(φjuj+(1φj)vj,UV)(1+η)(Fkj(uj,U")+Fkj(vj,V))+M(η)ujvjpLp(Ωkj)

    for every jjη. Hence appealing to (2.23), to the convergence φjuj+(1φj)vju, and to the obvious inequality FF, by taking the limit as j+, we get

    F(u,UV)(1+η)(F(u,U)+F(u,V)).

    Now letting η0, and then UU, VV in view of the inner-regularity of F we get

    F(u,UV)F(u,U)+F(u,V), (2.25)

    hence the subadditivity of F(u,) for uGSBVp(Ω)L(Ω).

    Now let uGSBVp(Ω) and, for every mN, set um:=(um)(m). Then, since F decreases by truncation (2.25) immediately gives

    F(um,UV)F(u,U)+F(u,V).

    Then, taking the limit as m+, in view of the convergence umu in L1(Ω) and the lower semicontinuity of F we obtain

    F(u,UV)lim infm+F(um,UV)F(u,U)+F(u,V),

    and thus the subadditivity of F(u,) for every uGSBVp(Ω).

    The proof of property 3 is achieved by showing that F" is inner-regular. Indeed, this is equivalent to F"=F", which by definition of F implies F"≤FF. Since clearly FF", we actually deduce that F is the Γ-limit of Fkj.

    The inner regularity of F" follows from the fundamental estimate Proposition 2.13. To see this, for every UA(Ω) let G(,U) be the localised version of the functional G defined in (2.8); i.e.,

    G(u,U):={c2U(1+|u|p)dx+c4Hn1(SuU)ifuGSBVp(U),+otherwise in L1(Ω). (2.26)

    Now fix WA(Ω) and uGSBVp(Ω); since G(u,) is the restriction to A(Ω) of a Radon measure, for every η>0 there exists a compact set ˜WW such that and MS(u,W˜W)<η.

    Now choose U,UA(Ω) satisfying ˜WU⊂⊂U⊂⊂W and set V:=W˜W. Recalling that F"(u,) is increasing, appealing to Proposition 2.13 easily gives

    F"(u,W)F"(u,UV)F"(u,U)+F"(u,V)=F"(u,U)+F"(u,W˜W).

    Recalling that F"≤G, by taking the sup on U⊂⊂W we get

    F"(u,W)F"(u,W)+G(u,W˜W)F"(u,W)+η.

    Hence, by the arbitrariness of η>0 we get F"(u,W)F"(u,W) for every WA(Ω) and every uGSBVp(Ω). Since the opposite inequality is always satisfied, we readily deduce the inner regularity of F"(u,), as desired.

    Eventually, the proof of property 4 is standard and follows as in, e.g., [9,Lemma 3.7].

    In the following theorem we show that the Γ-limit F can be represented in an integral form as a free-discontinuity functional. Moreover, thanks to [5,Theorem 4] the functional F turns out to be non-degenerate, unlike the functionals Fk.

    Theorem 2.15 (Integral representation of the Γ-limit). Let F be as in Theorem 2.14. Then, there exist a Carathéodory function f:Rn×Rn[0,+) and a Borel function g:Rn×Sn1(0,+) such that

    F(u,U)=Uf(x,u)dx+SuUg(x,[u],νu)dHn1 (2.27)

    for every uGSBVp(Ω) and every UA(Ω).

    Furthermore, the function f:Rn×Rn[0,+) satisfies the following properties:

    i) (convexity in ξ) for a.e. xRn, f(x,) is convex;

    ii) (p growth and coercivity) there exists ˜c1>0 such that for a.e. xRn and for every ξRn it holds

    ˜c1|ξ|pf(x,ξ)c2(1+|ξ|p), (2.28)

    where c2 is as in (2.1).

    The function g:Rn×R×Sn1[0,+) satisfies the following properties:

    iii) (monotonicity in t and symmetry) for a.e. xRn and for every νSn1, g(x,,ν) is nondecreasing on (0,+) and satisfies the symmetry condition g(x,t,ν)=g(x,t,ν) for every tR;

    iv) (subadditivity in t) for a.e. xRn and for every νSn1

    g(x,t1+t2,ν)g(x,t1,ν)+g(x,t2,ν),

    for every t1,t2R;

    v) (convexity in ν) for a.e. xRn and for every tR, the 1-homogeneous extension of g(x,t,) to Rn is convex. Equivalently, for a.e. xRn and for every tR the function g satisfies

    g(x,t,ν)λ1g(x,t,ν1)+λ2g(x,t,ν2),

    for every ν,ν1,ν2Sn1, λ1,λ20 such that λ1ν1+λ2ν2=ν;

    vi) (bounds) there exists ˜c3>0 such that for a.e. xRn, for every tR, and every νSn1 it holds

    ˜c3g(x,t,ν)c4, (2.29)

    where c4 is as in (2.2).

    Proof. Let Ek:L1(Ω)×A(Ω)[0,+] be the functionals defined as

    Ek(u,U):={c1Uk|u|pdx+c3Hn1(SuUk)ifuSBVp(U)+otherwise in L1(Ω), (2.30)

    with c_1 and c_3 as in (2.1) and (2.2), respectively. Appealing to [5,Theorem 4] and also noticing that the L^p -convergence in the statement can be equivalently replaced by the convergence in Definition 2.2, we deduce that \mathcal E_{k}(\cdot, U) \Gamma -converges to \mathcal E(\cdot, U) for every U\in \mathcal{A}(\Omega) , where

    \begin{equation*} \mathcal E(u,U) = \int_{U}\hat f(\nabla u) \,dx +\int_{S_u\cap U} \hat g(\nu)\,d\mathcal H^{n-1} \end{equation*}

    with \hat f and \hat g as in [5,Theorem 4] formulas (40) and (41), respectively. Moreover \hat f and \hat g satisfy

    \tilde c_1|\xi|^p \leq \hat f(\xi)\; \text{for every $\xi\in {\mathbb R}^n$} \quad \text{and}\quad \tilde c_3 \leq \hat g(\nu)\; \text{for every $\nu\in {\mathbb S}^{n-1}$},

    for some \tilde c_1, \tilde c_3 > 0 . Then, since \mathcal E_k \leq {\mathcal F}_k , we may deduce that for every u\in SBV^p(\Omega) and every U\in \mathcal{A}(\Omega) we have

    \begin{equation} \mathcal E(u,U) \leq {\mathcal F}(u,U). \end{equation} (2.31)

    We recall that for every u\in SBV^p(\Omega) and every U\in \mathcal{A}(\Omega) we also have

    \begin{equation} {\mathcal F}(u,U) \leq \mathcal G(u,U), \end{equation} (2.32)

    where \mathcal G is as in (2.26).

    Now let \sigma > 0 and for every u\in SBV^p(\Omega) and U\in \mathcal{A}(\Omega) set

    {\mathcal F}^\sigma(u,U): = {\mathcal F}(u,U)+\sigma \int_{S_u \cap U}|[u]|\,d\mathcal H^{n-1}.

    For every fixed \sigma > 0 the functional {\mathcal F}^\sigma satisfies properties 1, 2, and 4 in Theorem 2.14. Moreover, in view of (2.31)–(2.32) it holds

    \begin{multline*} \tilde c_1\int_{U} |\nabla u|^p \,dx +\int_{S_u \cap U}(\tilde c_3 +\sigma |[u]|)\,d\mathcal H^{n-1} \leq {\mathcal F}^\sigma(u,U) \\ \leq c_2\int_{U} (1+|\nabla u|^p) \,dx + \int_{S_u \cap U} (c_4 +\sigma |[u]|)\,d\mathcal H^{n-1} \end{multline*}

    Therefore, we can invoke the integral representation result [7,Theorem 1] to deduce that for every u\in SBV^p(\Omega) and every U\in \mathcal{A}(\Omega) we have

    {\mathcal F}^\sigma(u,U) = \int_U f^\sigma_{\infty}(x,\nabla u) \,dx +\int_{S_u\cap U}g^\sigma_{\infty}(x,[u],\nu_u)\, d\mathcal H^{n-1},

    where f^\sigma_\infty and g^\sigma_\infty are given by the following derivation formulas

    \begin{equation} f^\sigma_\infty (x,\xi): = \limsup\limits_{\rho \to 0^+}\frac{1}{\rho^n}\inf\big\{ {\mathcal F}^\sigma(u,Q_\rho(x)) \colon u\in SBV^p(Q_\rho(x)), u = u_\xi \; \text{near}\; \partial Q_\rho(x)\big\} \end{equation} (2.33)

    and

    \begin{equation} g^\sigma_\infty (x,t,\nu): = \limsup\limits_{\rho \to 0^+}\frac{1}{\rho^{n-1}}\inf\big\{ {\mathcal F}^\sigma(u,Q^\nu_\rho(x)) \colon u\in SBV^p(Q^\nu_\rho(x)), u = u^{t,\nu}_x \; \text{near}\; \partial Q^\nu_\rho(x)\big\}. \end{equation} (2.34)

    By (2.33) and (2.34) the sequences (f^\sigma_\infty)_{\sigma > 0} and (g^\sigma_\infty)_{\sigma > 0} are decreasing as \sigma decreases, therefore by setting f_\infty: = \lim_{\sigma\to0^+}f_\infty^\sigma and g_\infty: = \lim_{\sigma\to0^+}g_\infty^\sigma , by the pointwise convergence of ({\mathcal F}^\sigma)_{\sigma > 0} to {\mathcal F} and the Monotone Convergence Theorem, we get

    {\mathcal F}(u,U) = \int_U f_{\infty}(x,\nabla u) \,dx +\int_{S_u\cap U}g_{\infty}(x,[u],\nu_u)\, d\mathcal H^{n-1},

    for every u\in SBV^p(\Omega) and U\in \mathcal{A}(\Omega) . Eventually, a standard truncation and continuity argument allows to extend this integral representation to the whole space GSBV^p(\Omega) and thus to get exactly (2.27).

    The measurability properties of f_\infty and g_\infty follow from the derivation formulas (2.33) and (2.34), arguing as in the appendix of [12]. The convexity of f_\infty in \xi , the subadditivity of g_\infty in t , and the convexity in \nu of its 1 -homogeneous extension are immediate consequences of the L^1(\Omega) -lower semicontinuity of {\mathcal F} .

    To show that f_\infty and g_\infty satisfy, respectively, the lower bounds as in ii) and vi) we argue as follows. Set

    \Phi^\sigma(u,U): = \begin{cases} \tilde c_1\int_{U} |\nabla u|^p \,dx +\int_{S_u \cap U}(\tilde c_3 +\sigma |[u]|)\,d\mathcal H^{n-1} & \text{if}\; u\in SBV^p(U) \cr +\infty & \text{otherwise in }\, L^1( \Omega), \end{cases}

    and for every x\in {\mathbb R}^n and \xi\in {\mathbb R}^n define

    \phi^\sigma(x,\xi): = \limsup\limits_{\rho \to 0^+}\frac{1}{\rho^n}\inf\big\{\Phi^\sigma(u,Q_\rho(x)) \colon u\in SBV^p(Q_\rho(x)), u = u_\xi \; \text{near}\; \partial Q_\rho(x)\big\},

    while for every x\in {\mathbb R}^n , t\in {\mathbb R} , and \nu\in {\mathbb S}^{n-1} set

    \psi^\sigma(x,t,\nu): = \limsup\limits_{\rho \to 0^+}\frac{1}{\rho^{n-1}}\inf\big\{\Phi^\sigma(u,Q^\nu_\rho(x)) \colon u\in SBV^p(Q^\nu_\rho(x)), u = u^{\nu,t}_x \; \text{near}\; \partial Q^\nu_\rho(x)\big\}.

    Since \Phi^\sigma \leq {\mathcal F}^\sigma on SBV^p(\Omega) we clearly have both \phi^\sigma \leq f_\infty^\sigma and \psi^\sigma \leq g_\infty^\sigma . We now show that \phi^\sigma(x, \xi) = \tilde c_1 |\xi|^p for every x\in {\mathbb R}^n and every \xi\in {\mathbb R}^n and \psi^\sigma(x, t, \nu) = \tilde c_3 +\sigma t . To do so we notice that by the homogeneity in x of \Phi^\sigma , we have both \phi^\sigma(x, \xi) = \phi^\sigma(0, \xi) for every x\in {\mathbb R}^n and every \xi\in {\mathbb R}^n and \psi^\sigma(x, t, \nu) = \psi^\sigma(0, t, \nu) . We can now apply the integral representation result [7,Theorem 1] to \Phi^\sigma so that choosing u = u_\xi and U = Q we obtain

    \tilde c_1 |\xi|^p = \Phi^\sigma (u_\xi, Q) = \int_Q \phi^\sigma (y,\xi) \,dy = \phi^\sigma(0,\xi) = \phi^\sigma(x,\xi),

    while choosing u = u^{t, \nu}_0 and U = Q^\nu we obtain

    \tilde c_3 +\sigma t = \Phi^\sigma (u^{t,\nu}_0, Q^\nu) = \int_{\Pi^\nu \cap Q^\nu} \psi^\sigma (y,t,\nu)\,d\mathcal H^{n-1} = \psi^\sigma(0,t,\nu) = \psi^\sigma(x,t,\nu),

    and hence the desired equalities. Therefore we deduce

    \tilde c_1 |\xi|^p = \phi^\sigma(x,\xi) \leq f_\infty^\sigma(x,\xi) \quad \text{for every }\, x,\xi\in {\mathbb R}^n

    which immediately gives the lower bound of f_\infty ; moreover there holds

    \tilde c_3 \leq \tilde c_3 +\sigma t = \psi^\sigma(x,t,\nu) \leq g_\infty^\sigma(x,t,\nu) \quad \text{for every }\, x\in {\mathbb R}^n, t\in {\mathbb R}, \nu\in {\mathbb S}^{n-1}

    hence, taking the inf on \sigma > 0 yields the the lower bound on g_\infty .

    The upper bound in ii) immediately follows from (2.33) and the obvious inquality {\mathcal F}^\sigma(u_\xi, Q_\rho(x)) \leq \rho^n c_2(1+|\xi|^p) , while the upper bound in vi) follows from (2.34) and

    {\mathcal F}(u_x^{t,\nu}, Q^\nu_\rho(x))\leq {\mathcal F}^\sigma(u_x^{t,\nu}, Q^\nu_\rho(x)) \leq \rho^{n-1} (c_4+\sigma t),

    which holds true for every \sigma > 0 and hence also in the limit as \sigma \to 0^+ .

    Finally, the monotonicity in t and the symmetry of g_\infty easily follow from (2.34).

    Theorem 2.16 ( \Gamma -convergence). Let {\mathcal F}_k be the functionals defined in (2.3). Then, there exists a subsequence k_j \to +\infty such that ({\mathcal F}_{k_j}) \Gamma -converges to the functional {\mathcal F} given by (2.27), for some Carathéodory function f_\infty \colon {\mathbb R}^n \times {\mathbb R}^n\to [0, +\infty) and some Borel function g_\infty \colon {\mathbb R}^n \times {\mathbb S}^{n-1}\to (0, +\infty) satisfying properties (i)-(vi) as in Theorem 2.15.

    Proof. The proof is an immediate consequence of Theorem 2.14 and Theorem 2.15.

    Corollary 2.17 ( \Gamma -convergence of porous brittle materials). Let \alpha_ k = \beta _ k = 0 and let {\mathcal F}_k be the corresponding functionals given by (2.3). Then, there exists a subsequence k_j \to +\infty such that ({\mathcal F}_{k_j}) \Gamma -converges with respect to the L^1(\Omega) -convergence to the functional {\mathcal F} given by (2.27), for some Carathéodory function f_\infty \colon {\mathbb R}^n \times {\mathbb R}^n\to [0, +\infty) and some Borel function g_\infty \colon {\mathbb R}^n \times {\mathbb S}^{n-1}\to (0, +\infty) satisfying properties (i)-(vi) as in Theorem 2.15.

    Proof. Since the L^1(\Omega) -convergence implies the convergence in Definition 2.2, the proof of the liminf inequality is immediate from Theorem 2.16. Now let u\in GSBV^p(\Omega) , then by Theorem 2.16 there exists (u_j)\subset L^1(\Omega) such that u_j \leadsto u and \lim_j {\mathcal F}_{k_j}(u_j) = {\mathcal F}(u) . In view of Definition 2.2 this means that there exists a sequence (\tilde u_j) \subset L^1(\Omega) such that \tilde u_j = u_j a.e. in \Omega_{k_j} and \tilde u_j \to u in L^1(\Omega) . Then, since the choice \alpha _{k_j} = \beta _{k_j} = 0 implies the equality {\mathcal F}_{k_j}(\tilde u_j) = {\mathcal F}_{k_j}(u_j) , the sequence (\tilde u_j) is the desired recovery sequence.

    On account of the \Gamma -convergence result Theorem 2.14 in this section we establish a convergence result for minimisation problems associated to a suitable perturbation of the functionals {\mathcal F}_k . To this end, let h\in L^\infty(\Omega) and for every k set

    M_{k}: = \inf\big\{ {\mathcal F}_{k}(u)+ \|u-h\|^p_{L^p( \Omega_k)}: u\in L^1( \Omega) \big\}.

    By a standard truncation argument it is immediate to show that

    \begin{equation} M_{k} = \inf\big\{ {\mathcal F}_{k}(u)+ \|u-h\|^p_{L^p( \Omega_k)}: u\in SBV^p( \Omega), \; \|u\|_{L^\infty( \Omega)}\leq \|h\|_{L^\infty( \Omega)} \big\}. \end{equation} (2.35)

    Proposition 2.18. Let {\mathcal F} = \Gamma\hbox{-}\lim_{j} {\mathcal F}_{k_j} and let (u_j) \subset SBV^p(\Omega) be such that

    \begin{equation} \lim\limits_{j \to +\infty}\big( {\mathcal F}_{k_j}(u_j)+\|u_j-h\|^p_{L^p( \Omega_{k_j})}-M_j\big) = 0. \end{equation} (2.36)

    Then, up to subsequences (not relabelled), u_j converges in the sense of Definition 2.2 to a function \bar u \in SBV^p(\Omega)\cap L^\infty(\Omega) which solves

    M: = \min\big\{ {\mathcal F}(u)+ C(K)\|u-h\|^p_{L^p( \Omega)}: u\in SBV^p( \Omega), \; \|u\|_{L^\infty( \Omega)}\leq \|h\|_{L^\infty( \Omega)} \big\},

    where C(K): = \mathcal L^n(Q\setminus K) . Moreover it holds M_{j} \to M , as j\to +\infty .

    Proof. Let (u_j)\subset SBV^p(\Omega) be as in (2.36). Then, in view of (2.35), (H1), and (H3) we have

    \begin{equation*} \sup\limits_j \Big(\|u_j\|_{L^\infty( \Omega)}+\int_{ \Omega_{k_j}} |\nabla u_j|^p \,dx+\mathcal H^{n-1}(S_{u_j}\cap \Omega_{k_j})\Big) \lt +\infty. \end{equation*}

    Therefore Theorem 2.1 yields the existence of a function \bar u\in SBV^p(\Omega)\cap L^\infty(\Omega) and of a sequence (\tilde u_j) \subset SBV^p(\Omega) with \tilde u_j = u_j a.e. in \Omega_{k_j} such that (up to subsequences) \tilde u_j \to \bar u in L^p(\Omega) , moreover \|\bar u\|_{L^\infty(\Omega)}\leq \|h\|_{L^\infty(\Omega)} . We have

    C(K)\|\bar u -h\|^p_{L^p( \Omega)} = \lim\limits_{j \to +\infty}\|(\tilde u_j -h)\chi_{ \Omega_{k_j}}\|^p_{L^p( \Omega)} = \lim\limits_{j\to +\infty}\|u_j -h\|^p_{L^p( \Omega_{k_j})},

    thus by Theorem 2.14 we get

    \begin{eqnarray*} {\mathcal F}(\bar u)+ C(K)\|\bar u -h\|^p_{L^p( \Omega)}\leq \liminf\limits_{j \to +\infty} \Big( {\mathcal F}_{k_j}(u_j)+\|u_j -h\|^p_{L^p( \Omega_{k_j})}\Big). \end{eqnarray*}

    Therefore, by definition of u_j we obtain

    \begin{equation} {\mathcal F}(\bar u)+ C(K)\|\bar u -h\|^p_{L^p( \Omega)}\leq \liminf\limits_{j \to +\infty}M_j. \end{equation} (2.37)

    Now let w\in SBV^p(\Omega)\cap L^\infty(\Omega) be an arbitrary function such that \|w\|_{L^\infty(\Omega)}\leq \|h\|_{L^\infty(\Omega)} . Again appealing to Theorem 2.14 we can find (w_j) \subset L^1(\Omega) such that w_j \leadsto w and \lim_{j} {\mathcal F}_j(w_j) = {\mathcal F}(w) . Now let \tilde w_j be as in Definition 2.2, let m: = \|h\|_{L^\infty(\Omega)} and denote with (\tilde w_j^m) the sequence of truncated functions of (\tilde w_j) at level m . We clearly have \tilde w^m_j = w^m_j a.e. in \Omega_{k_j} and \tilde w^m_j \to w in L^p(\Omega) . Hence

    \lim\limits_{j\to +\infty}\|w_j^m -h\|^p_{L^p( \Omega_{k_j})} = \lim\limits_{j \to +\infty}\|(\tilde w_j^m -h)\chi_{ \Omega_{k_j}}\|^p_{L^p( \Omega)} = C(K)\|w -h\|^p_{L^p( \Omega)}.

    Moreover, since \limsup_{j} {\mathcal F}_j(w^m_j) \leq {\mathcal F}(w) , we immediately deduce

    \begin{eqnarray} \limsup\limits_{j \to +\infty} M_j \leq {\mathcal F}(w)+C(K)\|w -h\|^p_{L^p( \Omega)}. \end{eqnarray} (2.38)

    Finally, by gathering (2.37) and (2.38) we obtain

    \begin{eqnarray*} & {\mathcal F}(\bar u)+ C(K)\|\bar u -h\|^p_{L^p( \Omega)} \leq \liminf\limits_{j \to +\infty} M_j \leq \limsup\limits_{j \to +\infty} M_j & \\ & \leq {\mathcal F}(w)+ C(K)\|w -h\|^p_{L^p( \Omega)},& \end{eqnarray*}

    hence by the arbitrariness of w we deduce that \bar u is a minimiser for {\mathcal F}+C(K)\|\cdot -h\|^p_{L^p(\Omega)} . Finally, taking w = \bar u also implies M_j \to M . Since moreover this limit does not depend on the subsequence, the convergence holds true for the whole (M_j) .

    In this section we restrict the analysis to the case of \varepsilon_k -periodic integrands f_k and g_k . That is, we consider the functionals {\mathcal F}_k^{\alpha_ k, \beta _ k} \colon L^1(\Omega)\longrightarrow [0, +\infty] defined as

    \begin{equation} {\mathcal F}_k^{\alpha_ k,\beta _ k}(u): = \begin{cases} \int_{ \Omega_k} f\Big(\frac{x}{ \varepsilon_k},\nabla u\Big) \,dx + \alpha_ k\int_{ \Omega\setminus \Omega_k } f\Big(\frac{x}{ \varepsilon_k},\nabla u\Big) \,dx \\[10pt] \quad\quad\quad\quad\quad\quad +\int_{S_u \cap \Omega_k} g\Big(\frac{x}{ \varepsilon_k},\nu_u\Big)\,d {\mathcal H}^{n-1} + \beta _ k\int_{S_u \cap ( \Omega\setminus \Omega_k)} g\Big(\frac{x}{ \varepsilon_k},\nu_u\Big)\,d {\mathcal H}^{n-1} \\[8pt] \quad\quad\quad\quad\quad\quad\quad \quad\quad\quad\quad\quad\quad \quad\quad\quad\quad\quad\quad\quad\quad\quad \text{if}\; u\in SBV^p( \Omega) \cr \cr +\infty \quad\quad\quad\quad\quad\quad\quad \quad\quad\quad\quad\quad\quad \quad\quad\quad\quad\quad\quad\quad\quad \text{otherwise,} \end{cases} \end{equation} (3.1)

    where f and g are Q -periodic in the first variable and satisfy (H1)–(H2) and (H3)–(H4), respectively.

    With the help of some specific examples, which correspond to some specific choices of f , g , and \Omega_k , we show that the \Gamma -limit of {\mathcal F}_k^{\alpha_ k, \beta _ k} is highly sensitive both to the choice of the coefficients \alpha_ k , \beta _ k and to the asymptotic behaviour of \alpha_ k , \beta _ k compared to the period of the microstructure \varepsilon_k . The examples we are going to discuss are taken from Barchiesi and Focardi [5] (see also Cagnetti and Scardia [14] and Focardi, Gelli and Ponsiglione [17]), from Barchiesi, Lazzaroni and Zeppieri [6], and from Pellet, Scardia and Zeppieri [19]. For the corresponding proofs we refer the reader to the aforementioned papers.

    In this subsection we consider the limit case \alpha_ k = \beta _ k = 0 ; i.e., we consider the functionals

    \begin{equation} {\mathcal F}_k^{0,0}(u): = \begin{cases} \int_{ \Omega_k} f\Big(\frac{x}{ \varepsilon_k},\nabla u\Big) \,dx +\int_{S_u \cap \Omega_k} g\Big(\frac{x}{ \varepsilon_k},\nu_u\Big)\,d {\mathcal H}^{n-1} & \text{if}\; u\in SBV^p( \Omega), \cr \cr +\infty & \text{otherwise in $L^1( \Omega)$.} \end{cases} \end{equation} (3.2)

    Loosely speaking, in this case the soft or weak inclusions in the material are replaced by perforations [5,14,17].

    Theorem 3.1 (Homogenisation of periodic porous brittle materials). Let {\mathcal F}_k^{0, 0} be the functionals as in (3.2). Then ({\mathcal F}_k^{0, 0}) \Gamma -converges both with respect to the convergence in Definition 2.2 and with respect to the L^1(\Omega) -convergence to the functional {\mathcal F}^0 which is finite on GSBV^p(\Omega) and given by

    \begin{equation} {\mathcal F}^0(u) = \int_ \Omega f^{0}(\nabla u) \,dx +\int_{S_u}g^0(\nu_u)\, d\mathcal H^{n-1}, \end{equation} (3.3)

    where f^0 and g^0 are, respectively, given by the following homogenisation formulas

    \begin{equation} f^0(\xi) = \inf\left\{\int_{Q \cap E}f(y,\nabla u) \,dx \colon u\in W^{1,p}(Q\cap E), \; u = u_\xi \; \mathit{\text{near}}\; \partial Q \right\}, \end{equation} (3.4)

    for every \xi \in {\mathbb R}^n , whereas

    \begin{multline} g^0(\nu) = \lim\limits_{T \to +\infty} \frac{1}{T^{n-1}}\inf\bigg\{\int_{S_u \cap TQ^\nu \cap E}g(y,\nu_u)\,d\mathcal H^{n-1} \colon u\in \mathcal P(TQ^\nu\cap E), \; u = u^{\nu,1}_{0} \; \mathit{\text{near}}\; \partial TQ^\nu \bigg\}, \end{multline} (3.5)

    for every \nu \in {\mathbb S}^{n-1} .

    Proof. Theorem 2.16 and Corollary 2.17 yield the existence of a subsequence k_j \to +\infty such that the corresponding functionals {\mathcal F}_{k_j}^{0, 0} \Gamma -converge to {\mathcal F} as in (2.27), both with respect to the convergence in Definition 2.2 and to the L^1(\Omega) -convergence. Then, the homogenisation formulas (3.4) and (3.5) together with the identity {\mathcal F} = {\mathcal F}^0 follow from [5,Theorem 4]. Finally, since (3.4) and (3.5) are subsequence-independent, invoking the Urysohn property [15,Proposition 8.3] readily implies the \Gamma -convergence of the whole sequence ({\mathcal F}_k^{0, 0}) to {\mathcal F}^0 .

    The following result is an immediate consequence of Theorem 3.1 and of an adaptation of the Cagnetti and Scardia extension result [14,Theorem 1.3] to the case of a general exponent p > 1 .

    Corollary 3.2. Let \alpha_ k, \beta _ k \to 0 and let {\mathcal F}_k^{\alpha_ k, \beta _ k} be the corresponding functionals as in (2.3). Then, the sequence ({\mathcal F}_k^{\alpha_ k, \beta _ k}) \Gamma -converges to the functional {\mathcal F}^0 given by (3.3).

    Proof. By Theorem 2.16 (up to subsequences not relabelled) the functionals {\mathcal F}^{\alpha_ k, \beta _ k}_k \Gamma -converge to {\mathcal F} . Since {\mathcal F}_k^{0, 0} \leq {\mathcal F}_k^{\alpha_ k, \beta _ k} by Theorem 3.1 we immediately get {\mathcal F}^0\leq {\mathcal F} .

    We now prove the converse inequality. To this end let u\in SBV^p(\Omega) \cap L^\infty(\Omega) and let (u_k) \subset SBV^p(\Omega_k) be a recovery sequence for {\mathcal F}_k^{0, 0} . That is u_k \leadsto u and \lim_k {\mathcal F}_k^{0, 0}(u_k) = {\mathcal F}(u) . Since the functionals {\mathcal F}_k^{0, 0} decrease by truncation it is not restrictive to assume that \|u_k\|_{L^\infty(\Omega)}\leq \|u\|_{L^\infty(\Omega)} . Starting from (u_k) we now want to construct a sequence (v_k) which both satisfies v_k \to u in L^1(\Omega) and \lim_k {\mathcal F}_k^{\alpha_ k, \beta _ k}(v_k) = {\mathcal F}^0(u) . To this end, we start noticing that the bounds (2.1) and (2.2) readily imply

    \begin{equation} \sup\limits_k \Big(\int_{ \Omega_k}|\nabla u_k|^p \,dx+\mathcal H^{n-1}(S_{u_k}\cap \Omega_k)\Big) \lt +\infty. \end{equation} (3.6)

    For every fixed k let v_k: = T^k u_k \in SBV^p(\Omega) be the extended function of u_k to \Omega whose existence is given by [14,Theorem 1.3]; i.e., v_k is such that v_k = u_k a.e. in \Omega_k , \|u_k\|_{L^\infty(\Omega)}\leq \|u\|_{L^\infty(\Omega)} , and

    \begin{equation} \int_{ \Omega}|\nabla v_k|^p \,dx+\mathcal H^{n-1}(S_{v_k}\cap \Omega)\leq C\,\Big(\int_{ \Omega_k}|\nabla u_k|^p \,dx+\mathcal H^{n-1}(S_{u_k}\cap \Omega_k)\Big) \end{equation} (3.7)

    for some C > 0 independent of k . By definition of v_k , also invoking the Ambrosio compactness Theorem, it is immediate to check that v_k \to u in L^1(\Omega) .

    By (3.7) we get that

    \begin{eqnarray} & \alpha_ k\int_{ \Omega\setminus \Omega_k}|\nabla v_k|^p \,dx+\beta _ k\mathcal H^{n-1}(S_{v_k}\cap ( \Omega\setminus \Omega_k)) \leq \alpha_ k \int_{ \Omega}|\nabla v_k|^p \,dx+ \alpha_ k \mathcal H^{n-1}(S_{v_k}\cap \Omega)& \\ & \leq C\,\max\{\alpha_ k,\beta _ k\}\, \Big(\int_{ \Omega_k}|\nabla u_k|^p \,dx+\mathcal H^{n-1}(S_{u_k}\cap \Omega_k)\Big),& \end{eqnarray} (3.8)

    where (3.8) is infinitesimal thanks to (3.6), since \max\{\alpha_ k, \beta _ k\}\to 0 as k\to +\infty . Thus eventually

    \lim\limits_{k \to +\infty} {\mathcal F}_k^{\alpha_ k,\beta _ k}(v_k) = \lim\limits_{k \to +\infty} {\mathcal F}_k^{0,0}(u_k) = {\mathcal F}^0(u),

    hence (v_k) is the desired sequence. Therefore, by the \Gamma -convergence of {\mathcal F}_k^{\alpha_ k, \beta _ k} to {\mathcal F} we can deduce that for every u\in SBV^p(\Omega)\cap L^\infty(\Omega) it holds {\mathcal F}(u)\leq {\mathcal F}^0(u) .

    Now let u\in GSBV^p(\Omega) and denote with u^m its truncated function at level m > 0 . We clearly have {\mathcal F}(u^m)\leq {\mathcal F}^0(u^m)\leq {\mathcal F}^0(u) , hence the desired inequality follows by the L^1(\Omega) convergence u^m to u and by the lower semicontinuity of {\mathcal F} .

    The following remarks are in order.

    Remark 3.3. In view of Remark 2.10, both in Theorem 3.1 and in Corollary 3.2 assumption (H2) on f can be dropped.

    Remark 3.4 (On f^0 ). The homogenised volume energy density f^0 given by (3.4) is the same as that obtained by Acerbi, Chiadó-Piat, Dal Maso and Percivale [1] in the case of elastic perforated materials. Moreover, it is easy to check that if f is p -homogeneous then the corresponding f^0 given by (3.4) is also p -homogeneous.

    Remark 3.5 (Energy decoupling). In spite of the strong degeneracy of the integrands in (3.2) (resp. in (3.1)), which in this case are identically equal to zero (resp. both infinitesimal) in the \varepsilon_k -periodic set \Omega\setminus \Omega_k , Theorem 3.1 (resp. Corollary 3.2) shows that the functionals {\mathcal F}_k^{0, 0} (resp. {\mathcal F}_k^{\alpha_ k, \beta _ k} ) exhibit a limit behaviour which is qualitatively similar to that of free-discontinuity functionals with coercive integrands [9,13,18]. Namely, in the homogenised limit there is no interaction between bulk and surface term. As a consequence the homogenised surface energy density g^0 does not depend on t , and therefore the \Gamma -limit is of brittle type.

    In this section we show that if only one of the coefficients \alpha_ k and \beta _ k is infinitesimal (while the other stays uniformly bounded from below), then the asymptotic behaviour of the functionals {\mathcal F}^{\alpha_ k, \beta _ k}_k can be very different from that of {\mathcal F}^{0, 0}_k (or of {\mathcal F}^{\alpha_ k\to 0, \beta _ k\to 0}_k ). In particular, we show that in this case a volume-surface energy coupling cannot be excluded in general. To do so we exhibit coefficients \alpha_ k, \beta _ k , integrands f, g and a geometry for the periodic set E which give rise to the desired limit coupling. This is done by resorting to the analysis of Barchiesi, Lazzaroni and Zeppieri [6] and Pellet, Scardia and Zeppieri [19], which is briefly reviewed in Subsection 3.2.1 and Subsection 3.2.2, respectively.

    The functionals analysed in [6] and [19] are both of Mumford-Shah type and can be written in the form

    \begin{equation} \mathcal{MS}^{\alpha_k,\beta_k}_k(u) = \int_\Omega a_k\Big(\frac{x}{ \varepsilon_k}\Big)|\nabla u|^2\,dx+\int_{S_u}b_k\Big(\frac{x}{ \varepsilon_k}\Big)\,d\mathcal H^{n-1}, \quad u\in SBV^2( \Omega) \end{equation} (3.9)

    where a_k, b_k \colon {\mathbb R}^n \to [0, 1] are Q -periodic functions and in the periodicity cell Q are defined as

    \begin{equation} a_k(y) = \begin{cases} \alpha_ k & \text{if}\; y\in \overline Q_r\cr 1 & \text{if}\; y\in Q \setminus \overline Q_r \end{cases} \qquad b_k(y) = \begin{cases} \beta _ k & \text{if}\; y\in \overline Q_r\cr 1 & \text{if}\; y\in Q \setminus \overline Q_r \end{cases} \end{equation} (3.10)

    with r\in (0, 1) . From (3.9)–(3.10) we infer that in this case f = f(\xi) = |\xi|^2 , g\equiv 1 , and \Omega_k = \Omega \cap \varepsilon_k E with E = {\mathbb R}^n \setminus \bigcup_{i\in {\mathbb Z}^n} (\overline Q_r+i) .

    Remark 3.6 (Mumford-Shah functional in perforated domains). The choice \alpha_ k = \beta _ k = 0 corresponds to the Mumford-Shah functionals in perforated domain. The functional \mathcal{MS}^{0, 0}_k is a special instance of (3.2) and its homogenised limit is treated in [14,17] for general sets E . In this case the homogenised integrands (3.4) and (3.5) reduce, respectively, to

    \begin{equation} f^0(\xi) = \inf\left\{\int_{Q \setminus \overline Q_r}|\nabla u|^2 \,dx \colon u\in W^{1,2}(Q\setminus \overline Q_r), \; u = u_\xi \; \text{ near }\; \partial Q \right\}, \end{equation} (3.11)

    for every \xi \in {\mathbb R}^n , and to

    \begin{equation} g^0(\nu) = \lim\limits_{T \to +\infty} \frac{1}{T^{n-1}}\inf\left\{\mathcal H^{n-1}(S_u \cap TQ^\nu \cap E) \colon u\in \mathcal P(TQ^\nu\cap E), \; u = u^{\nu,1}_{0} \; \text{ near }\; \partial TQ^\nu \right\}, \end{equation} (3.12)

    for every \nu \in {\mathbb S}^{n-1} . From (3.11) and (3.12) it is easy to check that f^0(\xi) = A^0 \xi \cdot \xi , for some A^0 \in {\mathbb R}^{n\times n} which satisfies \tilde c_1 I \leq A^0 \leq I , in the sense of quadratic forms (cf. (2.28)). Hence, f^0 is a positive quadratic form. Moreover, it holds g^0(e_i) = 1-r^{n-1} , for every i = 1, \ldots, n .

    We consider the case \alpha_ k \to 0 and \beta _ k = 1 which models the situation where the periodic set \Omega \setminus \Omega_k is occupied by a brittle material with a very small elastic modulus. For this reason, we refer to the set \Omega \setminus \Omega_k as the set of soft inclusions. With this choice the functionals in (3.9) become

    \begin{equation} \mathcal{MS}^{\alpha_k,1}_k(u) = \int_{\Omega_k}|\nabla u|^2\,dx+\alpha_ k\int_{ \Omega\setminus \Omega_k}|\nabla u|^2\,dx+\mathcal H^{n-1}(S_u), \quad u\in SBV^2( \Omega). \end{equation} (3.13)

    In [6] Barchiesi, Lazzaroni and Zeppieri showed that the asymptotic behaviour of \mathcal{MS}^{\alpha_k, 1}_k heavily depends on the mutual vanishing rate of \alpha_ k and \varepsilon_k ; that is, it depends on the parameter

    \begin{equation} \ell: = \lim\limits_{k \to +\infty}\frac{\alpha_ k}{ \varepsilon_k} \in [0,+\infty]. \end{equation} (3.14)

    For the proof of Theorem 3.7 below we refer the reader to [6,Theorems 1,4,and Remark 6].

    Theorem 3.7 (Homogenisation of periodic brittle materials with soft inclusions). Let \mathcal{MS}_k^{\alpha_ k, 1} be the functionals defined in (3.13) and let \ell\in [0, +\infty] be as in (3.14). Then, up to subsequences not relabelled, (\mathcal{MS}_k^{\alpha_ k, 1}) \Gamma -converges to the functional {\mathcal F}^\ell which is finite on GSBV^2(\Omega) and given by

    \begin{equation} {\mathcal F}^\ell(u) = \int_ \Omega f^0(\nabla u) \,dx +\int_{S_u}g^\ell([u],\nu_u)\, d\mathcal H^{n-1}, \end{equation} (3.15)

    where f^0 is as in (3.11) and for every t\in {\mathbb R} , \nu\in {\mathbb S}^{n-1}

    g^\ell(t,\nu) = \begin{cases} g^0 (\nu) & \mathit{\text{if}}\; \ell = 0 \cr 1 & \mathit{\text{if}}\; \ell = +\infty. \end{cases}

    Moreover for every \ell \in (0, +\infty) it holds

    \begin{equation} \min\bigl\{g_0(e_i)+c_\ell t^2, \, 1\bigr\}\leq g^\ell(t,e_i) \leq \min\bigl\{g_0(e_i)+\hat c_\ell t, \, 1\bigr\} \end{equation} (3.16)

    for every t > 0 , i = 1, \dots, n , and for some c_\ell, \hat c_\ell > 0 , with \lim_{\ell \to 0+}c_\ell = \lim_{\ell \to 0+}\hat c_\ell = 0 .

    Remark 3.8. The following remarks are in order.

    (i) As far as the homogenised volume energy is concerned, the soft inclusions are (energetically) equivalent to the perforations in the material.

    (ii) For \ell = 0 , which corresponds to \alpha_ k \ll \varepsilon_k , the functionals \mathcal{MS}_k^{\alpha_ k, 1} are equivalent to the functionals \mathcal{MS}_k^{0, 0} , in the sense of \Gamma -convergence.

    (iii) For \ell\in (0, +\infty) the bounds in (3.16) imply that, along the coordinate directions, g^\ell depends on t . Moreover it becomes constant (and equal to 1 ) above a certain threshold t_0 > 0 ; i.e., g^\ell is of cohesive type. Being the microscopic energies \mathcal{MS}_k^{\alpha_ k, 1} of brittle type, the cohesive behaviour of g^\ell can only be explained as the result of a non trivial bulk-surface coupling by homogenisation. This interaction is particularly apparent from the upper-bound construction in [6] which we briefly illustrate here in the case n = 2 .

    For i = 1, 2 we have g^\ell(t, e_i) = {\mathcal F}^\ell(u_{0}^{e_i, t}, Q) , moreover it is immediate to check that g^\ell(t, e_1) = g^\ell(t, e_2) . Clearly g^\ell(t, e_2)\leq 1 for every t > 0 . Then, to get the upper bound in (3.16) it suffices to show that g^\ell(t, e_2)\leq g^0(e_2)+\hat c_\ell t for some \hat c_\ell > 0 . Let R\subset Q \subset {\mathbb R}^2 be the open rectangle defined as

    \begin{equation*} R: = (-\tfrac r2,\tfrac r2)\times(-\tfrac{\tau}{2},\tfrac{\tau}{{2}}), \end{equation*}

    with \tau \in (0, r) to be determined. Set

    R_k: = Q\cap \bigcup\limits_{i\in {\mathbb Z}}\Big( \varepsilon_k R+(i \varepsilon_k,0)\Big)

    and let (u_k)\subset SBV^2(Q) be the sequence of functions defined as

    \begin{equation*} u_k(x): = \begin{cases} t & \rm{if } x\in Q\setminus R_k \text{ and } x_2\geq0,\cr \frac{t}{2}+\frac{t}{\tau \varepsilon_k}\,x_n & \rm{if } x\in R_k,\cr 0 & \rm{if } x\in Q\setminus R_k \text{ and } x_2 \lt 0, \end{cases} \end{equation*}

    (see Figure 2). We clearly have u_k\rightarrow u_{0}^{e_2, t} in L^1(Q) ; moreover

    \begin{equation*} \int_{R_k}|\nabla u_k|^2 \,dx\leq \Big(\frac1{ \varepsilon_k}+1\Big)\frac{t^2}{\tau}\quad \text{and}\quad \mathcal H^1(S_{u_k})\leq \varepsilon_k\Big(\frac{1}{ \varepsilon_k}+1\Big)\bigl(1-r+2\tau\bigr), \end{equation*}
    Figure 2.  Construction of the recovery sequence in an \varepsilon_k -cell, across the interface x_2 = 0 .

    therefore

    \begin{equation*} g^\ell(t,e_2) = {\mathcal F}^\ell(u^{e_2,t}_{0},Q) \leq \limsup\limits_{k\to +\infty} \mathcal{MS}^{\alpha_ k,1}_k(u_k,Q). \leq 1-r+2\tau +\ell \frac{t^2}{\tau} \end{equation*}

    Hence, by optimising on \tau we get

    \begin{equation} g^\ell(t,e_2)\leq 1-r+2\sqrt{2 \ell}\,t \end{equation} (3.17)

    thus the desired estimate follows with \hat c_\ell = 2\sqrt{2 \ell} , by recalling that g^0(e_2) = 1-r .

    Loosely speaking, the construction as above shows that, the cost of an elastic deformation of the soft inclusions is of the same order of the energy spent to create a microscopic crack. Since the former depends linearly on t (while the latter is constant in t ) for small values of t , to approximate a macroscopic crack it can be convenient to combine microscopic deformations of the soft inclusions (with high gradients) and microscopic jumps.

    (iv) Even if not immediately apparent from the homogenisation formulas, a volume-surface interaction takes place for \ell = 0 , as well. Indeed, in this case g^\ell = g^0 whereas in \mathcal{MS}^{\alpha_ k, 1}_k the surface energy density is identically equal to one. In this case in fact, the cost of an elastic deformation of the soft inclusions is negligible (cf. (3.17) for \ell = 0 ) so that to approximate a macroscopic crack it is never convenient to introduce microscopic cracks inside the soft material. On the contrary, in the regime \ell = +\infty , which corresponds to \alpha_ k\gg \varepsilon_k , there is a complete volume-surface decoupling, as in the coercive case.

    We consider the case \alpha_ k = 1 and \alpha_ k \to 0 which models the situation where the periodic set \Omega \setminus \Omega_k is occupied by a brittle material with a very small fracture resistance. For this reason, we refer to the set \Omega \setminus \Omega_k as the set of weak inclusions. With this choice the functionals in (3.9) become

    \begin{equation} \mathcal{MS}^{1,\beta _ k}_k(u) = \int_{\Omega}|\nabla u|^2\,dx+ \mathcal H^{n-1}(S_u\cap \Omega_k)+ \beta _ k \mathcal H^{n-1}(S_u\cap ( \Omega\setminus \Omega_k)), \quad u\in SBV^2( \Omega). \end{equation} (3.18)

    In [19] Pellet, Scardia and Zeppieri showed that the asymptotic behaviour of \mathcal{MS}^{1, \beta _ k}_k heavily depends on the mutual vanishing rate of \beta _ k and \varepsilon_k , that is on the parameter

    \begin{equation} \ell': = \lim\limits_{k}\frac{\beta _ k}{ \varepsilon_k} \in [0,+\infty]. \end{equation} (3.19)

    For the proof of Theorem 3.9 below we refer the reader to [19].

    Theorem 3.9 (Homogenisation of periodic brittle materials with weak inclusions). Let \mathcal{MS}_k^{1, \beta _ k} be the functionals defined in (3.18) and let \ell'\in[0, +\infty] be as in (3.19). Then, up to subsequences not relabelled, (\mathcal{MS}_k^{1, \beta _ k}) \Gamma -converges to the functional {\mathcal F}^{\ell'} which is finite on GSBV^2(\Omega) and given by

    \begin{equation} {\mathcal F}^{\ell'}(u) = \int_ \Omega f^{\ell'}(\nabla u) \,dx +\int_{S_u}g^0(\nu_u)\, d\mathcal H^{n-1}, \end{equation} (3.20)

    where g^0 is as in (3.12) and for every \xi\in {\mathbb R}^n

    f^{\ell'}(\xi) = \begin{cases} f^0 (\xi) & \mathit{\text{if}}\; \ell' = 0 \cr |\xi|^2 & \mathit{\text{if}}\; \ell' = +\infty. \end{cases}

    Moreover for every \ell' \in (0, +\infty) it holds

    \begin{equation} f^0(\xi)\leq f^{\ell'}(\xi) \leq \min\bigl\{|\xi|^2, f^0(\xi)+C \ell' \bigr\} \end{equation} (3.21)

    for every \xi\in {\mathbb R}^n and for some C > 0 .

    Remark 3.10 The following remarks are in order.

    (i) As far as the homogenised surface energy is concerned, the weak inclusions are (energetically) equivalent to the perforations in the material.

    (ii) For \ell' = 0 , which corresponds to \beta _ k \ll \varepsilon_k , the functionals \mathcal{MS}_k^{1, \beta _ k} are equivalent to the functionals \mathcal{MS}_k^{0, 0} , in the sense of \Gamma -convergence. Indeed, "removing the weak inclusions from the material" has an infinitesimal cost of order \beta _ k/ \varepsilon_k given by the perimeter of the weak inclusions (proportional to \beta _ k \varepsilon_k^{n-1} ) multiplied by \varepsilon_k^{-n} (the number of \varepsilon_k -cells contained in \Omega ). In this case a volume-surface energy coupling takes place since the elastic energy can be lowered by introducing cracks in the materials.

    (iii) For \ell'\in (0, +\infty) the bounds in (3.21) hold true (see [19,Lemma 6.1]). The bound from below is immediate and it is a consequence of the trivial bound \mathcal{MS}_k^{0, 0}\leq \mathcal{MS}_k^{1, \beta _ k} . The bound from above shows that for large deformations; i.e., for large |\xi| , to approximate a macroscopic elastic deformation is energetically favourable to mix elastic deformations and jumps in the weak inclusions. Moreover, (3.21) implies that for |\xi| large it holds f^{\ell'}(\xi) < |\xi|^2 . The latter shows that a stiffness degradation occurs in the homogenised limit, and that the macroscopic energy {\mathcal F}^{\ell'} describes a damaged material (the same being true for \ell = 0 ).

    (iv) The bounds in (3.21) combined with an easy scaling argument show that in the regime \ell'\in(0, +\infty) the homogenised volume energy density f^{\ell'} is not 2 -homogeneous. Indeed, assume by contradiction that this is not the case and let \lambda \neq 0 . Taking into account that f^0 is 2 -homogeneous (see Remark 3.4), we can replace in (3.21) \xi with \lambda \xi and divide by \lambda^2 to get

    f^0 (\xi)\leq f^{\ell'}(\xi)\leq\min\bigg\{|\xi|^2,f^0(\xi)+\frac{C\ell'}{\lambda^2}\bigg\}.

    Therefore by letting |\lambda| \to +\infty we get f^{\ell'}\equiv f^0 which leads to a contradiction in view of [19,Proposition 6.10].

    (v) In the regime \ell' = +\infty , which corresponds to \beta _ k\gg \varepsilon_k , there is a complete volume-surface decoupling, as in the coercive case. Loosely speaking, in this case the fracture resistance of the weak inclusions is not small enough to make cracks energetically more convenient than (or at least comparable to) elastic deformations.

    This work was supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under the Germany's Excellence Strategy EXC 2044-390685587, Mathematics Münster: Dynamics–Geometry–Structure.

    The author declares no conflict of interest.



    [1] Acerbi E, Chiadó Piat V, Dal Maso G, et al. (1992) An extension theorem from connected sets, and homogenization in general periodic domains. Nonlinear Anal Theor 18: 481-496. doi: 10.1016/0362-546X(92)90015-7
    [2] Ambrosio L, Fusco N, Pallara D (2000) Functions of Bounded Variations and Free Discontinuity Problems. Oxford: Clarendon Press.
    [3] Barchiesi M (2018) Toughening by crack deflection in the homogenization of brittle composites with soft inclusions. Arch Ration Mech Anal 227: 749-766. doi: 10.1007/s00205-017-1173-5
    [4] Barchiesi M, Dal Maso G (2009) Homogenization of fiber reinforced brittle materials: The extremal cases. SIAM J Math Anal 41: 1874-1889. doi: 10.1137/080744372
    [5] Barchiesi M, Focardi M (2011) Homogenization of the Neumann problem in perforated domains: An alternative approach. Calc Var Partial Dif 42: 257-288. doi: 10.1007/s00526-010-0387-2
    [6] Barchiesi M, Lazzaroni G, Zeppieri CI (2016) A bridging mechanism in the homogenization of brittle composites with soft inclusions. SIAM J Math Anal 48: 1178-1209. doi: 10.1137/15M1007343
    [7] Bouchitté G, Fonseca I, Leoni G, et al. (2002) A global method for relaxation in W1,p and in SBVp. Arch Ration Mech Anal 165: 187-242. doi: 10.1007/s00205-002-0220-y
    [8] Braides A, Defranceschi A (1998) Homogenization of Multiple Integrals. New York: Oxford University Press.
    [9] Braides A, Defranceschi A, Vitali E (1996) Homogenization of free discontinuity problems. Arch Ration Mech Anal 135: 297-356. doi: 10.1007/BF02198476
    [10] Braides A, Garroni A (1995) Homogenization of periodic nonlinear media with stiff and soft inclusions. Math Mod Meth Appl Sci 5: 543-564. doi: 10.1142/S0218202595000322
    [11] Braides A, Solci M (2013) Multi-scale free-discontinuity problems with soft inclusions. Boll Unione Mat Ital 1: 29-51.
    [12] Cagnetti F, Dal Maso G, Scardia L, et al. (2019) Γ-convergence of free-discontinuity problems. Ann Inst H Poincaré Anal Non Linéaire 36: 1035-1079. doi: 10.1016/j.anihpc.2018.11.003
    [13] Cagnetti F, Dal Maso G, Scardia L, et al. (2019) Homogenisation of stochastic free-discontinuity problems. Arch Ration Mech Anal 233: 935-974. doi: 10.1007/s00205-019-01372-x
    [14] Cagnetti F, Scardia L (2011) An extension theorem in SBV and an application to the homogenization of the Mumford-Shah functional in perforated domains. J Math Pur Appl 95: 349-381. doi: 10.1016/j.matpur.2010.03.002
    [15] Dal Maso G (1993) An Introduction to Γ-Convergence. Boston: Birkhäuser.
    [16] Dal Maso G, Zeppieri CI (2010) Homogenization of fiber reinforced brittle materials: The intermediate case. Adv Calc Var 3: 345-370.
    [17] Focardi M, Gelli MS, Ponsiglione M (2009) Fracture mechanics in perforated domains: A variational model for brittle porous media. Math Mod Meth Appl Sci 19: 2065-2100. doi: 10.1142/S0218202509004042
    [18] Giacomini A, Ponsiglione M (2006) A Γ-convergence approach to stability of unilateral minimality properties. Arch Ration Mech Anal 180: 399-447. doi: 10.1007/s00205-005-0392-3
    [19] Pellet X, Scardia L, Zeppieri CI (2019) Homogenization of high-contrast Mumford-Shah energies. SIAM J Math Anal 51: 1696-1729. doi: 10.1137/18M1189804
    [20] Scardia L (2008) Damage as Γ-limit of microfractures in anti-plane linearized elasticity. Math Mod Meth Appl Sci 18: 1703-1740. doi: 10.1142/S0218202508003170
    [21] Scardia L (2010) Damage as the Γ-limit of microfractures in linearized elasticity under the non-interpenetration constraint. Adv Calc Var 3: 423-458.
  • This article has been cited by:

    1. L. De Luca, M. Ponsiglione, Variational models in elasticity, 2021, 3, 2640-3501, 1, 10.3934/mine.2021015
    2. A.G. Kolpakov, S.I. Rakin, Local stresses in the reinforced plate with orthogonal systems of fibers. The homogenized strength criterion, 2023, 275, 00457949, 106922, 10.1016/j.compstruc.2022.106922
    3. Xavier Pellet, Lucia Scardia, Caterina Ida Zeppieri, Stochastic homogenisation of free-discontinuity functionals in randomly perforated domains, 2024, 17, 1864-8258, 643, 10.1515/acv-2022-0052
  • Reader Comments
  • © 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3478) PDF downloads(375) Cited by(3)

Figures and Tables

Figures(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog