Sorghum is a versatile plant with various parts that can be utilized. However, information on the physical and mechanical properties of the sorghum plant is crucial for designing agricultural machinery for primary handling processes such as threshing and chopping. This study aimed to determine the technical characteristics of sorghum plants (Bioguma variety) including the physical and mechanical properties of the stems, leaves, panicles and seeds to design a configuration system concept for threshing and chopping machines capable of processing sorghum plants with high moisture content immediately after harvesting. The study used a descriptive method and included samples of sorghum plants randomly taken from fields in Majalengka and Bogor, Indonesia. The physical and mechanical properties were measured using several replications, ranging from 3 to 30 depending on the parameter. The sorghum plants were harvested at at ages 80, 90 and 108 days after transplanting (DAT). It was found that the moisture content of sorghum stem and seeds decreased with the increase of plant ages where stem ranged between 84.18–79.81 %wb and seeds ranged between 51.7–29.4 %wb. The result revealed that planting ages influenced its properties including stem properties and seed properties. Longer DAT increased the stem hardness from 290.64 ± 29.41 to 350.00 ± 0.81 N and seed hardness from 8.2 ± 1.7 to 44.9 ± 5.4 N but decreased the tensile force of seed form panicles from 16.7 + 3.2 to 6.0 ± 0.8 N. The data on stem strength and seed hardness provide important considerations for the development of several equipment for sorghum processing. The findings of this study can serve as a basis for designing effective and efficient threshing and chopping machines for sorghum plants at high moisture content.
Citation: Ana Nurhasanah, Wawan Hermawan, Tineke Mandang, Astu Unadi, Uning Budiharti, Suparlan, Herry Susanto, Anugerah Fitri Amalia, Teguh Wikan Widodo, Maria Josefine Tjaturetna Budiastuti, Ni Putu Dian Nitamiwati, Diang Sagita, Muhammad Hidayat, Arif Samudiantono. Engineering properties of sorghum bioguma-variety for designing appropriate thresher and chopper machine[J]. AIMS Agriculture and Food, 2023, 8(3): 720-735. doi: 10.3934/agrfood.2023039
[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 |
Sorghum is a versatile plant with various parts that can be utilized. However, information on the physical and mechanical properties of the sorghum plant is crucial for designing agricultural machinery for primary handling processes such as threshing and chopping. This study aimed to determine the technical characteristics of sorghum plants (Bioguma variety) including the physical and mechanical properties of the stems, leaves, panicles and seeds to design a configuration system concept for threshing and chopping machines capable of processing sorghum plants with high moisture content immediately after harvesting. The study used a descriptive method and included samples of sorghum plants randomly taken from fields in Majalengka and Bogor, Indonesia. The physical and mechanical properties were measured using several replications, ranging from 3 to 30 depending on the parameter. The sorghum plants were harvested at at ages 80, 90 and 108 days after transplanting (DAT). It was found that the moisture content of sorghum stem and seeds decreased with the increase of plant ages where stem ranged between 84.18–79.81 %wb and seeds ranged between 51.7–29.4 %wb. The result revealed that planting ages influenced its properties including stem properties and seed properties. Longer DAT increased the stem hardness from 290.64 ± 29.41 to 350.00 ± 0.81 N and seed hardness from 8.2 ± 1.7 to 44.9 ± 5.4 N but decreased the tensile force of seed form panicles from 16.7 + 3.2 to 6.0 ± 0.8 N. The data on stem strength and seed hardness provide important considerations for the development of several equipment for sorghum processing. The findings of this study can serve as a basis for designing effective and efficient threshing and chopping machines for sorghum plants at high moisture content.
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] | Hariprasanna K, Patil J (2015) Sorghum: Origin, classification, biology and improvement. In: Madhusudhana R, Rajendrakumar P, Patil JV (Eds.), Sorghum Molecular Breeding, 3–20. https://doi.org/10.1007/978-81-322-2422-8_1 |
[2] |
Rashwan AK, Yones HA, Karim N, et al. (2021) Potential processing technologies for developing sorghum-based food products: An update and comprehensive review. Trends Food Sci Technol 110: 168–182. https://doi.org/10.1016/j.tifs.2021.01.087 doi: 10.1016/j.tifs.2021.01.087
![]() |
[3] |
Mesfin AH, Girma F (2022) Understanding sorghum farming system and its implication for future research strategies in humid agro-ecologies in Western Ethiopia. J Agric Food Res 10: 100456. https://doi.org/10.1016/j.jafr.2022.100456 doi: 10.1016/j.jafr.2022.100456
![]() |
[4] | Kumari NN, Reddy YR, Blümmel M, et al. (2013) Effect of feeding sweet sorghum bagasse silage with or without chopping on nutrient utilization in deccani sheep. Anim Nutr Feed Technol 13: 243–249. |
[5] |
Ramatoulaye F, Mady C, Fallou S (2016) Production and use sorghum: A literature review. J Nutr Health Food Sci 4: 1–4. https://doi.org/10.15226/jnhfs.2016.00157 doi: 10.15226/jnhfs.2016.00157
![]() |
[6] |
Awad M, Fouda O, Fathy W, et al. (2022) A combined machine for collecting and chopping rice straw. Heliyon 8: e10412. https://doi.org/10.1016/j.heliyon.2022.e10412 doi: 10.1016/j.heliyon.2022.e10412
![]() |
[7] |
El Ghobashy H, Shaban Y, Okasha M, et al. (2023) Development and evaluation of a dual-purpose machine for chopping and crushing forage crops. Heliyon 9: e15460. https://doi.org/10.1016/j.heliyon.2023.e15460 doi: 10.1016/j.heliyon.2023.e15460
![]() |
[8] |
Ndukwu MC, Ejirika C (2016) Physical properties of the African walnut (Tetracarpidium conophorum) from Nigeria. Cogent Food Agric 2: 1232849. https://doi.org/10.1080/23311932.2016.1232849 doi: 10.1080/23311932.2016.1232849
![]() |
[9] |
Li Z, Miao F, Andrews J (2017) Mechanical models of compression and impact on fresh fruits. Compr Rev Food Sci Food Saf 16: 1296–1312. https://doi.org/10.1111/1541-4337.12296 doi: 10.1111/1541-4337.12296
![]() |
[10] |
Ndukwu MC, Ohia A, Anozie O (2019) Influence of moisture content and compression axis on mechanical, physical, and phytochemicals properties of Akuamma (Picralima nitida) fruits and seeds. J Inst of Eng (India): Ser A 100: 417–426. https://doi.org/10.1007/s40030-019-00375-x doi: 10.1007/s40030-019-00375-x
![]() |
[11] |
Mwithiga G, Sifuna MM (2006) Effect of moisture content on the physical properties of three varieties of sorghum seeds. J Food Eng 75: 480–486. https://doi.org/10.1016/j.jfoodeng.2005.04.053 doi: 10.1016/j.jfoodeng.2005.04.053
![]() |
[12] | Gely MC, Pagano AM (2017) Effect of moisture content on engineering properties of sorghum grains. Agric Eng Int: CIGR J 19: 200–209. |
[13] | Simonyan K, El-Okene A, Yiljep Y (2007) Some physical properties of Samaru Sorghum 17 grains. Int Comm Agric Eng 9: 1–15. |
[14] |
Abedi G, Abdollahpour S, Bakhtiari MR (2019) The physical and mechanical properties of potato (Solanum tuberosum L.) tubers as related to the automatic separation from clods and stones. Res Agr Eng 65:77–84. https://doi.org/10.17221/24/2018-RAE doi: 10.17221/24/2018-RAE
![]() |
[15] | Darmajana DA, Hidayat DD, Indriati A, et al. (2022) Comparative study of changes in physical, mechanical, and colour properties between arabica and robusta coffee beans prior to and after roasting. 2493: 040008. https://doi.org/10.1063/5.0109975 |
[16] |
Grundas S, Skubisz G (2008) Physical properties of cereal grain and rape stem. Res Agr Eng 54: 80–90. https://doi.org/10.17221/3/2008-RAE doi: 10.17221/3/2008-RAE
![]() |
[17] |
Sudaryanto A, Hidayat DD, Sagita D, et al. (2022) Engineering properties of the cashew nut in context of designing post-harvest handling and processing machinery. Res Agric Eng 68: 201–209. https://doi.org/10.17221/83/2021-RAE doi: 10.17221/83/2021-RAE
![]() |
[18] | Singh HJ, De D, Sahoo P (2014) Physical properties of soybean cultivated in NEH region of India. Agric Eng Int: CIGR J 16:55–59. |
[19] |
Chattopadhyay P, Pandey K (1999) Mechanical properties of sorghum stalk in relation to quasi-static deformation. J Agric Eng Res 73: 199–206. https://doi.org/10.1006/jaer.1999.0406 doi: 10.1006/jaer.1999.0406
![]() |
[20] |
Al-Hashemi HMB, Al-Amoudi OSB (2018) A review on the angle of repose of granular materials. Powder Technol 330: 397–417. https://doi.org/10.1016/j.powtec.2018.02.003 doi: 10.1016/j.powtec.2018.02.003
![]() |
[21] |
Sagita D (2019) Desain dan konstruksi mesin penyemai benih sayuran portabel tipe vakum untuk pembibitan pada talam semai. J Rekayasa Mesin 10: 265–275. https://doi.org/10.21776/ub.jrm.2019.010.03.7 doi: 10.21776/ub.jrm.2019.010.03.7
![]() |
[22] |
Sutejo A, Sutrisno MS, Wawan H, et al. (2020) Kajian karakteristik fisik, mekanik dan aerodinamik daun teh hasil petikan yang telah dilayukan. J Teknik Pertanian Lampung (J Agric Eng) 9: 171–183. https://doi.org/10.23960/jtep-l.v9i3.171-183 doi: 10.23960/jtep-l.v9i3.171-183
![]() |
[23] | Sriagtula R (2016) Evaluasi produksi, nilai nutrisi dan karakteristik serat galur sorgum mutan brown midrib sebagai bahan pakan ruminansia. Doctoral thesis. Bogor Agricultural University. |
[24] | Sale NA, Muhammed U, Dalha I, et al. (2016) An improved IAR sorghum thresher. Agric Eng Int: CIGR J 18: 119–126. |
[25] |
Lubis A, Mandang T, Hermawan W (2021) Study of the physical and mechanical characteristics of patchouli plants. AIMS Agric Food 6: 525–538. https://doi.org/10.3934/agrfood.2021030 doi: 10.3934/agrfood.2021030
![]() |
[26] | Handayani S (2009) Pengujian sifat mekanik bambu (metode pengawetan dengan boraks). J Teknik Sipil Dan Perencanaan 9: 43–53. |
[27] |
Zein El-den AM, Ahmed SF, Hanafy WM, et al. (2020) Review of some parameters related to the base-cutter of sugarcane harvesters. Misr J Agric Eng 37: 325–330. https://doi.org/10.21608/mjae.2020.47247.1012 doi: 10.21608/mjae.2020.47247.1012
![]() |
[28] |
Seifi MR, Alimardani R (2010) The moisture content effect on some physical and mechanical properties of corn (Sc 704). J Agric Sci 2: 125. https://doi.org/10.5539/jas.v2n4p125 doi: 10.5539/jas.v2n4p125
![]() |
[29] |
El Fawal Y, Tawfik M, El Shal A (2009) Study on physical and engineering properties for grains of some field crops. Misr J Agric Eng 26:1933–1951. https://doi.org/10.21608/mjae.2009.107579 doi: 10.21608/mjae.2009.107579
![]() |
[30] | Gorial B, O'callaghan J (1991) Separation of particles in a horizontal air stream. J Agric Eng Res 49: 273–284. |
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 |