Indonesian rubber farming has the largest area in the world, but its implementation faces various risks that decrease productivity and farm income. This study is designed to specify the risk perception, risk attitude and determinant factors for smallholder rubber farmers. The research location was in four subdistricts in West Kalimantan Province, with a sample size of 200 farmers. Data collection was carried out by interview using a structured questionnaire. The risk matrix, Holt and Laury's method and the logit model were used to identify risk perception, risk attitude and determinant factors. The study results showed that most rubber farmers were risk-averse and perceived climate change, plant diseases and price change as high risks. The logit model found that farmers' age, education, rubber plantation size, rubber age, distance and use of rubber clones had a positive and significant effect on farmers' risk perception, while the family size and farming experience had a negative effect. Regarding risk attitude, the logit model found that rubber age, distance and risk perception of price change had a positive and significant effect on farmers' risk aversion, while farmers' age and use of rubber clones had a negative effect. This study recommends providing informal education to the farmers through training and counseling, encouraging the farmers to replant old or damaged rubber trees and adopt rubber clones. Furthermore, it is also necessary to improve road facilities and infrastructure, communication and transportation access to facilitate farming activities.
Citation: Imelda, Jangkung Handoyo Mulyo, Any Suryantini, Masyhuri. Understanding farmers' risk perception and attitude: A case study of rubber farming in West Kalimantan, Indonesia[J]. AIMS Agriculture and Food, 2023, 8(1): 164-186. doi: 10.3934/agrfood.2023009
[1] | Filippo Gazzola, Elsa M. Marchini . The moon lander optimal control problem revisited. Mathematics in Engineering, 2021, 3(5): 1-14. doi: 10.3934/mine.2021040 |
[2] | Aleksandr Dzhugan, Fausto Ferrari . Domain variation solutions for degenerate two phase free boundary problems. Mathematics in Engineering, 2021, 3(6): 1-29. doi: 10.3934/mine.2021043 |
[3] | Valentina Volpini, Lorenzo Bardella . Asymptotic analysis of compression sensing in ionic polymer metal composites: The role of interphase regions with variable properties. Mathematics in Engineering, 2021, 3(2): 1-31. doi: 10.3934/mine.2021014 |
[4] | Hyeonbae Kang, Shigeru Sakaguchi . A symmetry theorem in two-phase heat conductors. Mathematics in Engineering, 2023, 5(3): 1-7. doi: 10.3934/mine.2023061 |
[5] | Lucio Boccardo, Giuseppa Rita Cirmi . Regularizing effect in some Mingione’s double phase problems with very singular data. Mathematics in Engineering, 2023, 5(3): 1-15. doi: 10.3934/mine.2023069 |
[6] | Stefano Biagi, Serena Dipierro, Enrico Valdinoci, Eugenio Vecchi . A Hong-Krahn-Szegö inequality for mixed local and nonlocal operators. Mathematics in Engineering, 2023, 5(1): 1-25. doi: 10.3934/mine.2023014 |
[7] | Sandro Salsa, Francesco Tulone, Gianmaria Verzini . Existence of viscosity solutions to two-phase problems for fully nonlinear equations with distributed sources. Mathematics in Engineering, 2019, 1(1): 147-173. doi: 10.3934/Mine.2018.1.147 |
[8] | Guido De Philippis, Filip Rindler . Fine properties of functions of bounded deformation-an approach via linear PDEs. Mathematics in Engineering, 2020, 2(3): 386-422. doi: 10.3934/mine.2020018 |
[9] | Cristiana De Filippis . Optimal gradient estimates for multi-phase integrals. Mathematics in Engineering, 2022, 4(5): 1-36. doi: 10.3934/mine.2022043 |
[10] | Filippo Gazzola, Gianmarco Sperone . Remarks on radial symmetry and monotonicity for solutions of semilinear higher order elliptic equations. Mathematics in Engineering, 2022, 4(5): 1-24. doi: 10.3934/mine.2022040 |
Indonesian rubber farming has the largest area in the world, but its implementation faces various risks that decrease productivity and farm income. This study is designed to specify the risk perception, risk attitude and determinant factors for smallholder rubber farmers. The research location was in four subdistricts in West Kalimantan Province, with a sample size of 200 farmers. Data collection was carried out by interview using a structured questionnaire. The risk matrix, Holt and Laury's method and the logit model were used to identify risk perception, risk attitude and determinant factors. The study results showed that most rubber farmers were risk-averse and perceived climate change, plant diseases and price change as high risks. The logit model found that farmers' age, education, rubber plantation size, rubber age, distance and use of rubber clones had a positive and significant effect on farmers' risk perception, while the family size and farming experience had a negative effect. Regarding risk attitude, the logit model found that rubber age, distance and risk perception of price change had a positive and significant effect on farmers' risk aversion, while farmers' age and use of rubber clones had a negative effect. This study recommends providing informal education to the farmers through training and counseling, encouraging the farmers to replant old or damaged rubber trees and adopt rubber clones. Furthermore, it is also necessary to improve road facilities and infrastructure, communication and transportation access to facilitate farming activities.
Scope of the paper In this paper, we aim at investigating several properties for a natural shape optimisation problem that arises in heterogeneous heat conduction: what is the optimal way to design the properties of a material in order to optimise its performance? This question has received a lot of attention from the mathematical community over the last decades [1,10,11,12,13,14,15,16,22,30] and our goal in this paper is to offer some complementary qualitative results. Mathematically, these problems are often dubbed two-phase problems and write, in their most general form, as follows: considering that the piece consists of a basic material, with conductivity α>0, we try to find the best location ω for the inclusion of another material having conductivity β>α. The resulting diffusive part of the equation under consideration writes
−∇⋅((α+(β−α)1ω))∇⋅). | (1.1) |
This diffusive part is supplemented with a source term, and can be considered in elliptic or parabolic models. We study some aspects of both cases in the present paper in the case of radial geometries. To state the generic type of question we are interested in, we write down the typical equation in the elliptic case: for a ball Ω, a source term f∈L2(Ω) (the influence of which is also discussed) and an inclusion ω⊂Ω, let uω be the unique solution of
{−∇⋅((α+(β−α)1ω)∇uω)=finΩ,uω=0on∂Ω. | (1.2) |
We consider a volume constraint, enforced by a parameter V1∈(0;Vol(Ω)), and we investigate the problem
supω⊂Ω,Vol(ω)=V1J(ω):=∫Ωj(uω), |
for a certain non-linearity j. More specifically, we consider the set
M(Ω):={a∈L∞(Ω):a=α+(β−α)1ωforsomemeasurableω⊂Ω,Vol(ω)=V1}, | (1.3) |
also called the set of bang-bang functions, as well as its natural compactification for the weak L∞−∗ topology,
A(Ω):={a∈L∞(Ω),α⩽a⩽β,∫Ωa=V0:=αVol(Ω)+(β−α)V1}. |
For a∈A(Ω), we define ua,f as the solution of (1.2) with α+(β−α)1ω replaced with a. We will be interested in two formulations: the initial (unrelaxed) one
supa∈M(Ω)∫Ωj(ua) |
as well as the relaxed one
supa∈A(Ω)∫Ωj(ua). |
It should be noted that we will also for some results have to optimise with respect to the source term f, but that the main difficulty usually lies in handling the term a. The two formulations of the problem have their interest, as it may be interesting to see when the two coincides. In other words, is a solution to the second problem a solution of the first one? Let us already underline several basic facts: first, as is customary in this type of optimisation problems (we detail the references later on and for the moment refer to [31]) we do not expect existence of solutions in all geometry, and the proper type of relaxation should rather be of the H-convergence type. Nevertheless, we offer some results about these two problems. Second, the type of problems we are considering are not energetic (in the sense that the criterion we aim at optimising can not a priori be derived from the natural energy associated with the PDE constraint). This leads to several difficulties, most notably in handling the adjoint of the optimisation problem and in the ensuing loss of natural convexity or concavity of the functional to optimise. Third, we distinguish between two types of results: the first type correspond to Talenti inequalities, where we rearrange both the coefficient a and the source term f. In the elliptic case, this follows from results of [6], and our contribution here is the application of these methods to the parabolic case. A second type of result, given in Theorem II, deals with a possible identification of the two formulations (i.e., if a solution to the relaxed problem exists then it is a solution of the unrelaxed one) in radial geometries, and we do not need for this second type of results to rearrange the source term f. This result is the main contribution of this article.
Informal statement of the results Our goal is thus threefold. For the sake of presentation we indicate to which case (i.e., optimisation with respect to a, f or both) each item corresponds. We write ua,f for the solution of the equation with diffusion a and source term f. We will need a comparison inequality provided in [6] and that we recall in Theorem A.
1) Existence and partial characterisation in radial geometries (optimisation with respect to a and f) This matter of existence and/or characterisation of optimal a in the case of radial geometries is the topic of the two first results. The Talenti inequality from [6] leads to a comparison principle, but leaves open the question of the existence of optimal shapes: if Ω is a centred ball, is there a radially symmetric solution a∗ to the optimisation problem in Ω of the form
a∗=α+(β−α)1ω∗ |
for some measurable subset ω∗⊂Ω? We prove in Theorem I that it is the case when the function j is convex and we also allow ourselves to rearrange the source term f. We use the ideas contained in [15] to do so and prove this theorem for the sake of completeness; we highlight that the main contribution here is to prove that the methods of [15] work for non-energetic functionals.
2) Weak bang-bang property under monotonicity assumption in radial geometries (optimisation with respect to a) The second result of the "elliptic problem" part, is the main result of this paper, Theorem II. In it, we give a weak bang-bang property that does not require convexity assumptions on the function j (and so no clear convexity on J). Namely, we prove that if, in a centred ball Ω, a solution a∗ to the optimisation problem exists, and if j is increasing, then this solution has to be of bang-bang type. It is notable that, in this theorem, we do not require the term f to be rearranged as well and that we can handle non-energetic problem. This is proved by introducing, for two-phase problems, an oscillatory method reminiscent of the ideas of [25].
3) Comparison results for parabolic models (optimisation with respect to time-dependent a and f) We provide, in Theorem III, a parabolic Talenti inequality. The proof is an adaptation of a result of [29], combined with the methods of [6].
Plan of the paper This paper is organised as follows: in Section 1.2 we present the models, the optimisation problems and give some elements about the Schwarz rearrangement. In Section 1.3 we state our main results. Section 1.4 contains the bibliographical references. The rest of the paper is devoted to the proofs of the main results. Finally, in the Conclusion, we state several open problems that we deem interesting.
Henceforth, Ω is a centred ball in IRd, and V0∈(0;Vol(Ω)) is a fixed parameter that serves as a volume constraint. As explained in the first paragraph, we are interested in both elliptic and parabolic models. This leads us to define two admissible classes: the first one, used for elliptic problems, is
A(Ω):={a∈L∞(Ω):α⩽a⩽βa.e.inΩ,∫Ωa=V0} | (1.4) |
while the second, defined for a certain time horizon T>0, is
A(Ω;T):={a∈L∞((0;T)×Ω):α⩽a⩽βa.e.in(0;T)×Ω,fora.e.t∈(0;T),∫Ωa(t,⋅)=V0} | (1.5) |
The set of admissible sources, on which we also place a volume constraint modelled via a constant F0∈(0;Vol(Ω)), is
F(Ω):={f∈L∞(Ω):0⩽f⩽1a.e.inΩ,∫Ωf=F0}. | (1.6) |
Similarly, we define, in the parabolic case,
F(Ω;T):={f∈L∞((0;T)×Ω):0⩽f⩽1a.e.in(0;T)×Ω,fora.e.t∈(0;T)∫Ωf(t,⋅)=F0}. | (1.7) |
Main equation in the elliptic case In the elliptic case, the main equation reads as follows: for any a∈A(Ω) and any f∈F(Ω), uell,a,f is the unique solution of
{−∇⋅(a∇uell,a,f)=finΩ,uell,a,f=0inΩ. | (1.8) |
The solution uell,a,f is the unique minimiser in W1,20(Ω) of the energy functional
Ea,f:W1,20(Ω)∋u↦12∫Ωa|∇u|2−∫Ωfu. | (1.9) |
Remark 1. Although for the classes A(Ω) and A(Ω;T) the lower bounds 0<α⩽a ensure coercivity of the associated energy, it may be asked whether the non-negativity constraint on the sources can be relaxed. It may be difficult, as we need in our proofs the following crucial fact: when Ω is the ball, when f and a are radially symmetric functions of A(Ω) and F(Ω) respectively, the solution ua,f is radially non-increasing in Ω. This may not be the case, for instance when f<0 close to the center of the ball. Thus we choose simplicity and assume f⩾0 almost everywhere.
In the elliptic case, the goal is to solve the following problem: let j∈C1(IR) be a given non-linearity, then the problem is
supa∈A(Ω),f∈F(Ω){Jell(a,f):=∫Ωj(uell,a,f)}⋅Pell,j |
In [6], a comparison result that we will use later on is proved; we recall it in Theorem A. This comparison result states roughly speaking, that if j is increasing, there exist two radially symmetric functions ˜a and f∗ such that Jell(a,f)⩽Jell(˜a,f), with f∗ still admissible; ˜a, however, may violate some constraints. Here, our main contribution is Theorem I, in which we prove it is possible to choose a radially symmetric ˜a that satisfies the constraints if we assume that j is convex and C2. This is done by adapting the methods of [15].
Second, in Theorem II, we are interested in the following alternative formulation: f∈F(Ω) being fixed, solve
supa∈A(Ω){Jell(a):=∫Ωj(uell,a,f)}.Pell,j,f |
We prove, using an oscillatory technique that, if a solution a∗ exists and if j is increasing, then we must have a∗∈M(Ω). We underline that this result does not require rearranging f.
Main equation in the parabolic case In the parabolic case, the main equation reads as follows: for any a∈A(Ω;T), any f∈F(Ω;T), uparab,a,f is the unique solution of
{∂uparab,a,f∂t−∇⋅(a∇uparab,a,f)=fin(0;T)×Ω,uparab,a,f=0on∂Ω,uparab,a,f(0,⋅)=0inΩ. | (1.10) |
The parabolic optimisation problem assumes the following form: for two given non-linearities j1 and j2 in C1(IR) we consider the optimisation problem
supa∈A(Ω;T),f∈F(Ω;T){Jparab(a,f):=∬(0;T)×Ωj1(uparab,a,f)+∫Ωj2(uparab,a,f(T))}.Pparab,j1,j2 |
The main result is Theorem III, in which a parabolic isoperimetric inequality (with respect to the coefficient a) is obtained: namely, it is better to have radially symmetric a and f.
In this section we recall the key points about the Schwarz rearrangement, which will be used constantly throughout this paper, and about the rearrangement of Alvino and Trombetti [6,7] that is crucial in dealing with two-phase isoperimetric problems.
Schwarz rearrangement: definitions, properties and order relations We refer to section 1.4 for further references, for instance for parabolic isoperimetric inequalities and for the time being we recall the basic definitions of the Schwarz rearrangement. We refer to [19,20,23] for a thorough introduction.
Definition 2 (Schwarz rearrangement of sets). For a given bounded connected open set Ω0, the Schwarz rearrangement Ω∗0 of Ω0 is the unique centred ball BΩ0=B(0;RΩ0) such that
Vol(BΩ0)=Vol(Ω0). | (1.11) |
For rearrangements of functions, we use the distribution function: for any p∈[1;+∞), for any function u∈Lp(Ω), u⩾0, its distribution function is
μu:IR+∋t↦Vol({u>t}). | (1.12) |
Definition 3 (Schwarz rearrangement of a function). For any function u∈Lp(Ω0),u⩾0, its Schwarz rearrangement is the unique radially symmetric function u∗∈Lp(Ω∗0) having the same distribution function as u. u# stands for the one-dimensional function such that u∗=u#(cd|⋅|d) where cd:=Vol(B(0;1)).
As a consequence of the equimeasurability of the function and of its rearrangement* there holds:
* Two functions are called equimeasurable if they have the same distribution functions.
∀p∈[1;+∞),∀u∈Lp(Ω0),u⩾0,∫Ωup=∫Ω∗0(u∗)p. | (1.13) |
Two results are particularly important in the study of the Schwarz rearrangement:
● Hardy-Littlewood inequality: for any two non-negative functions f,g∈L2(Ω),
∫Ω0fg⩽∫Ω∗0f∗g∗. | (1.14) |
● PólyaSzegö inequality: for any p∈[1;+∞), for any u∈W1,p0(Ω),u⩾0,
u∗∈W1,p0(Ω0)and∫Ω∗0|∇u∗|p⩽∫Ω0|∇u|p. | (1.15) |
Finally, we will rely on an ordering of the set of functions.
Definition 4. Let RΩ0>0 be the radius of the ball Ω∗0. For any two non-negative functions f,g∈L1(Ω0) we write
f≺g |
if
∀r∈[0;RΩ0],∫B(0;r)f∗⩽∫B(0;r)g∗. | (1.16) |
This ordering [17] provides the natural framework for comparison theorems in elliptic and parabolic equations [3,4,28,29,36,37,38]. The following property is proved in [4,Proposition 2]: for any non-decreasing convex function F such that F(0)=0, for any two non-negative functions f,g∈L1(Ω),
f≺g→F(f)≺F(g). | (1.17) |
We now pass to the definition of rearrangement sets:
Definition 5 (Rearrangement sets). For any non-negative function f∈L1(Ω0) we define
CΩ0(f):={φ∈L1(Ω0),φ∗=f∗} | (1.18) |
and
KΩ0(f):={φ∈L1(Ω0),φ⩾0a.e.,φ≺f,∫Ω0φ=∫Ω0f}. | (1.19) |
The following result can be found in [8,27,34]: for any non-negative f∈L1(Ω), KΩ(f) is a weakly compact, convex set; its extreme points are the elements of CΩ(f).
The Alvino-Trombetti rearrangement: definition and property The Alvino-Trombetti rearrangement is very useful when handling two-phase problems, and was introduced in [6,7] to establish some comparison principles for some elliptic equations with a diffusion matrix. The goal is the following: let u∈W1,20(Ω)∩L∞(Ω) be a non-negative function and let a∈A(Ω). We want to prove that there exists ˜a that is radially symmetric, such that ˜aId is uniformly elliptic and such that
∫Ωa|∇u|2⩾∫Ω∗˜a|∇u∗|2. | (1.20) |
One defines ˜a as the unique radially symmetric function such that
Fora.e.t∈(0;‖ | (1.21) |
It can be checked [15] that
\tilde a^{-1} \in \mathscr K_{{{\Omega}}^*}\left((a^*)^{-1}\right). |
Remark 6. In particular, if all the level-sets of u have Lebesgue measure zero and the gradient of u does not vanish on these level-sets, this definition rewrites as
\begin{equation} {\rm{ For\;a.e. }}\; t\in (0;\Vert u\Vert_{L^\infty})\, , \int_{\{u^* = t\}}\frac1{\tilde a |{{\nabla}} u^*|} = \int_{\{u = t\}}\frac1{a|{{\nabla}} u|}. \end{equation} | (1.22) |
This fact follows from the co-area formula, which states in particular that
\int_{\{u \geqslant c\}}\frac1a = \int_c^{\infty}\int_{\{u = t\}}\frac1{a|{{\nabla}} u|}. |
[7,Lemma 1.2] or [15,Proposition 4.9] assert that: for any u \geqslant 0\, , u\in W^{1, 2}_0({{\Omega}}) , for any a\in \mathcal A({{\Omega}}) , \tilde a being defined by (1.21), there holds
\begin{equation} \int_{{{\Omega}}}a|{{\nabla}} u|^2 \geqslant \int_{{{\Omega}}^*}\tilde a|{{\nabla}} u^*|^2. \end{equation} | (1.23) |
Talenti inequalities for the relaxed problem Let us start by recalling an application of the Alvino-Trombetti rearrangement to Talenti-like inequalities. Talenti inequalities originated in the seminal [36] and have, since then, been widely studied [2,3,4,8,9,28,29,37,38]. Roughly speaking, they amount to comparing, using the relation \prec , the solution u of an elliptic problem with the solution u' of a "symmetrised" elliptic equation. This first result [6] is the stepping stone to our main theorem and holds for the relaxed version of the problem:
Theorem A ([6], Comparison results, optimisation w.r.t. a and f ). Let {{\Omega}} be a centred ball. For any a\in \mathcal A({{\Omega}}) and any f\in \mathcal F({{\Omega}}) , \tilde a being defined by (1.21), there holds
\begin{equation} u_{{{\operatorname{ell}}}, a, f}^* \leqslant u_{{{\operatorname{ell}}}, \tilde a, f^*}. \end{equation} | (1.24) |
As a consequence, for any increasing function j ,
\int_\Omega j(u_{{{\operatorname{ell}}}, a, f}) \leqslant \int_\Omega j(u_{{{\operatorname{ell}}}, \tilde a , f^*}). |
Thus, it seems quite interesting to investigate whether the optimisation problem ({\bf{P}}_{{\rm{ell}}, j}) has a radial solution. This would seem natural given the equation above. However, the Alvino-Trombetti rearrangement only provides us with a rearranged coefficient \tilde a such that the inverse (\tilde a)^{-1}\in \mathscr K_{{{\Omega}}^*}(a^{-1}) . This last set is however different from \mathcal A({{\Omega}}) . The same problem arises when considering bang-bang functions a . We nonetheless obtain existence properties for the unrelaxed problem.
Theorem I (Existence and bang-bang property in radial geometry for convex integrand, optimisation w.r.t. a and f ). Assume j is a convex C^2 function. Let R > 0 . Let {{\Omega}} = \mathbb B(0;R) . The optimisation problem
\sup\limits_{a\in \mathfrak M({{\Omega}}), f\in \mathcal F({{\Omega}})}\int_{{{\Omega}}} j({{u_{{{\operatorname{ell}}}, a, f}}}) |
has a solution (\overline a, \overline f)\in \mathfrak M({{\Omega}})\times \mathcal F({{\Omega}}) .
The proof of this theorem is inspired by the proof of existence of optimal profiles for eigenvalue problems in [15].
Finally, the last result for elliptic problems deals with a bang-bang property when optimising only with respect to a : is it true that, if we just assume that j is increasing, if a solution \overline{a} of ({\bf{P}}_{{\rm{ell}}, j, f}) exists, then it is bang-bang? We can only partially answer this question, in the next theorem. It is the main result of our paper.
Theorem II (Weak bang-bang property for increasing cost functions, optimisation w.r.t radially symmetric a ). Assume j is an increasing function such that j' > 0 on \mathit{{\rm{I\kern-0.21emR}}\; }_+^* . Let R > 0 . Let {{{\Omega}}} = \mathbb B(0;R) and f\in \mathcal F({{{\Omega}}}) . Then, if the optimisation problem
\sup\limits_{a\in \mathcal A({{{\Omega}}}), a\;\mathit{radially\;symmetric}}\int_{{{{\Omega}}}} j({{u_{{{\operatorname{ell}}}, a, f}}})\qquad\qquad\qquad{\bf{P}}_{{{\operatorname{ell}}}, j, a} |
has a solution \overline a , there holds
\overline a\in \mathfrak M({{{\Omega}}})\mathit{or, in\;other\;words, }\overline a = \alpha+(\beta-\alpha)1_\omega\;\mathit{for\;some\;measurable}\omega\subset {{\Omega}} . |
The proof of this theorem is based on the development of an oscillatory method recently introduced in [25].
In this second part, we state our main result devoted to the parabolic optimisation problem ({\bf{P}}_{{\rm{parab}}, j1, j2}). The proof of the parabolic isoperimetric inequality is done by adapting the proofs of Theorem A and of [29,Theorem 2.1]. For the sake of clarity, for a function u of two variables u = u(t, x) , the notation u^*(t, \cdot) stands for the Schwarz rearrangement of u(t, \cdot) with respect to the space variable x .
Theorem III (Comparison results, optimisation w.r.t. a and f ). Let {{\Omega}} = \mathbb B(0;R) .Let a\in \mathcal M({{\Omega}}; T) and f\in \mathcal F({{\Omega}}; T) . Then there exists a radially symmetric function \tilde a defined on (0;T)\times {{\Omega}} such that \alpha \leqslant \tilde a \leqslant \beta almost everywhere and such that, for almost every t\in (0;T) and every r\in ({{\Omega}}; R) there holds
\int_{{{\mathbb B}}(0;r)}u^*_{{{\operatorname{parab}}}, a, f} \leqslant \int_{{{\mathbb B}}(0;r)}u_{{{\operatorname{parab}}}, \tilde a, f^*}. |
In particular, if j_1 and j_2 are convex increasing functions there holds
\mathcal J_{{{\operatorname{parab}}}}(a, f) \leqslant \mathcal J_{{{\operatorname{parab}}}}(\tilde a, f^*). |
Let us now offer some comments about this result, and about the method of proof.
Remark 7 (Comments on Theorem III). 1) The first thing that has to be noted is that, exactly as in the elliptic case, although the new weight \overline a satisfies the correct upper and lower bounds \alpha \leqslant \overline \beta , there is a priori no guarantee that \overline a\in \mathcal M({{\Omega}}) . Some other arguments would then be needed in order to conclude as to the integral constraint. It is not clear at this stage how one may go about this question.
2) The second remark has to do with the method of proof that is employed. The two main available approaches in the context of parabolic equations are, on the one hand, dealing with the parabolic problem directly, as is done in [29] and as we do here, and on the other hand by time-discretisation of the evolution problem, as in [3]. We believe the second of these approaches may prove more delicate. To see why, let us recall the main steps of the proof of [3]: the authors, which in particular try to prove a comparison result for rearrangement of the source f in the parabolic equation
\frac{\partial u}{\partial t}-\Delta u = f, |
approximate this equation by the discretisation with time step N\in {\rm{I\kern-0.21emN}}\; ^*
N({u_{k+1, N}-u_{k, N}})-\Delta u_{k+1, N} = f_{k, N} |
where f_{k, N} = N\int_{k/N}^{(k+1)/N}f. On each of these discretised problem, they use an elliptic Talenti inequality, yielding a comparison with the solution of the same system with f_{k, N}^* as a right-hand side. This Schwarz symmetrisation operation is independent of the time-step (in the sense that the definition of f^* does not depend on the time step N ). In our case however, since the Alvino-Trombetti rearrangement depends on the function u_{{{\operatorname{parab}}}, a, f} evaluated at the time t , this would translate, at the discretised level, as a rearrangement that would depend on both indexes k and N . This may lead to potential difficulties in passing to the limit.
Let us underline that this type of parabolic comparison results can be very useful when dealing with parabolic eigenvalue optimisation problems, as is done for instance in [32,Theorem 3.9].
In this paper, we offer contributions that may be viewed from several point of views, each of which stemming from very rich domains in mathematical analysis.
Two-phase spectral optimisation problems Two-phase optimisation problems have a rich history, and are deeply linked to homogenisation phenomenas. We refer, for instance, to [1,31] for a presentation of this rich theory, and we underline that one of the striking features of these problems is that there is often a lack of existence results. These results are typically obtained by proving that should an optimiser exist, then an overdetermined problem that can only solved in radial geometries should have a solution. This is done by using Serrin type theorems [35], and this phenomenon occurs in dimension d \geqslant 2 . A typical and famous example of such problems is the optimisation of the first Dirichlet eigenvalue of the operator -\nabla \cdot(a\nabla) under the constraint that a\in \mathcal A({{\Omega}}) . To the best of our knowledge, the proof of non-existence of an optimal a^*\in \mathcal A({{\Omega}}) when {{\Omega}} is not a ball was only recently completed in a series of papers by Casado-Diaz [10,11,12]. However, these negative results in the case of non-radial geometries do not allow to conclude as for the existence and/or characterisation of optimisers in radially symmetric domains. In this case, the same spectral optimisation problem being under consideration, the first proof of existence can be found in [15], using the Alvino-Trombetti rearrangement. We borrow from their ideas in the proof of Theorem I (and we highlight the fact that we do not consider here energetic problems). To underline the complexity of this spectral optimisation problem, let us also mention [22], in which it is shown that, in the ball, the qualitative features of the optimiser a^* strongly depend on the volume constraint. We also refer to [13,26] for the study of the spectral optimisation of operators with respect to a weight a\in \mathcal A({{\Omega}}) that appears both in the principal symbol -\nabla\cdot(a{{\nabla}}) and as a potential.
Elliptic and parabolic Talenti inequalities Talenti inequalities, which originate in the seminal [36] have been the subject of an intense research activity. For parabolic equations, the study of such inequalities started, as far as we are aware, in the works of Bandle [9], Vazquez [38] and were later deeply analysed by Alvino, Trombetti and Lions [2,3] on the one hand, and by Mossino and Rakotoson on the other [29]. We would like to mention that we have recently obtained a quantitative parabolic isoperimetric inequality for the source term in [24]. Alvino, Nitsch and Trombetti have recently proved an elliptic Talenti inequality under Robin boundary conditions, using a very fine analysis of the Robin problem [5]. This Robin Talenti inequality was then used in, for instance, [21,33].
Proof of Theorem I. For the first part of the theorem, we consider the case {{{\Omega}}} = \mathbb B(0;R) where R > 0 is a fixed constant. We work with functions a\in \mathfrak M({{\Omega}}) . In other words, there exists \omega\subset {{\Omega}} measurable such that
a = \alpha+(\beta-\alpha)1_\omega, |
and we aim at solving
\sup\limits_{a\in \mathfrak M({{{\Omega}}}), f \in \mathcal F({{{\Omega}}})} \mathcal J(f, a) = \int_{{{{\Omega}}}}j({{u_{{{\operatorname{ell}}}, a, f}}}), |
under the assumption that j is convex on {\rm{I\kern-0.21emR}}\; _+ .
Let us first note that for any a\in \mathfrak M({{\Omega}}) we have a^* = \alpha+(\beta-\alpha)1_{{{\mathbb B}}^*} where {{\mathbb B}}^* = {{\mathbb B}}(0;r^*) satisfies \beta {{\operatorname{Vol}}}({{\mathbb B}}^*)+{{\operatorname{Vol}}}({{{\Omega}}}\backslash {{\mathbb B}}^*) = V_0 . As a consequence, for any a_1, a_2\in \mathfrak M({{\Omega}}) ,
\mathscr K((a_1)^{-1}) = \mathscr K((a_2)^{-1}) |
where \mathscr K(\cdot) is the rearrangement class defined in definition 5. For the sake of notational convenience, we define
\mathfrak K: = \mathscr K(a^{-1})\;{\rm{ where\;a\;is\;any\;element\;of}}\;\mathfrak M({{{\Omega}}\;}) . |
By Theorem A, for any a\in \mathcal A ({{{\Omega}}}) there exists a radially symmetric \tilde a such that
\tilde a^{-1}\in \mathfrak K\, , J(\tilde a, \overline f) \geqslant J(a, f) |
where \overline f is simply the Schwarz rearrangement of f . By convexity of the functional with respect to f , \overline f is a bang-bang function. We henceforth consider it fixed and focus on optimisation with respect to a .
The problem with this reformulation is that there is a priori no guarantee that \tilde a\in \mathfrak M({{{\Omega}}}) , and it is in general false. To overcome this difficulty, we now focus on a slightly simplified version of our problem:
\begin{equation} \sup\limits_{ a\;{\rm{ radially\;symmetric\;s.t. }}\;\tilde a\in \mathfrak K}\int_{{{{\Omega}}}}j(u_{ a})\;{\rm{ where}}\;u_{ a}\;{\rm{solves}}\;\begin{cases}-\nabla \cdot ( a \nabla u_{ a}) = \overline f&{\rm{ in }}\;{{{\Omega}}}\, , \\ u_{ a} = 0&{\rm{ on }}\;\partial {{{\Omega}}}.\end{cases} \end{equation} | (2.1) |
Another refomulation, a priori encompassing a larger class, is
\begin{equation} \sup\limits_{\tilde \mu\;{\rm{ radially\;symmetric\;s.t. }}\;\mu\in \mathfrak K}\left(\mathcal H(\mu): = \int_{{{{\Omega}}}}j(v_\mu)\right)\;{\rm{ where}}\;v_\mu\;{\rm{solves}}\;\begin{cases}-\nabla \cdot \left(\frac1\mu \nabla v_\mu\right) = \overline f&{\rm{ in }}\;{{{\Omega}}}\, , \\ v_\mu = 0&{\rm{ on }}\;\partial {{{\Omega}}}.\end{cases} \end{equation} | (2.2) |
We now proceed in several steps, following the ideas of [15]:
1) Existence of solutions to (2.2): we first prove, in lemma 1, that there exists a solution \overline \mu to (2.2). This is done via the direct method in the calculus of variations.
2) The bang-bang property for \overline \mu : we then prove, in Lemma 9, that any solution \overline \mu of (2.2) is a bang-bang function. In other words, there exists a measurable subset \omega\subset {{{\Omega}}} such that
\overline\mu = \frac1{\alpha+(\beta-\alpha) 1_\omega}. |
This is done via a convexity argument. As a consequence, \frac1{\overline \mu}\in \mathfrak M({{{\Omega}}}) , hence concluding the proof.
Existence of solutions to (2.2) The main result of this paragraph is the following lemma:
Lemma 8. There exists a solution {{\overline \mu}} of the variational problem (2.2).
Proof of Lemma 8. We first note the following fact: if a sequence of radially symmetric functions \{\mu_k\}\in \mathfrak K^ {\rm{I\kern-0.21emN}}\; weakly converges in L^\infty-* to \mu_\infty (which is an element of \mathfrak K by the closedness of \mathfrak K for the weak convergence, see [27]) then, up to a subsequence,
\mathcal H(\mu_k)\underset{k\to\infty}\rightarrow \mathcal H(\mu_\infty). |
Obtaining this result boils down to proving that, for any sequence of radially symmetric functions \{\mu_k\}\in \mathfrak K^ {\rm{I\kern-0.21emN}}\; weakly converges in L^\infty-* to \mu_\infty there holds
\left(\mu_k\underset{k\to \infty}\rightharpoonup \mu_\infty\right)\Rightarrow \left(v_{\mu_k}\underset{k\to \infty}\rightarrow v_{\mu_\infty}\;{\rm{ a.e.\;up\;to\;a\;subsequence}}\;\right). |
Indeed, it then simply suffices to use the dominated convergence theorem to obtain the required result. Let us then prove that for any sequence of radially symmetric functions \{\mu_k\}\in \mathfrak K^ {\rm{I\kern-0.21emN}}\; weakly converging in L^\infty-* to \mu_\infty ,
\begin{equation} v_{\mu_k}\underset{k\to \infty}\rightarrow v_{\mu_\infty}{\rm{ in L^2({{\Omega}}\;}) }. \end{equation} | (2.3) |
However, since we are working with radially symmetric functions, this follows from explicit integration in radial coordinates of
\begin{cases}-\nabla \cdot \left(\frac1\mu \nabla v_{\mu_k}\right) = \overline f&{\rm{ in }}\;{{{\Omega}}}\, , \\ v_{\mu_k} = 0&{\rm{ on }}\;\partial {{{\Omega}}}\end{cases} |
which gives, for any k\in {\rm{I\kern-0.21emN}}\; (and with a slight abuse of notation),
r^{d-1}v_{\mu_k}'(r) = -\mu_k(r)\int_0^r s^{d-1}f(s)ds. |
Thus, from the Rellich-Kondrachov embedding
v_{\mu_k}\underset{k\to \infty}\rightarrow v_{\mu_\infty}\begin{cases}{\rm{ weakly\;in }}\;W^{1, 2}_0({{{\Omega}}})\, , \\ {\rm{ strongly\;in\;}}\;L^2({{{\Omega}}}).\end{cases} |
It suffices to extract a subsequence that is converging almost everywhere.
We turn back to the proof of the lemma: let \{\mu_k\}_{k\in {\rm{I\kern-0.21emN}}\; } be a maximising sequence for (2.2). Since the set \mathfrak K is weakly compact, and since for any k\in {\rm{I\kern-0.21emN}}\; \mu_k is radially symmetric, there exists a radially symmetric \mu_\infty\in \mathfrak K such that, up to a subsequence,
\mu_k\underset{k\to \infty}\rightarrow \mu_\infty \;{\rm{ weakly\;in}} \; L^\infty-* . |
Hence, up to a subsequence,
\mathcal H(\mu_k)\underset{k\to \infty}\rightarrow \mathcal H(\mu_\infty) |
so that \mu_\infty is a solution of (2.2).
The bang-bang property for \overline \mu We now present the key point of the proof of Theorem I, the bang-bang property.
Lemma 9. Any solution \overline \mu of (2.2) is of bang-bang type: there exists \omega\subset {{{\Omega}}} such that
\mu = \frac1{\alpha+(\beta-\alpha)1_\omega}. |
Proof of Lemma 9. We argue by contradiction and assume that there exists a solution \overline \mu of (2.2) that is not of bang-bang type. We will reach a conclusion using a second order information on the functional \mathcal H , namely, by using the first and second order Gâteaux-derivative of the functional \mathcal H . Let us first observe that it is standard [18] to see that the map \mathfrak K\ni \mu \mapsto v_\mu is Gâteaux-differentiable. Furthermore, for a given \mu \in \mathfrak K and an admissible perturbation h at \mu (i.e., such that \mu +th\in \mathfrak K for t > 0 small enough) the first order Gâteau-derivative of v_\mu in the direction h is the unique solution \dot v_\mu(h) of
\begin{equation} \begin{cases}-{{\nabla}} \cdot\left(\frac1\mu {{\nabla}} \dot v_\mu(h)\right)+{{\nabla}} \cdot\left(\frac{h}{\mu^2}{{\nabla}} v_\mu\right) = 0&{\rm{ in }}\;{{{\Omega}}}\, , \\ \dot v_\mu(h) = 0&{\rm{ on }}\;\partial {{{\Omega}}}, \end{cases} \end{equation} | (2.4) |
and the first Gâteaux-derivative of \mathcal H at \mu in the direction h is given by
\begin{equation} \dot{\mathcal H}(\mu)[h] = \int_{{{{\Omega}}}} j'(v_\mu)\dot v_{\mu}(h). \end{equation} | (2.5) |
This leads to introducing the adjoint state p_\mu as the unique solution of
\begin{equation} \begin{cases}-{{\nabla}} \cdot\left(\frac1\mu {{\nabla}} p_\mu\right) = j'(v_\mu)&{\rm{ in }}\;{{{\Omega}}}\, , \\ p_\mu = 0&{\rm{ on }}\;\partial {{{\Omega}}}, \end{cases} \end{equation} | (2.6) |
Remark 10. It should be noted that by explicit integration of the equation on v_\mu in radial coordinate, v_\mu\in L^\infty({{{\Omega}}}) ; as a consequence, j'(v_\mu) is an L^\infty function, so that p_\mu is well-defined.
Multiplying, on the one hand (2.4) by p_\mu , on the other hand (2.6) by v_\mu , and integrating by parts gives
\begin{equation} \dot{\mathcal H}(\mu)[h] = \int_{{{{\Omega}}}}\frac1\mu \langle {{\nabla}} p_\mu, {{\nabla}} \dot v_\mu(h)\rangle = -\int_{{{{\Omega}}}}\frac{h}{\mu^2}\langle {{\nabla}} v_\mu\, , {{\nabla}} p_\mu\rangle = -\int_{{{{\Omega}}}}h\langle \frac{{{\nabla}} v_\mu}{\mu}, \frac{{{\nabla}} p_\mu}{\mu}\rangle. \end{equation} | (2.7) |
We now compute the second order derivative of the criterion in a similar manner: the second order Gâteaux derivative of v_\mu at \mu in the direction h is zero. In other words, denoting by \ddot v_\mu this second order derivative, we have
\ddot v_\mu = 0. |
Indeed, this follows from the explicit computation of v_\mu as
r^{d-1}v_{\mu_k}'(r) = -\mu_k(r)\int_0^r s^{d-1}f(s)ds. |
Thus, it appears that \mu \mapsto v_\mu is linear. As a consequence we have that the second order Gâteaux derivative of \mathcal H at \mu in the direction h , henceforth abbreviated as \ddot{\mathcal H}(\mu) , is given by
\begin{equation} \ddot{\mathcal H}(\mu) = \int_{{{{\Omega}}}} j''(v_\mu)\left(\dot v_\mu\right)^2+\int_{{{{\Omega}}}} j'(v_\mu) \ddot v_\mu = \int_{{{{\Omega}}}}j''(v_\mu)\left(\dot v_\mu\right)^2. \end{equation} | (2.8) |
Hence, if j is convex, so is \mathcal H . Thus any solution \overline \mu of (2.2) is an extreme point of \mathfrak K . In other words
\overline \mu\in C(a^{-1})\; {\rm{ where\;a\;is\;any\;element\;of}}\;\mathfrak M({{{\Omega}}\;) }. |
It follows that \frac1{\overline\mu}\in \mathfrak M({{{\Omega}}}).
Conclusion of the proof As noted at the beginning of the proof, for any a\in \mathfrak M({{{\Omega}}}) there exists \tilde a such that \tilde a^{-1}\in \mathfrak K and such that \mathcal J(\tilde a, \overline f) \geqslant \mathcal J(a, f) . Since
\mathcal H\left(\frac1{\tilde a}\right) = \mathcal J(\tilde a, \overline f) |
it follows that
\mathcal J(\tilde a, \overline f) \leqslant \mathcal H(\overline \mu){\rm{ where}}\;\overline \mu\;{\rm{is \;the\; solution \;of \;(2.2)\; given\; by\;lemma\;8}}. |
From proposition 2 \overline\mu is a bang-bang function. As a consequence, \frac1{\overline \mu}: = \overline a is an element of \mathfrak M({{{\Omega}}}) . Thus
\mathcal J(a, f) \leqslant \mathcal J(\tilde a, \overline f) \leqslant \mathcal H(\overline \mu) = \mathcal J(\overline a, \overline f). |
Thus \overline a is a solution of the initial optimisation problem.
The proof of the theorem is now complete.
Proof of Theorem II. Throughout this proof we assume that we are given a radially symmetric solution \overline a of the optimisation problem
\sup\limits_{a\in \mathcal A({{{\Omega}}}), a {\rm{ radially symmetric}}\;}\int_{{{{\Omega}}}} j({{u_{{{\operatorname{ell}}}, a, f}}}) |
and we want to prove that \overline a \in \mathfrak M({{{\Omega}}}) . To reach the desired conclusion we argue by contradiction and we assume that \overline a\notin \mathfrak M({{{\Omega}}}) . We emphasise once again that this proof does not require rearranging the source term f . Since f is assumed to be fixed, we write \mathcal J(a) for \mathcal J(a, f) = \int_{{{{\Omega}}}}j({{u_{{{\operatorname{ell}}}, a, f}}}) and u_a for {{u_{{{\operatorname{ell}}}, a, f}}} .
Let us single out the following result, that follows from direct integration in radial coordinates of (1.10):
Lemma 11. For any radially symmetric a\in \mathcal A({{{\Omega}}}) and f\in \mathcal F({{{\Omega}}}) , u_a\in W^{1, \infty}({{{\Omega}}}) , u_a is radial and we furthermore have, with a slight abuse of notation, for a.e. r\in (0;R) ,
\begin{equation} u_a'(r) = -\frac1{a(r)r^{d-1}}\int_0^r \xi^{d-1} f(\xi)d\xi. \end{equation} | (3.1) |
In particular, u_a' is a non-positive function and, for any {{\varepsilon}} > 0 , \sup_{[{{\varepsilon}}; R]}u_a' < 0 . It is strictly decreasing if f > 0 in a neighbourhood of 0.
We now compute the Gateaux derivatives of both the maps a\mapsto {{u_{{{\operatorname{ell}}}, a, f}}} and of a\mapsto \mathcal J(a) (we note that the fact that both maps are Gateaux differentiable follow from standard arguments). We note that, to compute them, it is not necessary to assume that the coefficients a and f are radiallly symmetric.
The first-order Gateaux derivative of {{u_{{{\operatorname{ell}}}, a, f}}} at a in an admissible direction h (i.e., such that a+th\in \mathcal A({{{\Omega}}}) for t > 0 small enough), denoted by \dot u_a , is the unique solution to
\begin{equation} \begin{cases} -\nabla \cdot(a{{\nabla}} \dot u_a) = {{\nabla}} \cdot(h{{\nabla}} u_a)&{\rm{ in }}\;{{{\Omega}}}\, , \\ \dot u_a = 0&{\rm{ on }}\;\partial {{{\Omega}}}. \end{cases} \end{equation} | (3.2) |
The Gateaux derivative of \mathcal J at a in the direction h is given by
\begin{equation} \dot{\mathcal J}(a)[h] = \int_{{{{\Omega}}}}j'(u_a) \dot u_a. \end{equation} | (3.3) |
This leads to introducing the adjoint state p_a as the unique solution to
\begin{equation} \begin{cases} -{{\nabla}} \cdot(a{{\nabla}} p_a) = j'(u_a)&{\rm{ in }}\;{{{\Omega}}}\, , \\ p_a = 0&{\rm{ on }}\;\partial {{{\Omega}}}. \end{cases} \end{equation} | (3.4) |
Multiplying (3.4) by \dot u_a and (3.2) by p_a and integrating by parts gives
\begin{equation} \dot{\mathcal J}(a)[h] = \int_{{{{\Omega}}}}j'(u_a)\dot u_a = \int_{{{{\Omega}}}}a \langle {{\nabla}} p_a\, , {{\nabla}} \dot u_a\rangle = -\int_{{{{\Omega}}}} h\langle {{\nabla}} u_a, {{\nabla}} p_a\rangle. \end{equation} | (3.5) |
In the same way, the second order Gateaux derivative of u_a at a in the direction h , denoted by \ddot u_a , is the unique solution to
\begin{equation} \begin{cases} -\nabla \cdot(a{{\nabla}} \ddot u_a) = 2{{\nabla}} \cdot(h{{\nabla}}\dot u_a)&{\rm{ in }}\;{{{\Omega}}}\, , \\ \ddot u_a = 0&{\rm{ on }}\;\partial {{{\Omega}}}, \end{cases} \end{equation} | (3.6) |
and the second order Gateaux derivative of \mathcal J at a in the direction h is given by
\begin{equation} \ddot{\mathcal J}(a)[h, h] = \int_{{{{\Omega}}}} j''(u_a)(\dot u_a)^2+\int_{{{{\Omega}}}} j'(u_a)\ddot u_a. \end{equation} | (3.7) |
However, multiplying (3.6) by p_a , integrating by parts and using the weak formulation of (3.4) yields
\begin{equation} \int_{{{{\Omega}}}} j'(u_a)\ddot u_a = \int_{{{{\Omega}}}}a \langle {{\nabla}} p_a\, , {{\nabla}} \ddot u_a\rangle = -2\int_{{{{\Omega}}}} h\langle {{\nabla}} \dot u_a\, , {{\nabla}} p_a\rangle. \end{equation} | (3.8) |
Plugging (3.8) in (3.7) gives
\begin{equation} \ddot{\mathcal J}(a)[h, h] = \int_{{{{\Omega}}}} j''(u_a)(\dot u_a)^2-2\int_{{{{\Omega}}}} h\langle {{\nabla}} \dot u_a\, , {{\nabla}} p_a\rangle. \end{equation} | (3.9) |
We now use the radial symmetry assumption: since h, a and f are radially symmetric, and since u_a'(0) = \dot u_a'(0) = 0 , (3.2) implies, in radial coordinates, as
\begin{equation} -a \dot u_a' = h u_a'. \end{equation} | (3.10) |
Furthermore, we have the following lemma:
Lemma 12. If j' > 0 on \mathit{{\rm{I\kern-0.21emR}}\; }_+^* then p_a is a radially symmetric decreasing function:
p_a' < 0\;in\;(0;R). |
Proof of Lemma 12. The fact that p_a is decreasing simply follows from, first, the strong maximum principle which implies that
u_a > 0\;{\rm{ in }}\;[0;R) |
and, second, from explicit integration of the equation on p_a in radial coordinates, which gives
p_a'(r) = -\frac1{a(r)r^{d-1}}\int_0^r j'(u_a)r^{d-1}dr < 0 \;{\rm{ for}}\;r > 0 .\; |
The radiality of p_a implies
\langle {{\nabla}} \dot u_a\, , {{\nabla}} p_a\rangle = \dot u_a' p_a'. |
As a consequence of (3.10), we have, for a constant M_d > 0
\begin{align*} -\int_{{{{\Omega}}}} h \langle {{\nabla}} \dot u_a\, , {{\nabla}} \dot p_a\rangle& = -M_d \int_0^R r^{d-1} h(r) \dot u_a'(r)p_a'(r)dr \\& = M_d\int_0^R a(r) \left(\dot u_a'\right)^2 \frac{p_a'}{u_a'}(r)dr \\& = M_d\int_0^R a(r) \left(\dot u_a'\right)^2 \frac{ap_a'}{au_a'}(r)dr. \end{align*} |
Let us first define
{{\varphi}}:r\in (0;R]\mapsto \frac{ap_a'}{au_a'}(r). |
We observe that {{\varphi}} > 0 in (0;R] and that, as r\to 0 ,
\begin{equation} \underset{r\to 0}{\lim \inf} {{\varphi}}(r) \geqslant \underset{r\to 0}{\lim \inf}\frac{j'(u_a(0))}{f(r)} > 0. \end{equation} | (3.11) |
If f > 0 in a neighbourhood of 0 then we can extend {{\varphi}} by \frac{j'(u_a(0))}{f(0)} > 0 in 0. If on the other hand f = 0 in a neighbourhood of 0, then {{\varphi}} \to +\infty when r\to 0 . In any case, there exists a constant A > 0 such that
{{\varphi}} \geqslant\frac{ A_0}2 > 0\;{\rm{ in }}\;[0;R]. |
We define the function
\Phi:{{{\Omega}}}\ni x\mapsto {{\varphi}}\left(|x|\right), |
and we thus have
\begin{align*} -\int_{{{{\Omega}}}} h \langle {{\nabla}} \dot u_a\, , {{\nabla}} \dot p_a\rangle& = -2M_d \int_0^R r^{d-1} h(r) \dot u_a'(r)p_a'(r)dr \\& = 2M_d\int_0^R a(r) \left(\dot u_a'\right)^2 \frac{p_a'}{u_a'}(r)dr \\& = 2\int_{{{{\Omega}}}}\Phi a|{{\nabla}} \dot u_a|^2 \\& \geqslant A_0\int_{{{{\Omega}}}} a|{{\nabla}} \dot u_a|^2. \end{align*} |
Finally, as j\in C^2({\rm{I\kern-0.21emR}}\; _+) and u_a\in L^\infty({{{\Omega}}}) there exists a constant B > 0 such that
\begin{equation} j''(u_a) \geqslant -B\;{\rm{ in }}\;{{{\Omega}}}. \end{equation} | (3.12) |
We end up with the following estimate on \ddot{\mathcal J}(a)[h, h] :
\begin{align} \ddot{\mathcal J}(a)[h, h]& = \int_{{{{\Omega}}}} j''(u_a)(\dot u_a)^2-2\int_{{{{\Omega}}}} h\langle {{\nabla}} \dot u_a\, , {{\nabla}} p_a\rangle \end{align} | (3.13) |
\begin{align} & \geqslant A_0\int_{{{{\Omega}}}}a |{{\nabla}} \dot u_a|^2-B\int_{{{{\Omega}}}}\dot u_a^2. \end{align} | (3.14) |
Remark 13. It should be noted that, at this level, we recover the convexity of the functional \mathcal J if we assume that j'' \geqslant 0 . Indeed, in that case we can take B = 0 .
Let us now turn back to the core of the proof: we have a maximiser \overline a \notin \mathfrak M({{{\Omega}}}) . Let
\tilde\omega: = \{\alpha < a < \beta\}. |
Since \overline a\notin \mathfrak M({{{\Omega}}}) ,
{{\operatorname{Vol}}}(\tilde \omega) > 0. |
Furthermore, for any h\in L^\infty(\tilde \omega) , extended by 0 outside of \tilde \omega and such that \int_{\tilde \omega}h = 0 , we have (since both h and -h are admissible perturbations at \overline a )
\begin{equation} \dot{\mathcal J}(a)[h] = 0. \end{equation} | (3.15) |
To reach a contradiction, it suffices to build h\in L^\infty(\tilde \omega)\backslash\{0\} such that
\begin{equation} \boldsymbol{\int_{{{{\Omega}}}}h 1_{{{\tilde{\omega}}}} = 0}\;{\rm{ and }}\; \boldsymbol{\ddot{\mathcal J}(a)[h1_{{\tilde{\omega}}}, h1_{{\tilde{\omega}}}] > 0.} \end{equation} | (3.16) |
Actually, by approximation it suffices to build h\in L^2({{\tilde{\omega}}}) such that (3.16) holds. For the sake of notational simplicity, for any h\in L^2({{\tilde{\omega}}}) , we identify h with h1_{{\tilde{\omega}}}\in L^2({{{\Omega}}}) . To obtain the existence of such a perturbation we single out the estimate
\begin{equation} \ddot{\mathcal J}(a)[h, h] \geqslant A_0\int_{{{{\Omega}}}}a |{{\nabla}} \dot u_a|^2-B\int_{{{{\Omega}}}}\dot u_a^2. \end{equation} | (3.17) |
We introduce the sequence \{\sigma_k\, , \psi_k\}_{k\in {\rm{I\kern-0.21emN}}\; ^*} of eigenvalues of the operator -\nabla\cdot(a{{\nabla}}) . We pick a non-decreasing of the eigenvalue sequence:
0 < \sigma_0 \leqslant \sigma_1 \leqslant \dots \leqslant \sigma_k\underset{k\to \infty}\rightarrow +\infty. |
The eigenequations are given by
\begin{equation} \forall k\in {\rm{I\kern-0.21emN}}\;^*\, , \begin{cases}-\nabla \cdot(a{{\nabla}} \psi_k) = \sigma_k\psi_k&{\rm{ in }}\;{{{\Omega}}}\, , \\ \psi_k = 0&{\rm{ on }}\;\partial {{{\Omega}}}\, , \\ \int_{{{{\Omega}}}} \psi_k^2 = 1.\end{cases} \end{equation} | (3.18) |
For any admissible perturbation h at a , we decompose \dot u_a in this basis as
\begin{equation} \dot u_a = \sum\limits_{k = 1}^\infty\alpha_k(h) \psi_k, \end{equation} | (3.19) |
where the coefficients \left\{\alpha_k(h)\right\}_{k\in {\rm{I\kern-0.21emN}}\; ^*} are determined by equation (3.2). If we assume that, for an integer K large enough, we have
\begin{equation} \forall k \leqslant K\, , \alpha_k(h) = 0\, , \sum\limits_{k = K}^\infty\alpha_k(h)^2 > 0 \end{equation} | (3.20) |
then we obtain, by expanding the right hand-side of (3.17),
\begin{equation} \boldsymbol{\ddot{\mathcal J}(a)[h1_{{\tilde{\omega}}}, h1_{{\tilde{\omega}}}] \geqslant A\sum\limits_{k = K}^\infty\sigma_k {\alpha_k(h)^2}-B\sum\limits_{k = K}^\infty{\alpha_k(h)^2} \geqslant \left(A\sigma_K-{B}{}\right)\sum\limits_{k = K}^\infty \frac{\alpha_k(h)^2}{\sigma_k} > 0.} \end{equation} | (3.21) |
As a consequence it remains to construct a perturbation h such that
\begin{equation} \nabla\cdot (h{{\nabla}} u_a) = \sum\limits_{k \geqslant K}\eta_k \psi_k. \end{equation} | (3.22) |
We need however to be careful, since \nabla \cdot \left(h{{\nabla}} u_a\right) merely lies in W^{-1, 2}({{{\Omega}}}) . To overcome this difficulty, we define, for any h\in L^2({{\tilde{\omega}}}) , the coefficient
\begin{equation} \eta_k(h): = \int_{{{{\Omega}}}} h \langle {{\nabla}} u_a\, , {{\nabla}} \psi_k\rangle. \end{equation} | (3.23) |
Since \nabla \cdot \left(h{{\nabla}} u_a\right)\in W^{-1, 2}({{{\Omega}}}) , each of this quantities is well-defined. Furthermore, setting, for any k\in {\rm{I\kern-0.21emN}}\; ,
\begin{equation} \alpha_k(h): = \int_{{{{\Omega}}}} \psi_k \dot u_a \end{equation} | (3.24) |
we have
\begin{equation} \alpha_k(h): = \int_{{{{\Omega}}}} \psi_k \dot u_a = \frac1{\sigma_k}\int_{{{{\Omega}}}} a\langle {{\nabla}} \psi_k\, , {{\nabla}} \dot u_a\rangle = \frac{\eta_k(h)}{\sigma_k}. \end{equation} | (3.25) |
Since \dot u_a\in W^{1, 2}_0({{{\Omega}}}) we have, in W^{1, 2}_0({{{\Omega}}}) , the decomposition
\begin{equation} \dot u_a = \sum\limits_{k = 1}^\infty \frac{\eta_k(h)}{\sigma_k}\psi_k. \end{equation} | (3.26) |
As a consequence, to ensure a decomposition of the form (3.19) it suffices to find, for K large enough, an h\in L^2({{\tilde{\omega}}}) such that
● h is radially symmetric,
● \Vert h\Vert_{L^2({{\tilde{\omega}}})} = 1\, ,
● For any k = 1, \dots, K-1 , \eta_k(h) = 0 ,
● There holds \int_{{{{\Omega}}}}h1_{{\tilde{\omega}}} = \int_{{\tilde{\omega}}} h = 0 .
We define L^2_{{\operatorname{rad}}}(\tilde \omega) as the space of radially symmetric functions in L^2(\tilde \omega) . Let us first note that for any k\in \{1, \dots, K\} the linear maps \eta_k: = h\mapsto \eta_k(h) are continuous on L^2_{{\operatorname{rad}}}(\tilde \omega) . This continuity property is a consequence of the radial symmetry assumption on the coefficients, which from Lemma 11 implies that {{\nabla}} u_a\in L^\infty({{{\Omega}}}) . Indeed, we can then simply write
\begin{align*} \left| \eta_k(h)\right|& \leqslant \Vert {{\nabla}} u_a\Vert_{L^\infty({{{\Omega}}})} \int_{{{{\Omega}}}} \vert h\vert \cdot\vert {{\nabla}} \psi_k\vert \\& \leqslant \Vert {{\nabla}} u_a\Vert_{L^\infty({{{\Omega}}})} \Vert {{\nabla}} \psi_k\Vert_{L^2({{{\Omega}}})}\Vert h\Vert_{L^2({{{\Omega}}})} \\& = \Vert {{\nabla}} u_a\Vert_{L^\infty({{{\Omega}}})} \Vert {{\nabla}} \psi_k\Vert_{L^2({{{\Omega}}})}\Vert h\Vert_{L^2({{\tilde{\omega}}})}, \end{align*} |
whence the continuity.
Defining R_0\in L^2_{{\operatorname{rad}}}({{\tilde{\omega}}})' as
R_0(h): = \int_{{\tilde{\omega}}} h, |
which is obviously continuous on L^2({{\tilde{\omega}}}) we are hence looking for h_k such that
\begin{equation} \Vert h_K\Vert_{L^2({{\tilde{\omega}}})}\, , h_K\in \ker(R_0)\cap \left(\bigcap\limits_{k = 1}^{K-1} \ker(\eta_k)\right). \end{equation} | (3.27) |
However, L^2_{{\operatorname{rad}}}({{\tilde{\omega}}}) is an infinite dimensional Hilbert space, \ker(R_0) and \bigcap_{k = 1}^{K-1} \ker(\eta_k) are closed subspaces of finite co-dimension, hence \ker(R_0)\cap \left(\bigcap_{k = 1}^{K-1} \ker(\eta_k)\right) has finite co-dimension. In particular it is non empty, so there exists h_K such that (3.27) holds. The conclusion follows.
Proof of Theorem III. For the proof of the parabolic Talenti inequalities we follow the main ideas of [29,Theorem 2.1] and of [7]. Since the proof is very similar we mostly present the main steps. To alleviate notations, we simply write u_{a, f} for {{u_{{{\operatorname{parab}}}, a, f}}} . For a fixed a\in \mathcal M({{\Omega}}; T) , we define, for almost every t\in (0;T) , \tilde a(t, \cdot) as the Alvino-Trombetti rearrangement of a(t, \cdot) with respect to {{u_{a, f}}}(t, \cdot) . In other words, for almost every t\in (0;T) and almost every s\in (0;\Vert {{u_{a, f}}}(t, \cdot)\Vert_{L^\infty}) ,
\int_{\{u_{a, f}^*(t, \cdot) \geqslant s\}}\frac1{\tilde a (t, \cdot)} = \int_{\{u_{a, f}(t, \cdot) \geqslant s\}}\frac1{a(t, \cdot)}. |
From [7,Proof of Lemma 1.2] we have, with S_d = d\mathrm{Vol}({{\mathbb B}}(0, 1))^{\frac1d} ,
S_d^2 \mu_{{u_{a, f}}}(t, s)^{2-\frac2d} \leqslant\left(-\frac{d}{ds}\int_{\{u_{a, f}(t, \cdot) > s\}}\frac1{a}\right)\left(-\frac{d}{ds}\int_{\{u_{a, f}(t, \cdot) > s\}} a|{{\nabla}} {{u_{a, f}}}|^2\right). |
Let \overline a is the one-dimensional counterpart of \tilde a (i.e., \tilde a = \overline a(c_d|\cdot|^d) ). Then, as in [7], there holds, almost everywhere,
-\frac{d}{ds}\int_{\{u_{a, f}(t, \cdot) > s\}}\frac1{a} = -\frac{\frac{\partial \mu_{{u_{a, f}}}}{\partial s}}{\overline a}(t, s). |
On the other hand the same arguments as in [3,Proof of Theorem 1] (see also [36]) show that, almost everywhere
-\frac{d}{ds}\int_{\{u_{a, f}(t, \cdot) > s\}}a|{{\nabla}} u_{a, f}|^2 = \int_{\{u_{a, f}(t, \cdot) > s\}}\left(f-\frac{\partial u_{a, f}}{\partial t}\right). |
We can hence conclude that
\begin{equation} S_d^2\overline a \mu_{{u_{a, f}}}(t)^{2-\frac2d} \leqslant -\frac{\partial \mu_{{u_{a, f}}}}{\partial s}(t, s)\int_{\{{{u_{a, f}}} > t\}}\left(f-\frac{\partial {{u_{a, f}}}}{\partial t}\right). \end{equation} | (4.1) |
We now rewrite
\int_{\{{{u_{a, f}}}(t, \cdot) > s\}}\frac{\partial {{u_{a, f}}}}{\partial t} = \int_0^{\mu_{{u_{a, f}}}(t, s)}\frac{\partial u^\#_{a, f}}{\partial t}. |
Introducing as in [29] the function k defined as
k(t, \xi): = \int_0^\xi u^\#(t, \cdot) |
we hence obtain
\begin{equation} \int_{\{{{u_{a, f}}}(t, \cdot) > s\}}\frac{\partial {{u_{a, f}}}}{\partial t} = \frac{\partial k}{\partial t}(t, \mu(t, s)). \end{equation} | (4.2) |
By the Hardy-Littlewood inequality we have
\begin{equation} \int_{\{{{u_{a, f}}} > t\}}f \leqslant \int_0^{\mu_{{{u_{a, f}}}}(t, s)}f^\#. \end{equation} | (4.3) |
Combining these estimates we are left with
\begin{equation} 1 \leqslant -S_d^{-2}\overline a^{-1} \mu_{{u_{a, f}}}(t, s)^{-2+\frac2d} \frac{\partial \mu_{{u_{a, f}}}}{\partial s}(t, s)\left( \int_0^{\mu_{{{u_{a, f}}}}(t, s)}f^\#-\frac{\partial k}{\partial t}(t, \mu(t, s))\right) \end{equation} | (4.4) |
which, after integration, gives
\begin{equation} 0 \leqslant -\frac{\partial^2 k}{\partial \xi^2} \leqslant S_d^{-2}\overline a^{-1} \xi^{-2+\frac2d} \left( \int_0^{\xi}f^\#-\frac{\partial k}{\partial t}(t, \xi)\right). \end{equation} | (4.5) |
We denote by k_* the function obtained by replacing {{u_{a, f}}} by u_{\tilde a, f^*} in the definition of k . Since all the previous inequalities become equalities in this case it follows that the function K: = k-k^* satisfies
\begin{equation} \begin{cases} \frac{\partial k}{\partial t}-S_d^2\overline a\xi^{2-\frac{2}d}\frac{\partial^2 K}{\partial^2 \xi^2} \leqslant 0&{\rm{ in }}\;(0;{{\operatorname{Vol}}}({{\Omega}}))\times (0;T) \\ K(0, \cdot) = 0\, , \\ K(t, 0) = 0 = \frac{\partial K}{\partial \xi}(t, {{\operatorname{Vol}}}({{\Omega}})). \end{cases} \end{equation} | (4.6) |
From the maximum principle, we have K \leqslant 0 , so that the conclusion follows. If j_1 and j_2 are convex non-decreasing functions, the second conclusion of the theorem follows from [4,Proposition 2].
In this paper, we have undertaken the study of certain non-energetic two-phase optimisation problems. Of course, our results are partial, and we now present some open problems that we think are worth investigating.
Open problem I: rearrangements for the time-independent case The first crucial question has to do with the parabolic problem. Indeed, since the Alvino-Trombetti rearrangement we use is defined differently for every time t , the question of time-independent a remains completely open, and we believe it may be fruitful to investigate in the future.
Open problem II: possible relaxations of the problem, bang-bang property for the parabolic optimisation problem The second problem has to do with the conclusion of Theorem I. A more satisfying conclusion that we could not reach would have been a weak bang-bang property, namely that, a profile a\in \mathcal A({{\Omega}}) being given, there exists \tilde a\in \mathfrak M({{\Omega}}) that improves the criterion. Usually, this type of property is obtained using the convexity or concavity of the functional. However, here, what we proved in Theorem II was that the second-order derivative of the functional is positive on an infinite dimensional subspace of the space of admissible perturbations. It is unclear whether this weaker information may be sufficient.
Open problem III: Robin boundary conditions Finally, let us note that, following the recent progresses in the study of Robin Talenti inequalities [5], it may be very interesting to try and understand which type of rearrangement of the weight a may be suitable to obtain Talenti inequalities for two-phases problems under Robin boundary conditions.
I. Mazari was partially supported by the Austrian Science Fund (FWF) through the grant I4052-N32. This work was also partially funded by the French ANR Project ANR-18-CE40-0013 - SHAPO on Shape Optimization and by the Project "Analysis and simulation of optimal shapes - Application to lifesciences" of the Paris City Hall. The author would like to warmly thank the referee for his or her suggestions.
The author declares no conflict of interest.
[1] |
Adnan KMM, Ying L, Sarker SA, et al. (2020) Simultaneous adoption of risk management strategies to manage the catastrophic risk of maize farmers in Bangladesh. GeoJournal 86: 1981–1998. https://doi.org/10.1007/s10708-020-10154-y doi: 10.1007/s10708-020-10154-y
![]() |
[2] |
Fahad S, Wang J, Khan AA, et al. (2018) Evaluation of farmers' attitude and perception toward production risk: Lessons from Khyber Pakhtunkhwa Province, Pakistan. Hum Ecol Risk Assess 24: 1710–1722. https://doi.org/10.1080/10807039.2018.1460799 doi: 10.1080/10807039.2018.1460799
![]() |
[3] |
Oyetunde-Usman Z, Olagunju KO, Ogunpaimo OR (2021) Determinants of adoption of multiple sustainable agricultural practices among smallholder farmers in Nigeria. Int Soil Water Conserv Res 9: 241–248. https://doi.org/10.1016/j.iswcr.2020.10.007 doi: 10.1016/j.iswcr.2020.10.007
![]() |
[4] |
Fahad S, Wang J (2018) Farmers' risk perception, vulnerability, and adaptation to climate change in rural Pakistan. Land Use Policy 79: 301–309. https://doi.org/10.1016/j.landusepol.2018.08.018 doi: 10.1016/j.landusepol.2018.08.018
![]() |
[5] |
Komarek AM, De Pinto A, Smith VH (2020) A review of types of risks in agriculture: What we know and what we need to know. Agric Syst 178: 102738–102747. https://doi.org/10.1016/j.agsy.2019.102738 doi: 10.1016/j.agsy.2019.102738
![]() |
[6] | Hardaker JB, Lien G, Anderson JR, et al. (2004) Coping with risk in agriculture, CAB International. https://doi.org/10.1079/9780851998312.0000 |
[7] | The World Bank (2016) Agricultural Sector Risk Assessment : Methodological Guidance For Practitioners, Washington DC, World Bank Group. |
[8] |
Tudor K, Spaulding A, Roy KD, et al. (2014) An analysis of risk management tools utilized by Illinois farmers. Agric Financ Rev 74: 69–86. https://doi.org/10.1108/AFR-09-2012-0044 doi: 10.1108/AFR-09-2012-0044
![]() |
[9] |
Jankelova N, Masar D, Moricova S (2017) Risk factors in the agriculture sector. Agric Econ (Czech Republic) 63: 247–258. https://doi.org/10.17221/212/2016-AGRICECON doi: 10.17221/212/2016-AGRICECON
![]() |
[10] |
Gebreegziabher K, Tadesse T (2014) Risk perception and management in smallholder dairy farming in Tigray, Northern Ethiopia. J Risk Res 17: 367–381. https://doi.org/10.1080/13669877.2013.815648 doi: 10.1080/13669877.2013.815648
![]() |
[11] |
Akhtar S, LI G cheng, Ullah R, et al. (2018) Factors influencing hybrid maize farmers' risk attitudes and their perceptions in Punjab Province, Pakistan. J Integr Agric 17: 1454–1462. https://doi.org/10.1016/S2095-3119(17)61796-9 doi: 10.1016/S2095-3119(17)61796-9
![]() |
[12] |
Lebel L, Lebel P (2018) Emotions, attitudes, and appraisal in the management of climate-related risks by fish farmers in Northern Thailand. J Risk Res 21: 933–951. https://doi.org/10.1080/13669877.2016.1264450 doi: 10.1080/13669877.2016.1264450
![]() |
[13] |
Mgale YJ, Yunxian Y (2021) Price risk perceptions and adoption of management strategies by smallholder rice farmers in Mbeya region, Tanzania. Cogent Food Agric 7: 1919370. https://doi.org/10.1080/23311932.2021.1919370 doi: 10.1080/23311932.2021.1919370
![]() |
[14] |
Ahmad D, Afzal M, Rauf A (2020) Environmental risks among rice farmers and factors influencing their risk perceptions and attitudes in Punjab, Pakistan. Environ Sci Pollut Res 27: 21953–21964. https://doi.org/10.1007/s11356-020-08771-8 doi: 10.1007/s11356-020-08771-8
![]() |
[15] |
Iqbal MA, Ping Q, Abid M, et al. (2016) Assessing risk perceptions and attitude among cotton farmers: A case of Punjab province, Pakistan. Int J Disaster Risk Reduct 16: 68–74. https://doi.org/10.1016/j.ijdrr.2016.01.009 doi: 10.1016/j.ijdrr.2016.01.009
![]() |
[16] |
Saqib SE, Ahmad MM, Panezai S, et al. (2016) An empirical assessment of farmers' risk attitudes in flood-prone areas of Pakistan. Int J Disaster Risk Reduct 18: 107–114. https://doi.org/10.1016/j.ijdrr.2016.06.007 doi: 10.1016/j.ijdrr.2016.06.007
![]() |
[17] |
Khan I, Lei H, Shah IA, et al. (2020) Farm households' risk perception, attitude and adaptation strategies in dealing with climate change: Promise and perils from rural Pakistan. Land Use Policy 91: 1–11. https://doi.org/10.1016/j.landusepol.2019.104395 doi: 10.1016/j.landusepol.2019.104395
![]() |
[18] |
Huet EK, Adam M, Giller KE, et al. (2020) Diversity in perception and management of farming risks in southern Mali. Agric Syst 184: 1–15. https://doi.org/10.1016/j.agsy.2020.102905 doi: 10.1016/j.agsy.2020.102905
![]() |
[19] |
De Salvo M, Capitello R, Gaudenzi B, et al. (2019) Risk management strategies and residual risk perception in the wine industry: A spatial analysis in Northeast Italy. Land Use Policy 83: 47–62. https://doi.org/10.1016/j.landusepol.2019.01.022 doi: 10.1016/j.landusepol.2019.01.022
![]() |
[20] | Sjöberg L, Moen EB, Rundmo T (2004) Explaining risk perception. An evaluation of the psychometric paradigm in risk, Rotunde. |
[21] |
Begho T, Daubry TP, Ebuka IA (2022) What do we know about Nigerian farmers' attitudes to uncertainty and risk? A systematic review of the evidence. Sci Afr 17: e01309. https://doi.org/10.1016/j.sciaf.2022.e01309 doi: 10.1016/j.sciaf.2022.e01309
![]() |
[22] | Directorate General of Estate Crops (2022) Statistic of National Leading Estate Crops Comodity 2020–2022, Jakarta, Directorate General of Estate Crops, Ministry of Agriculture. |
[23] | Sdoodee S, Rongsawat S (2012) Impact of Climate Change on Smallholders ' Rubber Production in Songkhla Province, Southern Thailand. Proceeeding of The 2012 International and National Conference For The Sustainable Community Development of "Local Community : The Foundation of Development in the ASEAN Economic Community (AEC), 81–86. |
[24] |
Hazir MHM, Kadir RA, Gloor E, et al. (2020) Effect of agroclimatic variability on land suitability for cultivating rubber (Hevea brasiliensis) and growth performance assessment in the tropical rainforest climate of Peninsular Malaysia. Clim Risk Manag 27: 100203–100221. https://doi.org/10.1016/j.crm.2019.100203 doi: 10.1016/j.crm.2019.100203
![]() |
[25] |
Ali MF, Aziz AA, Sulong SH (2020) The role of decision support systems in smallholder rubber production: Applications, limitations and future directions. Comput Electron Agric 173: 105442. https://doi.org/10.1016/j.compag.2020.105442 doi: 10.1016/j.compag.2020.105442
![]() |
[26] |
Abdulla I, Arshad FM (2017) Exploring relationships between rubber productivity and R & D in Malaysia. Outlook Agric 46: 28–35. https://doi.org/10.1177/0030727016689731 doi: 10.1177/0030727016689731
![]() |
[27] |
Qi D, Zhu J, Huang Y, et al. (2021) Factors affecting technology choice behaviour of rubber smallholders: a case study in central Hainan, China. J Rubber Res 1–11. https://doi.org/10.1007/s42464-021-00096-6 doi: 10.1007/s42464-021-00096-6
![]() |
[28] |
Min S, Waibel H, Cadisch G, et al. (2017) The economics of smallholder rubber farming in a mountainous region of Southwest China: Elevation, ethnicity, and risk. Mt Res Dev 37: 281–293. https://doi.org/10.1659/MRD-JOURNAL-D-16-00088.1 doi: 10.1659/MRD-JOURNAL-D-16-00088.1
![]() |
[29] |
Jin S, Huang J, Waibel H (2020) Location and economic resilience in rubber farming communities in southwest China. China Agric Econ Rev 13: 367–396. https://doi.org/10.1108/CAER-06-2020-0153 doi: 10.1108/CAER-06-2020-0153
![]() |
[30] |
Goh HH, Tan KL, Khor CY, et al. (2016) Volatility and market risk of rubber price in Malaysia: Pre- and post-global financial crisis. J Quant Econ 14: 323–344. https://doi.org/10.1007/s40953-016-0037-4 doi: 10.1007/s40953-016-0037-4
![]() |
[31] |
Azis M, Dermoredjo SK, Sayaka B, et al. (2021) Economic perspective of indonesian rubber on agroforestry development, E3S Web Conf 305: 02007. https://doi.org/10.1051/e3sconf/202130502007 doi: 10.1051/e3sconf/202130502007
![]() |
[32] |
Zaw ZN, Sdoodee S, Lacote R (2017) Performances of low frequency rubber tapping system with rainguard in high rainfall area in Myanmar. Aust J Crop Sci 11: 1451–1456. https://doi.org/10.21475/ajcs.17.11.11.pne593 doi: 10.21475/ajcs.17.11.11.pne593
![]() |
[33] | Ismail T, Gohet E (2020) Impact of climate change on latex harvesting, Natural rubber systems and climate change (FTA Working Report), 16–18. |
[34] |
Rivano F, Maldonado L, Simbaña B, et al. (2015) Suitable rubber growing in Ecuador: An approach to South American leaf blight. Ind Crops Prod 66: 262–270. https://doi.org/10.1016/j.indcrop.2014.12.034 doi: 10.1016/j.indcrop.2014.12.034
![]() |
[35] |
Olaniyi ON, Szulczyk KR (2022) Estimating the economic impact of the white root rot disease on the Malaysian rubber plantations. For Policy Econ 138: 102707. https://doi.org/10.1016/j.forpol.2022.102707 doi: 10.1016/j.forpol.2022.102707
![]() |
[36] |
Febbiyanti TR, Fairuza Z (2019) Identification of causes of rubber leaves outbreak in Indonesia. Indones J Nat Rubber Res 37: 193–206. https://doi.org/10.22302/ppk.jpk.v37i2.616 doi: 10.22302/ppk.jpk.v37i2.616
![]() |
[37] |
Jin S, Min S, Huang J, et al. (2021) Falling price induced diversification strategies and rural inequality: Evidence of smallholder rubber farmers. World Dev 146: 105604. https://doi.org/10.1016/j.worlddev.2021.105604 doi: 10.1016/j.worlddev.2021.105604
![]() |
[38] |
Nicod T, Bathfield B, Bosc PM, et al. (2020) Households' livelihood strategies facing market uncertainties: How did Thai farmers adapt to a rubber price drop? Agric Syst 182: 1–11. https://doi.org/10.1016/j.agsy.2020.102846 doi: 10.1016/j.agsy.2020.102846
![]() |
[39] | Penot E, Chambon B, Wibawa G (2017) An history of Rubber Agroforestry Systems development in Indonesia and Thailand as alternatives for a sustainable agriculture and income stability, IRRDB International Rubber Conference, Bali. |
[40] |
Penot E, Thériez M, Michel I, et al. (2022) Agroforestry rubber networks and farmers groups in Phatthalung area in Southern Thailand: A potential for an innovation platform? For Soc 6: 503–526. https://doi.org/10.24259/fs.v6i2.12481 doi: 10.24259/fs.v6i2.12481
![]() |
[41] |
Ali MF, Akber MA, Smith C, et al. (2021) The dynamics of rubber production in Malaysia: Potential impacts, challenges and proposed interventions. For Policy Econ 127: 102449. https://doi.org/10.1016/j.forpol.2021.102449 doi: 10.1016/j.forpol.2021.102449
![]() |
[42] |
Andriesse E, Tanwattana P (2018) Coping with the end of the commodities boom: Rubber smallholders in Southern Thailand oscillating between near-poverty and middle-class status. J Dev Soc 34: 77–102. https://doi.org/10.1177/0169796X17752420 doi: 10.1177/0169796X17752420
![]() |
[43] |
Khan NA, Gao Q, Iqbal MA, et al. (2020) Modeling food growers ' perceptions and behavior towards environmental changes and its induced risks : evidence from Pakistan. Environ Sci Pollut Res 27: 20292–20308. https://doi.org/10.1007/s11356-020-08341-y doi: 10.1007/s11356-020-08341-y
![]() |
[44] |
Ahmad D, Afzal M, Rauf A (2019) Analysis of wheat farmers' risk perceptions and attitudes : evidence from Punjab, Pakistan. Nat Hazards 95: 845–861. https://doi.org/10.1007/s11069-018-3523-5 doi: 10.1007/s11069-018-3523-5
![]() |
[45] |
Lebel P, Whangchai N, Chitmanat C, et al. (2015) Climate risk management in river-based tilapia cage culture in northern Thailand. Int J Clim Chang Strateg Manag 7: 476–498. https://doi.org/10.1108/IJCCSM-01-2014-0018 doi: 10.1108/IJCCSM-01-2014-0018
![]() |
[46] |
Alam A, Guttormsen AG (2019) Risk in aquaculture : Farmers' perceptions and management strategies in Bangladesh. Aquac Econ Manag 23: 359–381. https://doi.org/10.1080/13657305.2019.1641568 doi: 10.1080/13657305.2019.1641568
![]() |
[47] |
Meraner M, Finger R (2017) Risk perceptions, preferences and management strategies: Evidence from a case study using German livestock farmers. J Risk Res 22: 110–135. https://doi.org/10.1080/13669877.2017.1351476 doi: 10.1080/13669877.2017.1351476
![]() |
[48] |
Bishu KG, O'Reilly S, Lahiff E, et al. (2016) Cattle farmers' perceptions of risk and risk management strategies: Evidence from Northern Ethiopia. J Risk Res 21: 579–598. https://doi.org/10.1080/13669877.2016.1223163 doi: 10.1080/13669877.2016.1223163
![]() |
[49] | Thomas G (2018) Risk Management in Agriculture, Scotlandia, The Scottish Parliament. |
[50] |
Durrani H, Syed A, Khan A, et al. (2021) Understanding farmers' risk perception to drought vulnerability in Balochistan, Pakistan. AIMS Agric Food 6: 82–105. https://doi.org/10.3934/agrfood.2021006 doi: 10.3934/agrfood.2021006
![]() |
[51] |
Yarong L, Minpeng C (2021) Farmers' perception on combined climatic and market risks and their adaptive behaviors : a case in Shandong Province. Environ Dev Sustain 23: 13042–13061. https://doi.org/10.1007/s10668-020-01198-8 doi: 10.1007/s10668-020-01198-8
![]() |
[52] |
Onyeneke RU, Nwajiuba CA, Emenekwe CC, et al. (2019) Climate change adaptation in Nigerian agricultural sector: A systematic review and resilience check of adaptation measures. AIMS Agric Food 4: 967–1006. https://doi.org/10.3934/agrfood.2019.4.967 doi: 10.3934/agrfood.2019.4.967
![]() |
[53] |
He R, Jin J, Kuang F, et al. (2020) Farmers' risk cognition, risk preferences and climate change adaptive behavior: A structural equation modeling approach. Int J Environ Res Public Health 17: 85–97. https://doi.org/10.3390/ijerph17010085 doi: 10.3390/ijerph17010085
![]() |
[54] |
Huong NTL, Bo YS, Fahad S (2017) Farmers' perception, awareness and adaptation to climate change: evidence from northwest Vietnam. Int J Clim Chang Strateg Manag 9: 555–576. https://doi.org/10.1108/IJCCSM-02-2017-0032 doi: 10.1108/IJCCSM-02-2017-0032
![]() |
[55] |
Iqbal MA, Ping Q, Zafar MU, et al. (2018) Farm risk sources and their mitigation: A case of cotton growers in Punjab. Pakistan J Agric Sci 55: 683–690. https://doi.org/10.21162/PAKJAS/18.7070 doi: 10.21162/PAKJAS/18.7070
![]() |
[56] |
Le TC, Cheong F (2010) Perceptions of risk and risk management in vietnamese catfish farming: An empirical study. Aquac Econ Manag 14: 282–314. https://doi.org/10.1080/13657305.2010.526019 doi: 10.1080/13657305.2010.526019
![]() |
[57] | Hillson D, Murray-webster R (2005) Understanding and Managing Risk Attitude, United States of America. |
[58] | Pindyck RS, Rubinfeld DL (2013) Microeconomics, United States of America, Pearson. |
[59] | Theuvsen L (2013) Risks and Risk Management in Agriculture. |
[60] | Directorate General of Estate Crops (2019) Tree Crop Estate Statistics of Indonesia (Rubber) 2018–2020, Jakarta, Directorate General of Estates Crops Ministry of Agriculture. |
[61] |
Winarni B, Lahjie AM, Simarangkir BDAS, et al. (2017) Forest gardens management under traditional ecological knowledge in West Kalimantan, Indonesia. Biodiversitas 19: 77–84. https://doi.org/10.13057/biodiv/d190113 doi: 10.13057/biodiv/d190113
![]() |
[62] | Statistics of Kalimantan Barat Province (2022) Kalimantan Barat Province in Figures 2022, Pontianak, Statistics of Kalimantan Barat Province. |
[63] | Moser S, Mußhoff O (2017) Comparing the use of risk-influencing production inputs and experimentally measured risk attitude: Do the decisions of indonesian small-scale rubber farmers match? Ger J Agric Econ 66: 124–139. |
[64] |
Holt CA, Laury SK (2002) Risk aversion and incentive effects. Am Econ Rev 92: 1644–1655. https://doi.org/10.1257/000282802762024700 doi: 10.1257/000282802762024700
![]() |
[65] |
Michailova J, Mačiulis A, Tvaronavičienė M (2017) Overconfidence, risk aversion and individual financial decisions in experimental asset markets. Econ Res 30: 1119–1131. https://doi.org/10.1080/1331677X.2017.1311234 doi: 10.1080/1331677X.2017.1311234
![]() |
[66] | Hosmer DW, Lemeshow S, Sturdivant RX (2013) Applied Logistic Regression., United States of America, Wiley. https://doi.org/10.1002/9781118548387 |
[67] | Pedace R (2013) Econometrics For Dummies, United States of America, John Wiley & Sons Inc. |
[68] | Glorya MJ, Nugraha A (2019) Private Sector Initiatives To Boost Productivity of Cocoa, Coffee, and Rubber in Indonesia, Jakarta. https://doi.org/10.35497/291601 |
[69] | Ruangsri K, Makkaew K, Sdoodee S (2015) The impact of rainfall fluctuation on days and rubber productivity in Songkhla Province. Int J Agric Technol 11: 181–191. |
[70] |
Thaochan N, Pornsuriya C, Chairin T, et al. (2020) Roles of systemic fungicide in antifungal activity and induced defense responses in rubber tree (Hevea brasiliensis) against leaf fall disease caused by Neopestalotiopsis cubana. Physiol Mol Plant Pathol 111: 101511. https://doi.org/10.1016/j.pmpp.2020.101511 doi: 10.1016/j.pmpp.2020.101511
![]() |
[71] |
Go WZ, Chin KL, H'ng PS, et al. (2021) Virulence of rigidoporus microporus isolates causing white root rot disease on rubber trees (Hevea brasiliensis) in Malaysia. Plants 10: 2123. https://doi.org/10.3390/plants10102123 doi: 10.3390/plants10102123
![]() |
[72] |
Siri-udom S, Suwannarach N, Lumyong S (2017) Applications of volatile compounds acquired from Muscodor heveae against white root rot disease in rubber trees (Hevea brasiliensis Müll. Arg.) and relevant allelopathy effects. Fungal Biol 121: 573–581. https://doi.org/10.1016/j.funbio.2017.03.004 doi: 10.1016/j.funbio.2017.03.004
![]() |
[73] |
Syahza A, Backe D, Asmit B (2018) Natural rubber institutional arrangement in efforts to accelerate rural economic development in the province of Riau. Int J Law Manag 60: 1509–1521. https://doi.org/10.1108/IJLMA-10-2017-0257 doi: 10.1108/IJLMA-10-2017-0257
![]() |
[74] | Syarifa LF, Agustina DS, Nancy C, et al. (2016) Low Rubber Prices Impact on Socio Economic Condition of Rubber Smallholders in South Sumatera. Indones J Nat Rubber Res 34: 119–126. |
[75] |
Ullah R, Shivakoti GP, Ali G (2015) Factors effecting farmers' risk attitude and risk perceptions: The case of Khyber Pakhtunkhwa, Pakistan. Int J Disaster Risk Reduct 13: 151–157. https://doi.org/10.1016/j.ijdrr.2015.05.005 doi: 10.1016/j.ijdrr.2015.05.005
![]() |
[76] |
Islam MD Il, Rahman A, Sarker MNI, et al. (2021) Factors influencing rice farmers' risk attitudes and perceptions in bangladesh amid environmental and climatic issues. Polish J Environ Stud 30: 177–187. https://doi.org/10.15244/pjoes/120365 doi: 10.15244/pjoes/120365
![]() |
[77] |
Sun Y, Han Z (2018) Climate change risk perception in Taiwan: Correlation with individual and societal factors. Int J Environ Res Public Health 15: 91. https://doi.org/10.3390/ijerph15010091 doi: 10.3390/ijerph15010091
![]() |
[78] |
Rizwan M, Ping Q, Saboor A, et al. (2020) Measuring rice farmers' risk perceptions and attitude: Evidence from Pakistan. Hum Ecol Risk Assess 26: 1832–1847. https://doi.org/10.1080/10807039.2019.1602753 doi: 10.1080/10807039.2019.1602753
![]() |
[79] |
Fauzi IR, Syarifa LF, Herlinawati E, et al. (2014) Performance of tapping premium system in some rubber plantation enterprises. Indones J Nat Rubber Res 32: 157–180. https://doi.org/10.22302/ppk.jpk.v32i2.162 doi: 10.22302/ppk.jpk.v32i2.162
![]() |
[80] | Oduro-ofori E, Aboagye Anokye P, Acquaye NEA (2014) Effects of education on the agricultural productivity of farmers in offinso municipality. Int J Dev Res 4: 1951–1960. |
[81] |
Paltasingh KR, Goyari P (2018) Impact of farmer education on farm productivity under varying technologies: case of paddy growers in India. Agric Food Econ 6: 7. https://doi.org/10.1186/s40100-018-0101-9 doi: 10.1186/s40100-018-0101-9
![]() |
[82] |
Asadullah MN, Rahman S (2009) Farm productivity and efficiency in rural Bangladesh: The role of education revisited. Appl Econ 41: 17–33. https://doi.org/10.1080/00036840601019125 doi: 10.1080/00036840601019125
![]() |
[83] |
Tongkaemkaew U, Chambon B (2018) Rubber plantation labor and labor movements as rubber prices decrease in Southern Thailand. For Soc 2: 18–27. https://doi.org/10.24259/fs.v2i1.3641 doi: 10.24259/fs.v2i1.3641
![]() |
[84] | Tanakasempipat P (2019) Reuters, Top rubber producer Thailand hit by fungal disease outbreak. Available from: https://www.reuters.com/article/us-thailand-rubber-idUSKBN1X0118. |
[85] |
Razar RM, Hamid NRA, Ghani ZA (2021) GxE effect and stability analyses of selected rubber clones (Hevea brasiliensis) in Malaysia. J Rubber Res 24: 475–487. https://doi.org/10.1007/s42464-021-00115-6 doi: 10.1007/s42464-021-00115-6
![]() |
[86] |
Ouattara PD, Kouassi E, Egbendéwé AYG, et al. (2019) Risk aversion and land allocation between annual and perennial crops in semisubsistence farming: a stochastic optimization approach. Agric Econ (United Kingdom) 50: 329–339. https://doi.org/10.1111/agec.12487 doi: 10.1111/agec.12487
![]() |
[87] |
Djanibekov U, Villamor GB (2017) Market-based instruments for risk-averse farmers: Rubber agroforest conservation in Jambi Province, Indonesia. Environ Dev Econ 22: 133–155. https://doi.org/10.1017/S1355770X16000310 doi: 10.1017/S1355770X16000310
![]() |
[88] |
Sitepu MH, McKay A, Holt RJ (2019) An approach for the formulation of sustainable replanting policies in the Indonesian natural rubber industry. J Clean Prod 241: 118357–118368. https://doi.org/10.1016/j.jclepro.2019.118357 doi: 10.1016/j.jclepro.2019.118357
![]() |
[89] |
Penot E, Ilahang (2021) Rubber Agroforestry Systems (RAS) in West Kalimantan, Indonesia: An historical perspective. E3S Web Conf 305: 02001. https://doi.org/10.1051/e3sconf/202130502001 doi: 10.1051/e3sconf/202130502001
![]() |
1. | Idriss Mazari, A note on the rearrangement of functions in time and on the parabolic Talenti inequality, 2022, 68, 0430-3202, 137, 10.1007/s11565-022-00392-y | |
2. | Dario Mazzoleni, Benedetta Pellacci, Calculus of variations and nonlinear analysis: advances and applications, 2023, 5, 2640-3501, 1, 10.3934/mine.2023059 | |
3. | Idriss Mazari-Fouquer, Optimising the carrying capacity in logistic diffusive models: Some qualitative results, 2024, 393, 00220396, 238, 10.1016/j.jde.2024.02.007 |