Skip to main content


Front. Chem., 15 May 2019
Sec. Chemical and Process Engineering
Volume 7 - 2019 |

Rational Design of Packaging: Toward Safer and Ecodesigned Food Packaging Systems

Yan Zhu1,2 Bruno Guillemat2 Olivier Vitrac1*
  • 1UMR GENIAL, INRA, AgroParisTech, Université Paris-Saclay, Massy, France
  • 2Pernod Ricard—Breakthrough Innovation Group, Paris, France

It is necessary and urgent to reduce the impact of single-use plastics and in particular those used for beverages and similar products. Fewer plastics or reduced shelf-life need to be explored simultaneously while keeping open the possibilities of bringing on the innovative market designs and consumer experiences. This work details a general computational framework bridging simulation of coupled mass transfer and non-linear thermodynamics with 3D geometry optimization, ecodesign and other mass transfer-related risks (collapse, migration of chemicals). The three-dimensions (impacts, design, shelf-life) of the storage of alcoholic beverages in plastic containers are presented as a representative case-study of the multi-objective ecodesign problem and as an illustration of new trends to accelerate prototyping and innovation. The capacity to evaluate, maximize product shelf-life and minimize packaging waste was tested on 450 miniatures served usually in plane and train. Four alcoholic strengths (including pure water) and four storage conditions were considered. Several alternatives (900 different 3D designs) were explored by simulation and constrained optimization to reach the targeted shelf-life. The large set was used to establish the Pareto optimal design enabling to minimize weights while maximizing shelf-life. The approach is sufficiently versatile and flexible to adapt arbitrary designs and additional phenomena controlling shelf-life (aroma scalping, oxidation, migration of non-intentionally added substances).


During the last decade, the perception of plastics by the public and authorities has fundamentally evolved from beneficial to significantly negative (Hamaide et al., 2014). They have been successfully associated with potential sources of harmful substances, wastes, greenhouse gases, etc. The negative perception has been amplified by the absence of harmonized regulations (see an overview of regulations in the world in Rijk and Veraart, 2010) and the lack of good manufacturing and design practices for materialsused in association in packaging components (laminates, printed and coated surfaces, closure with gaskets or varnishes, plastics associated with paper and board…). Among developed countries, the EU is equipped with the strictest arsenal of rules. EU framework regulation 1935/2004/EC (EC, 2004b) makes a risk assessment and risk management compulsory for the introduction of any new substance, material or industrial practice (active/intelligent packaging system, recycling) regardless the material is in plastics or not. While not mentioning design, EU regulation 2023/2006/EC (EC, 2006) encourages the developments of good manufacturing practices and quality assurance systems at all stages of the production and handling for each of the 17 groups materials and their combinations accepted for food contact. The reduction of the environmental impact of packaging and its wastes has been enforced since 25 years ago in EU with the Directive 94/62/EC (EEC, 1994) and its successive amendments (EC, 2004a, 2005, 2015). Its annex II imposes that any “packaging shall be so manufactured that the packaging volume and weight be limited to the minimum adequate amount to maintain the necessary level of safety, hygiene, and acceptance for the packed product and for the consumer.” Since foods are not specifically mentioned, no association has been made between barrier properties, shelf-life, and amount of plastic materials and finally wastes. Food packaging represents in western countries most of the plastic wastes (Plastics Europe, 2018) with a growing contribution of ready-to-eat foods (Silberbauer and Schmid, 2017). A recent study (Wilson et al., 2017) showed that the willingness of consumers to waste food was affected not only by product size and expiration dates but also by date labels themselves. In simple words, reducing the shelf-life of food products would enable the use of less barrier materials, but will also increase food waste and will reduce access to global markets. FAO (Gustavsson et al., 2011) concluded that well-designed packaging was strongly effective to reduce food losses and waste for all food categories and all along the food supply chain until the final consumer. As a result, the design of a new packaging for a new product or market should be envisioned as a trade-off between several issues: safety, shelf-life and environmental impact.

This study addresses food packaging design as a multicriteria optimization problem, which can be tackled with the tools and concepts of Chemical Engineering used to reduce material and energy requirements (Rangaiah and Bonilla-Petriciolet, 2013). For food packaging, the chief question is how to maximize shelf-life while minimizing packaging weight in worst-case storage and transport conditions (min-max problem). Additional problems of minimization, such as mitigating migration issues of intentionally-added substances, can also be considered during packaging design as shown by Nguyen et al. (2013). The considered degrees of freedom are various, they include the choice of the materials, their formulation and industrial practices along the supply chain; but only the packaging shape and wall thicknesses are common to all mass transfer problems: mass loss, migration issues, aroma scalping, oxidation, etc. The proposed approach is complementary to life cycle assessment (LCA) codified in standard ISO 14040-2006 (ISO, 2006), which can be used to recommend materials, but which lacks enough granularity to orient packaging design. Bottles in polyethylene terephthalate (PET) offer the best illustrative example of environmental challenges, as they represent 7.4 % of all packaging (Plastics Europe, 2018) and have been associated with the growing marine litter in oceans and seas (Jambeck et al., 2015). For comparison, one kilogram of glass is approximately required to produce two bottles of one-liter capacity when 28 bottles can be processed with one kilogram of PET. When PET bottles are produced from recycled material, as authorized in EU regulation 282/2008/EC (EC, 2008a), they offer low environmental footprints (Toniolo et al., 2013; Accorsi et al., 2015; Manzardo et al., 2016). Since the shelf-life of pasteurized, sterilized, and liquid food products tend to be shortened in plastic bottles, the food product and its convenience should be introduced in the holistic approach. Similar conclusions were drawn by Wikström et al. (2014) by noticing that that food production has a higher impact on the environment than the packaging. For single-use beverage bottles, the draft of EU directive 2018/0172/EC (EC, 2018) on the reduction of the impact of certain plastic products on the environment favors product design requirements too restrictive and awareness measures. Original LCA methodology is mainly forensic and retrospective. It has therefore little applicability at early stages of product development when the concepts do not enable a clear determination of material use, energy consumption, and impacts (Holdway et al., 2002). Despite 29 methods available for sustainable product development (Buchert et al., 2017), current global indices remain too rough and too uncertain to finalize design details (Yu et al., 2017). Rezaei et al. (2019) proposed a pairwise preference comparison methodology recently to select the best tradeoff between the consumer acceptability and environment impacts. The approach is qualitative and does not integrate the food explicitly. During the last decades, life cycle design and ecodesign have been promoted by authorities (Keoleian and Menerey, 1993; EC, 2018) and international organizations (ISO, 2002; UNEP, 2009). They are considered as pioneering concepts for food packaging design (Karlsson and Luttropp, 2006; González-García et al., 2016). The methods and tools of ecodesign are, however, still at a primitive stage and can be considered as “vectors of learning, i.e., allowing designers to improve their individual and collective expertise in eco-design” (Vallet et al., 2013). In details, they exhibit some similarities with the TRIZ-problem solving method (Ilevbare et al., 2013) by reusing a knowledge base to self-similar problems: renewable or recycled materials, low material weight, efficient transport, reusability, recyclability or compostability, low waste emissions, etc. Innovative solutions integrating original design, consumer convenience, and new functionalities cannot be explored, because mainly engineering aspects are not explicitly considered such as packaging weight, shape, barrier properties, food shelf-life, cost, safety, etc. By taking the example of alcoholic beverages stored in PET bottles as a representative case study, a tailored and more systematic methodology is proposed to solve, holistically, several mass transfer problems impacting the amount of wastes, the shelf-life of the product and food safety issues. The developments must be seen as an extension of safe-by-design concepts already applied to food packaging (Nguyen et al., 2013). The proposed multiscale modeling, simulation, and optimization framework reuses indeed several mass transfer simulation components of the open-source project FMECAengine (Vitrac, 2018). The chained simulation engine and the procedural solver enable was initially developed to offer concurrent engineering on contamination issues along the supply chain and between packaging components for one to several hundreds of substances. The new extensions focus on packaging geometry optimization in the presence of coupled mass transfer (the diffusion of one substance is affecting the others). Additionally, the new capability to integrate mass transfer equations over three dimensional geometry can be seen complementary to 3D simulation-optimization strategies aiming at maximizing mechanical resistance of bottles with less plastic material (Masood and Keshavamurthy, 2005; Demirel and Daver, 2009; Daver et al., 2012; Hu et al., 2012; Huang et al., 2018). Mass transfer are, indeed, rarely considered in such approaches (Abbès et al., 2010) and never while taking into account the impact on the shelf-life of the food product inside.

The paper is organized as follows. The principles of the description of coupled mass transfer of water and ethanol across polymer walls are described in section theory. Special attention is devoted to the thermodynamical coupling between the chemical potential of water and ethanol in the case of alcoholic beverages. The consequences of mass transfer on the total pressure inside the packaging, the primary cause of leaks and packaging collapse, are also discussed according to packaging filling conditions. Section materials and methods details studied experimental conditions and our strategy of 3D simulation under geometry and shelf-life constraints. The thermodynamical diagrams and the global validation of our predictive approaches are presented in section results and discussion for PET miniatures, which are small plastic bottles containing fruit juices, liquors, or other alcoholic beverages, usually served in trains and planes, and also available in hotel minibars. Comparatively to full-size bottles, their design, and their conditions of storage are very challenging both (i) for the shelf-life of the beverage, as they have the highest surface-to-volume ratios packaging in particularly, and (ii) for the constraints of weight at the board of an aircraft. This example is used to demonstrate that conflicting constraints packaging weight and shelf-life can be managed in a consensual manner if the first is correctly adjusted to the second and if new shapes are introduced. Though only coupled water and ethanol mass transfer from the inside to the outside of the packaging are considered in this study, the approach remains transposable to any permeation, and migration problem met in beverages, pharmaceutical, cosmetic, and biotechnological products. The main findings are summarized in the last section.


Shelf-Life of Alcoholic Beverages

Packaging for long shelf-life products is designed to offer significant barrier properties to gases from the ambiance (water vapor, oxygen) and to food constituents. Aside from oxidation issues in oxygen-sensitive products, the reduction of shelf-life of beverages is primarily caused by mass loss (usually water content) and by the loss of both aroma and sapid compounds. Weight loss occurs as the balance of all mass transfer across container walls and the closure system. As an example, water bottles stored in a dry place are losing weight due to a net permeation rate of water from the inside to the outside. The loss of aroma is usually slower and associated to reversible sorption in the walls, coined “scalping” process (Ducruet et al., 2007; Dombre et al., 2015); and which dominates over mass loss by permeation. Predicting the rates of each mass transfer requires a proper characterization of diffusion and sorption properties in the considered polymer, usually a polyester material such as polyethylene terephthalate (PET), at the different stages of transportation, retailing, storage, and final consumption. In this study, the non-linear behavior of water-ethanol mixtures stored in PET bottles is thought to represent the properties of real alcoholic beverages such as ciders, beer, wines, and spirits realistically. Since only water and ethanol are considered, the type of beverage is defined by its alcoholic strength by volume, denoted abv. It is legally defined in EU as “the ratio of the volume of pure alcohol present in the product in question at 20°C to the total volume of that product at the same temperature” (see Annex I of EU Regulation 110/2008; EC, 2008b). The case of bottled water can be seen as a special case corresponding to abv = 0. Tolerances setting the shelf-life of alcoholic beverages are also relative to the variation of abv during storage and transportation, denoted |Δabv|. EU regulation 1169/2011/EC (EC, 2011) sets out a stringent tolerance of ± 0.3% for non-beer related beverages with abv values larger than 1.2%. Finally, the EU applies a tolerable negative error of 1.5% for the weight of prepackaged liquids, Δw, with volumes equal or larger than 1 L (see Annex I of Directive 76/211/EEC; EEC, 1976). In the context of plastic containers, it should be interpreted as the maximum allowable variation from the weight labeled on the package according to unavoidable variations in weighing, measuring, and mass transfer across packaging walls.

Common Assumptions

Weight and abv tolerances are governed by the mass transfer of both water and ethanol, with possible antagonist effects. For instance, losing water in a dry storage place will concentrate ethanol and increase consequently abv. Losing simultaneously water and ethanol will, by contrast, keep constant abv or will make it decrease. The effects are highly non-linear and require (i) a proper description of the thermodynamics of water-ethanol mixtures during the shelf-life of the beverage and (ii) a predictive description of mutual sorption and diffusion phenomena across the walls of the container. In a real container, the liquid mixture is in equilibrium with a headspace containing air or an inert gas. Its volume and its pressure are complex functions of temperature (thermal expansion of the liquid and gas in fact), of mechanical constraints applied on the semi-flexible packaging, but also of the resulting mass transfer between the liquid and its headspace. This section introduces a comprehensive and consistent description for arbitrary water-ethanol mixtures below the boiling point of ethanol, enabling to predict overpressure and collapse issues. Mass transfer across the walls of the container are described in simplified but not oversimplified calculations, which do not require any steady state for the fluxes for water and ethanol, but which neglect the modifications of the polymer induced by the sorption of water and ethanol. In the case of PET, polymer relaxation effects are triggered by critical temperatures and relative humidities (see the pseudo state diagram 5 in Dubelley et al., 2017), but they do not modify the diffusion coefficients of water in PET significantly (see Figure 14 in Burgess et al., 2014a). Although such non-linear behaviors could have been integrated, it has been preferred to use robust simulations based on worst-case assumptions, which enable the full exploration of the geometry space for containers intended to preserve the properties of alcoholic beverages under strict tolerance constraints (|Δabv| ≤ 0.3% and Δ w ≥ −1.5%), while minimizing the amount of plastic material to produce the container itself.

Thermodynamics of Alcoholic Beverages

Water and ethanol mass transfer, denoted with subscripts w and e, are controlled by two driving forces: partial pressure gradients across walls and the difference of total pressure on both sides of the closure system. By neglecting internal mass transfer in the beverage and the heat transfer required for the vaporization or the condensation of water and ethanol at the liquid-gas interface, the driving potentials on the beverage side are set by the total pressure and the partial pressures in the headspace.

Practical Approximations of Partial Pressures: {pi(T)}i=e,w

In this study, the state of pure-components is chosen as reference state of water and ethanol. As a result, the partial pressures, {pi(T,abv)}i=e,win the headspace, can be inferred from the binary activity coefficients, {γi(T,abv)}i=e,w, relative to the molar fraction of water or of ethanol in the beverage, {xi}i = e, w and from their pure vapor saturation pressures {pi,sat(T)}i=e,w:

pi(T,abv)=xiγi(T,abv)pi,sat(T)for i=e,w    (1)

with T the absolute temperature in Kelvin, but expressed in °C in the description of storage conditions.

Equation (1) is exact when the interactions between water (and also ethanol in a less extent) and the other beverage constituents (pectins, proteins, aroma…) can be neglected. When the real alcoholic beverage is replaced by an equivalent hydroalcoholic solution, one gets:

xe=1-xw for i=e,w    (2)

By noticing that abv is defined at 20°C, molar fractions are related to the density of the mixture with the same composition at the same temperature, denoted ρw+e(20°C,abv):

xe=we/Mewe/Me+(1-we)/Mw            =1Me(ρe20°Cρw+e(20°C,abv)abv)1Me(ρe20°Cρw+e(20°C,abv)abv)+1Mw(1-ρe20°Cρw+e(20°C,abv)abv)    (3)

with we=ρe20°Cρw+e(20°C,abv)abv being the weight fraction in ethanol, {Mi}i = e, w the molecular weights and {ρi20°C}i=e,w the densities of pure components at 20°C.

Practical Approximation of the Total Pressure in the Headspace in Equilibrium With the Beverage.

The total pressure in the headspace Phead results from several equilibria: (i) the chemical equilibrium between the liquid and the headspace, (ii) the thermal equilibrium with the gas and liquid phase and (iii) the balance of forces on the sides of the wall. All equilibria are connected, changing temperature or changing the volume of the headspace affect both the total pressure and the composition of the headspace (and corollary the composition of the liquid). The presence of air or an inert gas affect these equilibria. An efficient procedure to calculate the final total pressure and composition of the headspace has been devised. It takes into account the internal volume of the bottle Vbottle, the initial temperature T0, the initial total pressure P0, the initial composition of the beverage abv0 and the initial headspace volume Vheadt=0. Without a loss of generality, some simplifications are introduced. The deformations of the wall are not considered (the total volume of the gas and liquid is constant), the leaks and mass transfer across the walls are discarded (the equilibration is faster than permeation), the headspace behaves as an ideal gas, air or inert gas cannot dissolve in the liquid and is incondensable.

At any time t, the number of solutes i = w, e in the gas phase, {nig}i=w,e, is derived by combining Equation and the ideal gas law as:

nig=piVheadRT=γi(T,abv)xiPi,sat(T)VheadRT.    (4)

where Vhead is the volume of headspace and R is the ideal gas constant.

Since air (or the inert gas) is assumed to be incondensable, insoluble and initially dry (no water inside), the number of air molecules, na, is given by the initial amount of air in the headspace:

na=P0Vheadt=0RT0.    (5)

Equation (5) and subsequent mass balance on air assumes that nitrogen and oxygen are not exchanged between the beverage and the headspace. The assumption is reasonable with nitrogen (low solubility), but questionable for oxygen. At atmospheric conditions, the solubility of oxygen in water is very low, about eight parts of oxygen per million at 25°C (Truesdale and Downing, 1954). Its solubility increases slightly with ethanol content to reach a value of 3.5% higher than in pure water for abv = 0.15(see Table 1 and Figure 1 of Kutsche et al., 1984). The total desorption or sorption of oxygen in a headspace representing 1:15 of the volume beverage will cause a variation of the total pressure lower than ±0.02 %. Comparatively to the risk associated with ethanol and water, the contribution of oxygen on the risk of collapse was therefore neglected. In the absence of an inert gas, the beverage is additionally assumed to be already at or close to the equilibrium with the earth atmosphere, so that significant desorption, or sorption of oxygen is unlikely.

By assuming that the headspace does not contain ethanol initially before the bottle is filled, the total number of ethanol molecules in the bottle is hence given by:

net=0=ρe20°CMe(Vbottle-Vheadt=0)abv0.    (6)

For water, the number of water molecules in the bottle is derived from the density of the initial water-ethanol mixture, ρw+e(T0,abv0):

nwt=0=ρw+e(T0,abv0)(Vbottle-Vheadt=0)-net=0MeMw.    (7)

For a given volume of the headspace Vhead, the theoretical density of the liquid mixture, ρw+etheoretical(Vhead), is given:

ρw+etheoretical(Vhead)=(net=0-neg)Me+(nwt=0-nwg)MwVbottle-Vhead.    (8)

Equilibrium is reached when the density of the mixture matches the equilibrium density of hydroalcoholic mixture as reported in handbooks and denoted ρw+e(T,abv). Finding Vhead is equivalent to solve the equality:

ρw+etheoretical(Vhead)=ρw+eT(xe(Vhead)).    (9)


xe(Vhead)=net=0-neg(Vhead)net=0+nwt=0-neg(Vhead)-nwg(Vhead).    (10)

Finally, the total pressure in the headspace at equilibrium is given by:

Phead=na+nwg+negVheadRT    (11)

The full procedure reads:

1: calculate na(Eq.7),net=0(Eq.8),nwt=0(Eq.9)2: initialization VheadVheadt=0 at T03: initialization xexet=04:xeoldxet=0, calculate {nig}i=w,e(Eq.6)with the previous guess of xe,1-xe5:(fix point iteration): update xe(Eq.12)6: go to step 4 until |xe-xeold|7: find abv from xe(Eq.3) and read ρw+e(T,abv)8: calculate ρw+etheoretical(Vhead)(Eq.10)9: calculate Δρ=|ρw+e(T,abv)-ρw+etheoretical(Vhead)|10: update Vhead according to the Gold section search method11: go to step 3 until Δρ is minimized12: calculate Phead(Eq.13)    (12)

Transport Equations at Bottle Walls

The internal side of the container (e.g., bottle) is exposed to the vapors of the water-ethanol mixture whereas the external side is exposed only to the relative humidity of the storage place or room at the same temperature. In this description, the beverage is assumed to be macroscopically at thermal equilibrium with the surroundings and that the storage place does not accumulate ethanol (i.e., good ventilation and extraction).

Without distinguishing whether the mass transfer was initiated in the liquid or in the headspace compartment (they are both at thermodynamical equilibrium), the diffusive fluxes at the interfaces in contact with polymer walls are described satisfactorily within the thin film approximation. On the internal side, denoted ∂Ωint and its normal vector nΩint, an effective mass transfer resistance can be considered to account for the vaporization-condensation process associated with the sorption of water and ethanol. By denoting {Cix,y,z,t}i=e,w the solute concentration (SI units in kg·m−3) in the walls at a position (x, y, z) at the time t and {fi,PT}i=e,wthe binary sorption isotherms in the walls, the mass flux densities at the internal boundary condition reads:

[je(T,abv)|Ωintjw(T,abv)|Ωint]=[De(T)Cex,y,z,t|Ωint· nΩintDw(T)Cwx,y,z,t|Ωint· nΩint] =[hepe,sat(T)MeRT(xeγe(T,abv)fe,PT1(Cex,y,z,t|Ωint))hwpw,sat(T)MwRT((1xe)γw(T,abv)fe,PT1(Cwx,y,z,t|Ωint))]    (13)

where {fi,PT-1}i=e,wis the inverse function of the function {fi,PT}i=e,wgiving the amount of solute i absorbed in the walls according to its activity at temperature T; {hi}i = e, w is the equivalent mass transfer conductance across the boundary layer (SI units in m·s−1). The molar fraction of ethanol in the mixture xe is derived from its weight fraction me(t)me(t)+mw(t) (see Equation 3) and the global mass balance in water and ethanol over all internal interfaces of the packing. For a bottle of height H and cylindrical symmetry, the rates of variation of the mass in water and ethanol, denoted mw(t) and me(t), reads:

ddt[me(t)mw(t)]=Ωint[je(T,abv,x,y,z)|Ωintjw(T,abv,x,y,z)|Ωint]dA                         =z=0H[je(T,abv,z)jw(T,abv,z)]2πr(z)dz    (14)

where r(z) is the radial profile of the bottle along the vertical coordinate z.

Transport Equations Within Bottle Walls and Strategy of Resolution

Since the thickness of bottle walls is several magnitude orders smaller than the height of the bottle, mass transfer through bottle walls can be therefore assumed one-dimensional without a significant loss of accuracy. By neglecting the contribution of the polymer relaxation as driving force, the mutual diffusion of water and ethanol in polymer walls reads:

t[Ce(r,z,t)Cw(r,z,t)]=1rr(r[De(T,[Ce,Cw],z)00Dw(T,[Ce,Cw],z)]r[Ce(r,z,t)Cw(r,z,t)])    (15)

where {Di(T,Ci,z)}i=w,e is the diffusion coefficient in the normal direction of the walls. The dependence with the vertical position zand with the local solute concentration accounts for non-uniform drawing and plasticizing effects.

Equation (15) can be solved efficiently for any arbitrary initial solution [Ce(r,z,t=0),Cw(r,z,t=0)] along with boundary condition by decomposing the bottle into vertical sections with similar thickness and transport properties. In other words, mass transfer across container walls can be factorized into n ≥ 1 boundary conditions between the beverage and n equivalent sections coupled via mass balance Equation. By assuming that container walls are at equilibrium with the storage atmosphere at a relative humidity RH0 and a temperature T0, and which is not contaminated by ethanol vapors, the corresponding initial condition is:

[Ce(r,z,t=0)Cw(r,z,t=0)]=[0fw,PT0(RH0)]    (16)

Equations (1–3) (13–6) have been implemented by following the finite-volume formulation of Nguyen et al. (2013) and the opensource project FMECAengine (Vitrac, 2018).

Materials and Methods

Mass Transfer Challenge Test on Real Bottles

Miniature PET bottles with a capacity of 55 mL (injected-mold, supplier Pernod Ricard, France) and a weight of 10.4 ± 0.2 g were used to carry out a long-term challenge test mimicking real storage conditions (up to 7 months) and real alcoholic typical beverages. The bottles were produced by one step injection-blowing. The crystallinity of main bottle walls was determined by the differential scanning calorimetry (model DSC 3, METTLER TOLEDO) and estimated to 32 ± 1.2% (based on three bottles), with an enthalpy of fusion of 140.1 J·g−1 as reported by Wunderlich (2005). The geometry of the bottles and the profile of wall thickness are shown in Figure 1. The thickness profile was measured in more than fifty different random locations per bottles on 10 randomly sampled bottles with the help of a magnetic wall thickness gauge (model MiniTest 7,200 FH with probe FH10, ElektroPhysik GmbH, Germany).


Figure 1. Bottle geometry: (a) radius profile (the arrows indicate the top and bottom position of the sleeve and the vertical line indicate the level of liquid for pure water); (b) 3D-representation assuming a revolution geometry; (c) thickness profile (the shadow indicates the variation range of measurements among 10 bottles); (d) filled and sealed bottle; (e) bottle equipped with a shrink sleeve label.

Four hundred twenty-five flasks were tested in conditions simulating four typical beverages listed in Table 1 and severe storage conditions reported in Table 2. Two temperature levels (35 and 50°C), three levels of relative humidity, RH, (10, 20, and 75%) and four abv levels (0, 0.15, 0.4, 0.7) were covered. Temperatures ± 0.5°C and relative humidities ± 2% were controlled and monitored in three climatic chambers with renewed air. Bottles were filled with 44 ± 0.5 mL at room temperature and atmospheric pressure to keep the headspace volume independent of the hydroalcoholic solution considered. They were subsequently and without delay sealed with 21.5 mm diameter pre-cut foils (aluminum-polyester lid, supplier Embatherm, France) as shown in Figure 1d. A lab-scale pneumatic-controlled heat sealer (model TIME 160, Embatherm, France) was used with a set temperature of 180°C and a set pressure of 200 kPa. Large headspace was preferred to minimize the overpressure and minimize the risk of leaking. Under these specifications, the number of leaking bottles was ranged from 40 (for the most severe condition 1) down to 0 (for pure water stored in all conditions). The initial test was prolonged for beverages B3 and B4 stored in condition S4 by adding a 60 μm thick shrinkable PET sleeve label (supplier Fuji Seal France SAS, France) on remaining bottles as illustrated in Figure 1d. The sleeve did not cover the whole surface area of the bottle (positions are indicated in Figure 1a), but it was sufficient to induce a mass transfer reduction via the thinner surface areas and to confirm the physical description of mass transfer by mutual sorption/permeation across the walls of the bottles.


Table 1. Hydroalcoholic solutions simulating real alcoholic beverages.


Table 2. Studied long-term storage conditions.

Mass transfer was monitored regularly by weighing the whole bottles (bottle, content and sealing) with an electronic precision balance (±0.001 g, PE160, Mettler Toledo, USA). At regular times, some bottles were sacrificed (opened) randomly to enable ethanol content determinations. The contents of the tested bottles were transferred to glass measurement tubes. Alcohol strength by volume of each sample was subsequently determined at three different depths in the tube using an automatic densimeter (model DMA 5000, Anton Paar GmbH, Austria; accuracy of 10−5 % with pure water-ethanol mixtures) equipped with an accurate temperature controller and an autosampler. Ultrapure water was used as reference.

Databases and Properties

The accuracy of the whole calculations is highly sensitive to the quality of thermodynamic data over the range of temperatures and abv values of interest.

Vapor saturation pressures of pure components. Saturation pressures of water and ethanol were approximated with a good accuracy via the modified Goff and Gratch equation (see Equation 17) (Goff and And Gratch, 1946; Goff, 1957) and with the Antoine equation (see Equation 18) (Ddbst, 2018) for water and ethanol, respectively.

With pressures in Pa and temperatures in K, pw,sat(T)reads:

log10pw,sat(T)=-7.90298(Tw,satT-1)+5.02808log10Tw,satT-1.3816·10-7(1011.344(1-TTw,sat)-1)+8.1328·10-3(103.49149*(1-Tw,satT)-1)+log10Pstd    (17)

with all temperatures in Kelvin and the water saturation temperature Tw,sat=373.15 K chosen at the standard pressure Pstd = 101324.6 Pa.

Similarly, pe,sat(T) reads:

log10pe,sat(T)=2.12490+A-BC+T    (18)

with A = 8.20417, B = 1642.89, C = 230.3 for −57 ≤ T ≤ 80°C.

Density of water-ethanol mixtures. Equation (12) requires accurate determinations ofρw+e(T,abv). Values in vacuum have been initially tabulated by OIML (1975) ranging from −20°C up to 40°C, with factual errors discussed in Chanson (2015) to include the revised formula proposed by Bettin and Spieweck (1990). The upper limit of 40°C was imposed by the decrease of the boiling point of ethanol at low pressures. In this study, an extensive database was compiled by assembling the OIML data and by extending them with the predictions of the commercial software AlcoDens (version 3.3, Katmar Software, USA) above 40°C. The slight compressibility of the mixture at atmospheric pressure was based on the compressibility values was corrected from the values reported by the US National Bureau of Standards [see Table 6 of CFR (2006) originated from Osborne et al. (1913)].

Activity coefficients of water and ethanol in hydroalcoholic mixtures. Vapor-liquid equilibria of hydroalcoholic mixtures were calculated from the UNIFAC contribution method (Wittig et al., 2003).

Binary diffusion and sorption isotherms of water and ethanol in PET. Amorphous PET exposed to water (Burgess et al., 2014a) and ethanol (Chandra and Koros, 2009) above 35°C exhibits complex behaviors on long-time scales. Polymer and relaxation and densification are associated to non-local effects and because they are slower than normal diffusion in thick materials, they were neglected in the current study as well as the possible cooperative interactions between water and ethanol. As reported by Burgess et al. (2014a) and Dubelley et al. (2017), Dw(T) was considered independent of the amount of absorbed water and uniform across the wall thickness and along the contour of the bottle. The effect of orientation during blowing was not considered in agreement with observations of Swaroop and Gordon (1980). The sorption isotherm of water was based on a conservative, but realistic, approximation fw,PT0(RH)=RH100Cw,sat(T), where is the saturation concentration Cw,sat(T) of PET. A similar approximation was applied for ethanol based on the experimental values of Ce,sat(T)in soaking experiments. Binary diffusion and sorption properties for all tested conditions are summarized in Table 3 along with their strategy of estimation. Since all properties were inferred from literature or independent measurements, the results were compared to experimental measurements without any adjustment.


Table 3. Thermodynamic and transport properties of water and ethanol in PET.

Simulation Optimization of Bottle Shapes and Geometry

Simulation optimization refers to the optimization of the shape, wall thicknesses according to a set of constraints, including internal capacity, one or several shelf-life criteria. The problem itself is slightly different from algebraic model-based mathematical programming as it couples physical modeling with several optimization loops and calculations of constraints, which may be not algebraic (i.e., not related to algebraic operations, namely, addition, subtraction, multiplication and division). The objective function defining the quantity to maximize or minimize (e.g., the weight of a packaging unit, shelf-life) and constraints (e.g., capacity, shape ratios) need to be calculated iteratively from independent numeric simulations. According to the considered problem, the number of objectives to consider simultaneously can be reduced by treating some objectives by constraints.

We propose a global framework to resolve the coupling between 3D geometries and mutual mass transfer problems (sorption, diffusion, desorption) mathematically. At the expense of adding more solutes (from the food, the packaging or the surrounding), the framework is enough flexible to manage one or multiple linked optimization problems associated to mass transfer phenomena and thermodynamics: shelf-life maximization, packaging weight reduction, minimization of the risk of collapse, minimization of migration issues (for one or several substances). The degrees of freedom are the geometry itself and the shape constraints defined by the designer. As a result, an algorithm executed by a computer can explore iteratively several generations of new designs and provide a much smaller number of feasible solutions meeting goals and adhering closely to several constraints fixed by the end-user, current regulations. The simulation-optimization strategy implemented in this work is illustrated in Figure 2. The global task is split into three independent tasks, denoted [E][D][S]: [E]valuation of mass transfer, [D]ecision, respectively, to targets and [S]olving coupled constraints and goals. Taken individually, the method to resolve each task has been explored and justified by their respective field: chemical engineering, risk, and impact assessment, applied mathematics. As an illustration, mass transfer controlling shelf-life and the safety of food contact materials have been detailed in section two, as well as by Vitrac and Hayert (2005) and by Nguyen et al. (2013). Acceptable thresholds are defined either by regulations or by technological considerations (see an overview in the seminal books of Piringer and Baner, 2008; Robertson, 2016). [S]olving the optimization problem, where some objective functions or constraints are nonlinear, can be carried out by extending well-known techniques used in convex optimization and linear programming techniques, such as the simplex method (Strongin and Sergeyev, 2000; Borwein and Lewis, 2006). Once packaged in a single software or library, the whole [E][D][S] loop is intended be combined with tools commonly used by industry, such as Life Cycle Assessment (ISO, 2006), the data-driven Six Sigma approach (Pyzdek, 2003), the TRIZ method or “the theory of inventive problem solving” (Ilevbare et al., 2013). A notable feature of the approach is to use 3D representations, making it virtually interoperable with computer-aided design and drafting tools (solid and surface modelers), engineering and finite element tools (mechanical resistance of bottles), manufacturing tools (e.g., extrusion blowing of preforms).


Figure 2. Overview of the simulation-optimization framework and its integration with product-design tools used commonly by industry. Extensions are applications, which are covered by the proposed approach but not applied in the presented case-study.

Nonlinear and non-convex problems such as shelf-life and weight optimization problems remain intrinsically more difficult to [S]olve. The optimal solution may be not unique, it may occur at an interior point of the feasible region, at its boundary or at one the extreme points of the feasible region. Integrating uncertainty inherent to any [E]valuation step complicates the definition of a robust termination condition. The best solution lies within a tolerance from the last best optimum and at an acceptable distance of critical constraints to be considered “safe” or “compliant” ([Decision step]). In this work, the internal domain is mapped with a regular grid to offer to enable a sensitivity analysis close to constraints. An instance of the [S]olving step is detailed for a packaging design problem with a prescribed capacity (150 mL). The shape of the bottle is homothetic to the 3D shape presented in Figure 3 (non-optimized) while matching exactly an internal volume of 160 mL for a nominal capacity of Vbeverage = 150 mL and a shelf-life of at least 180 days at 25°C for aa vodka-type beverage (abv0 = 0.4).


Figure 3. Example of initial bottle design to be optimized.

Based on the geometry parameters defined in Figure 3, the weight optimization problem was split into a two-step optimization process: i) setting macroscopic dimensions H, W, D, d, h to match the prescribed internal volume and ii) optimizing lwall, lbottom, ltop, lneck to adjust the barrier properties of the bottle. Without a loss of generality, the general problem is solved using a multiresolution method [see Amaran et al. (2014) for a review of applicable algorithms]. At the lowest level, all parameters in the group are frozen except one. The elementary problem with shelf-life and capacity constraints reads hence:

lwall(X)=argmin[massbottle(X,lwall)with shelf-life(lwall,X,abv0,Vbeverage,T)180 days]where X=[Vbottle,W,D,d,h,ltop,lbottom,lneck]and Vhead=Vbottle-Vbeverage0(19)    (19)

where the shelf-life function is defined as:

shelf-life(lwall,X,abv0,Vbeverage,T)=min[argmin(|abv(t,lwall,X,abv0,Vbeverage,T)|0.003),argmin(mbeverage(t,lwall,X,abv0,Vbeverage,T)1.5%)]    (20)

and where the height of the bottle is found by adding the complementary constraint:

H(Vbeverage,lwall)=arg min [|Vinternal(H,lwall,Vbottle)Vbeverage|]    (21)

As a result, the internal loop manages implicitly design and product constraints in a conventional simulation framework. All input parameters are set conventionally except H and lwall, which are optimized by resolving successively the capacity and the shelf-life problem via a golden-section search method (see section 10.1 in Press, 1992). Additional degrees of freedom controlling the shape of the bottle are introduced at a higher tier by resolving again the elementary problem. The current implementation is sufficiently efficient to screen hundreds or thousands of geometries/designs under constraints. The global optimization when W and D are released has been explored on a 30 × 30 grid (i.e., 900 bottle geometries were optimized).

Results and Discussion

The problem of mass transfer between an alcoholic beverage and the surrounding ambient is described via quantities, which are causally related to the shelf-life of the product: mass loss (Δm) and the residual alcoholic strength (abv). Mass loss may be, however, ambiguous whether it is defined as the net variation of the amount of the consumable product (i.e., legal definition) or as the net variation of the beverage mass plus its packaging (i.e., permeation definition). The second one was used in this study. The difference between the two definitions is related to the amounts of water and ethanol present in the walls. A correction by subtracting the initial mass of the container was avoided as the walls (bottles were stored at 50% RH prior experiments) were also containing water. With this respect, the parameter abv is only associated to the content of the package (destructive measurement), but its variation may evidence ethanol loss (by sorption or permeation), water loss (by sorption or permeation) or more likely a combination of both. Detailed simulations were used to reconstruct the behavior of an average container exposed to controlled conditions. The reconstruction for individual bottles was out of each as they depended on factors out of our control and which are listed in decreasing order of importance: leaks from the lid, volume filling errors, ethanol evaporation before sealing, distribution of PET material in the bottle (thickness variation between bottles), the mass of the bottle.

Driving Forces Controlling Shelf-Life

By noting that PET is not a porous material, water and ethanol mass transfer across bottle walls occur prior dissolution in the polymer and obey to the general sorption-diffusion-desorption model. The main driving forces are the partial pressure differences between the beverage and the surrounding. The variations of binary properties of water-ethanol mixtures with abv between 10 and 70°C: mixture density ρw+e(T,abv), activities {ai(T,abv)}i=e,w, partial pressures{pi(T,abv)}i=e,wand the theoretical pressure of the headspace at constant volume (no air) ptotaltheoretical=pw(T,abv)+pe(T,abv) in equilibrium with the mixture are presented in Figure 4.


Figure 4. Binary properties of water-ethanol mixtures: (A) activities, (B) partial pressures, (C) total pressure in vacuum, (D) liquid density. The symbols depict the tested beverages and storage conditions. Horizontal lines in (A) show the variation of RH when air temperature is increased from 20°C up to 70°C.

By neglecting the amounts transferred to the headspace, binary equilibrium data offer a first approximation of the driving forces based on the alcohol strength (abv) of the considered beverages and the conditions of storage. Activities of water and ethanol do not evolve significantly with temperature. Water activities at 35°C of typical alcoholic beverages correspond to 0.977 (abv = 0.08), 0.960 (abv = 0.15), 0.911 (abv = 0.4) and 0.819 (abv = 0.7) as lower limits for beer, wine, vodka, rum, respectively. In non-tropical regions (for RH < 0.8), alcoholic beverages stored in plastic bottles lose water and therefore weight. Ethanol activities are comparatively much lower 0.153 (abv = 0.08), 0.243 (abv = 0.15), 0.408 (abv = 0.4) and 0.541 (abv = 0.7). The effective driving forces aw(abv)-RH(T)100 and ae(T,abv)-0 across bottle walls tend to be of the same magnitude order in dry conditions. Temperature affects differently the partial pressure of water and ethanol. Increasing the temperature in a closed storage room will cause a dramatic decrease in RH(T) and therefore of water transfer. In our experiments, condition S1 is the most extreme and corresponds to conditions met in container transport. The cumulated pressures ptotaltheoretical offer a lower bound for the total pressure reached in the headspace under vacuum (no air). It can reach up to 25 kPa for liquors at 50°C (point B1(abv = 0.7, ptotaltheoretical = 25 kPa) in Figure 4B).

Bottles are filled with air at atmospheric pressure and at a temperature lower than the storage temperature. These conditions will provoke significant evaporation of ethanol and in a less extent of water, as well as thermal expansion of the mixture (see Figures 4B,C). According to the volume of headspace, Vhead, the pressure difference between the headspace, Phead, and the surrounding, Patm = 100 kPa, can be much greater than ptotaltheoretical. Calculations were carried out for the 55 mL bottles depicted in Figure 1 with a headspace volume varying from 0.5 to 35 mL and by assuming that the walls of the bottles are rigid. When the filling temperature T0 is equal to the storage temperature T, the internal overpressure, PheadPatm, is commensurable to ptotaltheoretical and independent of Vhead. Using storage temperatures greater than T0 provokes strong overpressures strongly dependent on Vhead and on abv. In our experimental conditions Vhead≈11 mL and T0=20°C, Phead reaches 135, 145, 152 kPa at T = 50°C and 116, 121, 123 kPa at T = 35°C for abv = 0.15, 0.4, 0.7 shown in Figure 5, respectively. Vheaddecreases from 11 mL down to 10.43, 9.99, 9.67 mL at T = 50°C and 10.76, 10.51, 10.36 mL at T = 35°C for abv = 0.15, 0.4, 0.7, respectively. Dividing Vheadby 2 and 4 in the worst-case conditions (T = 50°C, abv = 0.7) increases the overpressure by 17.4 and 91%. Negative pressures (below Patm) associated with important risk of collapse occurs in presence of headspace volumes lower than 5 mL and filling temperatures 10°C above the storage temperature (T0>T+10°C).


Figure 5. Calculated iso-headspace pressures (Phead) and iso-headspace volume (Vhead) at equilibrium in the bottle depicted in Figure 1 according to the headspace volume (Vhead0) and temperature (T0) at filling time. The results are presented for two storage temperature 50°C (A–C) and 35°C (D–F) three abv values (0,15, 0.4, 0.7).

A significant total pressure difference across the lid existed for all tested storage conditions and beverages including when the bottles were filled with pure water. For abv > 0.4, the internal pressure was estimated 86 % higher in the presence of ethanol than without and made the risk of leaking through the more likely. Leaking was detected by a rapid mass loss during the week of storage. The number of leaking bottles reached 48% at 50°C and 23% at 35°C for abv≥0.4 (see Table 1). For lowabvvalues (0.15), the number of leaking bottles dropped down to 9.4% and was reduced to zero in the case of pure water.

Coupled Water-Ethanol Mass Transfer From 50 mL Bottles

After discarding leaking bottles, 1889 mass measurements and 286 abv determinations were collected on 358 bottles for all tested conditions (see Table 1). Bottle weights were monitored during a period of up to 5 months with bottles regularly sacrificed to allow abv determinations. The same experimental conditions were simulated using an axisymmetric bottle geometry averaged over ten bottles (see Figure 1) and using an averaged filling volume for each tested condition. Temperature and relative humidities introduced in boundary layers were based on the average of log values. Experimental and simulated mass losses and abv values are reported in Figure 6 with the individual fluxes of water and ethanol reconstructed by simulation detailed in Figure 7. Experimental losses were also standardized to the same filling weight to remove the initial dispersion of data. To facilitate the interpretation of coupled water-ethanol mass transfer, the permeation curve of water alone is shown as reference for all tested conditions. The variations of abv corrected from water loss are also presented by assuming that the permeation of water was proportional to the applied driving forces.


Figure 6. Comparison of experimental (exp) and simulated (sim) mass loss and abv variations for storage conditions S1 (A), S2 (B), S3 (C–E) and beverages B1 (A–C), B2 (D), B3 (E), B4 (A–E). The dashed lines plot the thresholds used to calculate shelf-life (see text). Empty symbols (triangles) represent the theoretical variation of abv when concentration effects due to water permeation are corrected.


Figure 7. Reconstruction of mass losses by sorption and permeation for (A,B) water and (C,D) ethanol.

Weight differences between bottles filled with water and hydroalcoholic mixtures are little and controlled by the differences in water activities, aw(abv)-RH100, on both sides of bottle walls. Though water losses dominated in alcoholic beverages, the shelf-life was determined mainly by the variation of the abv criterion in hydroalcoholic mixtures. The maximum storages times fulfilling weight loss and abv criteria are summarized in Table 4 for the strictest conditions associated to containers larger than 1L and extrapolated down to containers of 50 mL. The rapid permeation of water trough bottle walls caused a concomitant fast increase in ethanol strength (abv) in conditions S1 and S3 (see Figures 6A,C,D). Weight losses exceed faster limits only for mixtures with low abv values (Figure 6E) or when the difference of water activities was low (2b). For both observed criteria (abv and weight loss), a steady regime was reached within 3 to 6 days. Simulations enabled to reconstruct all aspects of mass transfer for all tested conditions and provided a mechanistic interpretation of coupled mass transfer. Figure 7 decomposes mass losses for water and ethanol according to sorption and permeation. Contact times were much longer than permeation lag-times (Figure 7B) and justified the apparent steady state of weight losses and abv variations. The water content in bottle walls evolved in a non-linear manner with desorption dominating in very dry conditions (condition S1). At high (S2) and intermediate (S3) relative humidities, the rapid water sorption is followed by a slow drying when the surrounding RHduring storage is lower than the equilibration RHbefore the experimentation starts (Figure 7A). Ethanol diffuses slowly through PET without reaching a steady state. Lag-times beyond which permeation of ethanol may reach a steady state are evaluated ca. 1.4 and 3.6 years, at 35 and 50°C, respectively.


Table 4. Comparison of shelf-lives extrapolated from experiments and calculated from simulations.

The results presented a significant variability between bottles, which could not be reduced experimentally, and which could not be reproduced by simulation. The sealing was identified as the source of major leaks and defects, but the contribution of a weak gas leak, which could be confused with permeation, could not be tested directly along with other sources of variabilities: variable plasticizing of PET and variable distribution of PET material at injection-blowing time. It is worth noticing that possible heterogeneities in hydrodynamic, temperature and humidity conditions around bottles cannot be invoked to explain variabilities, as all climatic chambers were equipped with an efficient air recirculation and renewal. Further insights were brought by adding a shrinking 60 μm biaxially oriented PET film (sleeve) around the thinnest part of the bottle (cylindrical section of the bottle as seen in Figure 1), which was thought to contribute the most to the overall permeation rate. The sleeve was added to bottles of B3 (abv = 0.15) and B4 (water only) at the end of storage S3 (see Table 2); it covered ca. 80% of the total surface area. Since ethanol could not significantly cross the ca. Eight hundred micrometer hick bottle walls, the external film reduced mainly the flux of water. Mass losses before and after the addition of the sleeve are analyzed individually for 20 bottles containing ethanol (10) or not (10) in Figure 8. The results were not normalized to highlight the different evolutions of the bottles. Before the addition of the sleeve, the “trajectories” followed by the bottles diverged (different permeation rates). Adding a sleeve inflected all permeation rates in a very similar way and made all permeation kinetics looking very similar with almost parallel slopes. This trend confirmed the prevalence of the body of the bottle and of water permeation on overall mass transfer. Since the thickness of the sleeve was commensurable to the thickness fluctuations in the body region (see Figure 1c), the dispersion of permeation rates between bottles was associated to the slight variations of the distribution of weight after injection-blowing (see similar analyses in Figures 3–6, 10 of Daver et al., 2012; Demirel and Daver, 2012, respectively). The flux reduction followed approximately the rules of serial association of mass transfer resistances: (lwall-1+lsleeve-1)-1, with lwall and lsleeve the thicknesses of the wall and the sleeve, respectively.


Figure 8. Modifications of mass transfer when a sleeve is added to the bottles containing the beverage B3 and B4 after several months of storage: mass loss during condition S3 (dash lines, one line per considered bottle) and S4 (continuous lines, one line per considered bottle).

Optimal Design of Packaging Systems for Liquors

The design of a new packaging can be seen either as a problem of a relationship between wall thickness and shelf-life or as a more global problem seeking the maximization of shelf-life while minimizing the mass of plastic material. As previously shown on the tested 55 mL bottles, wall thickness and shelf-life are not independent parameters as soon as a steady regime is reached. Under a strict permeation control, doubling shelf-life requires twice the initial thickness and doubles the weight of plastics. Substantial gain can be achieved without damaging weight only if some non-linarites can be introduced in the original engineering problem. They appear spontaneously if the shape and, in particular, the surface-to-volume ratio is also a degree of freedom. Its contribution was investigated by analyzing the effects of parameters W and D on the weight of bottles matching the design presented in Figure 3. For each value of W and D, the wall thickness of the main sides was optimized to grant a minimum shelf-life of 6 months for a vodka-type product stored in tempered conditions (25°C at 50% RH). The full explored domain covered a 30 × 30 grid describing a uniform variation of both parameters between 30 and 120 mm. All other properties were recovered by simulation-optimization. It is worth noticing that only lwall was updated for each bottle. The other constraints were managed by mathematical constraints applied to the bottle geometry. For efficiency and robustness in particular close to edges, all bottles were modeled as cylindrical B-spline surfaces rather than swept surfaces. All integrations (masses, fluxes) were integrated accordingly. The corresponding space of explored shapes is shown in Figure 9 with optimized values summarized in Figure 10. Shapes vary from almost symmetric cylindrical bottles, flask shapes, long tubes, flat bottles. Some shapes were exotic, but they could be examined at no extra cost by calculations.


Figure 9. Space of geometries explored and optimized for a 160 mL bottle (capacity 150 mL) containing a vodka-type beverage: (A) shapes corresponding to a 30 × 30 combination of W and D, (B) 5 × 5 combination, (C) weights of bottles optimized to fulfill a shelf-life of 180 days at 25°C.


Figure 10. (A) iso-weight and (B) iso-thickness contours of optimized bottles (minimum shelf-life of 180 days) shown in Figure 9. The filled symbols locate the sampled Pareto front. A selection of designs sampled along the Pareto front (from left to right) is presented in (C).

The optimization procedure generates two families of solutions for the same final capacity 150 mL and shelf-life but with weights varying in a ratio 1:2. Flat geometries associated with large W and D values led to heavy bottles and large weights. On the contrary square sections with intermediate W and D values led to very similar weight reductions, while keeping the same shelf-life. The transition between the two regions appears noisy with several maxima. The transition region corresponded to a shift of the shelf-life criterion controlled by abv for heavy bottles and by the mass loss for light bottles. For a given WorDvalue, a unique bottle minimizes the weight of the PET bottle. The continuous curve connecting all points represents the Pareto front for all the criteria considered: capacity, shelf-life, weight, shape factor. Changing Wor Dwill at the expense of either the shape of the bottle or of its weight. A selection of bottles along the Pareto front is shown in Figure 10C. The result can be seen as a shape interpolation process from one shape where the curvature is on the smallest side to a final shape where the curvature dominates the largest side. As the bottle is not symmetric, the transformation is not equivalent to a swap of sides. The optimal shape is closer to the cylinder shape and offers as expected the minimum surface-to-volume ratio.

Conclusions and Perspectives

Minimizing plastic wastes is an urgent priority for our planet. New prototyping tools to explore the consequences (positive or negative) of the packaging design on the shelf-life of sensitive products is a chief priority. It could be though naively than alternative shapes could reduce the weight of packaging dramatically while keeping unchanged shelf-life. For beverages, it is usually not verified as the shapes of commercial bottles tend already to minimize surface-to-volume ratios, at least for those with large capacities. The alternative scenario offers, however, broader perspectives: what is the loss of shelf-life induced by a modification of the format, of the shape or by slimming the walls? The presented computational framework is generic and flexible. It can tackle both design problems, but also analyze the conditions of shelf-life loss.

Integrating simulation and semi-supervised decision making in the ecodesign process is highly desirable. Ecodesign was theorized in the groundbreaking manual of the United Nations Environment Programme (Brezet and Van Hemel, 1997). The strategy is split into seven sub-problems or steps, which could be translated as follows for packaging design: preferring materials with low impacts (S1), minimizing packaging weight (S2), optimizing shape (S3), mitigating transportation and retailing conditions including the risk of overpressure and collapse (S4), updating the packaging format to the food volume consumed (S5), optimizing food shelf-life (S6), optimizing collection, recovery and recycling (S7). Steps S1 and S7 require a global evaluation and should be oriented according to the indices “first impacts” and “end-of-life impacts” common to life cycle analysis. Steps S2 to S5 correspond to a detailed optimization involving the packaging, the food, retailing, and consumer practices or preferences; they can be streamlined to target many optimal or sub-optimal solutions. For instance, the format of the packaging can be optimized by minimizing the amount of packaging waste generated by an average consumer over 1 year (mass of a single packaging unit × number of units consumed in 1 year). The solution should accept shelf-life as a free-parameter bounded by a maximum commercial value and a minimum value combining the retailing time, consumption time and the maximum durability of the product after opening. The different levels of optimization associated to the full problem “less plastics, more shelf-life, exact or better internal capacity, minimum collapse/overpressure, interpolation of the right shape” are resolved iteratively via a three-step scheme, denoted [E][D][S], involving an [E]valuation of mass transfer and thermodynamics; a [D]ecision based on a comparison with acceptable thresholds; and a [S]olving procedure adapted for non-linear problems. Step [E] can be extended to a broad range of phenomena including chemical reactions (oxidation and hydrolysis) and additional mass fluxes (aroma scalping, migration of plastic additives). With this respect, safe-by-design and ecodesign approaches can be treated in parallel without substantial additional cost and by adding rows to Equations (13–16). The only limit is the availability of diffusion and partition coefficients. Step [D] aims at adjusting walls the thickness of the walls, shelf-life or concentrations to fulfill the chosen goals. When many independent phenomena are reducing shelf-life, the latter should be defined as the lowest value imposed by all phenomena considered (e.g., mass loss, scalping, migration, collapse) in the most severe scenario. Finally, the entire geometry can be modified and resolved [S]. Parametric modeling (with many additional geometric constraints) enables to preserve the design intent between dimensions, sections, parts, and assemblies while exploring thousands of combinations. Very few primitives enable to test almost any imaginable bottles in almost real time and to integrate them in the packaging workflow. For example, the approach can resolve what the optimal bottles for given preform are? How to distribute the mass during blow molding.

The detailed case-study demonstrate that very small formats (miniatures) cannot compete with larger PET bottles for storing alcoholic beverages. Similar tolerance on mass loss and alcohol strength by combining two strategies: increasing wall thicknesses and changing their shapes. Additionally, the main stressors have been identified and are associated contrary to intuition to the difference of partial pressure of water across the walls and not the one of ethanol. As the thickening of the wall is beneficial only for water permeation, keeping miniatures in a humid place or even immersed in water can contribute to extending shelf-life without producing additional waste. Beyond brute-force optimization, the [E][D][S] approach can bring a comprehensive description of causes and stimulate the production of alternative designs and industrial practices for alcoholic beverages, but virtually for any packaged food. The gains might be highly substantial. Revised miniatures with triplicate capacity exhibit thus double shelf-life with half-weight. The acceptability of the calculated solutions can be tested in real time in augmented reality approaches and 3D printing solutions.

Future works aim at integrating seamlessly the [E][D][S] iterations within a new multi-purpose open-source project integrating the safe-by-design approaches already developed by INRA (essentially 1D, multiple steps, multiple materials) with these new optimization strategies targeting environmental efficiency and fully 3D. The benefits are expected to be mutual and greater than the sum of the two parts. Mass transfer are integrated over the exact thickness profile of the bottle and surface areas are not anymore overestimated. Shelf-life can be estimated according to various criteria: weight loss, collapse, specific migration. In the case of PET bottles, the capacity to estimate the best ratio of virgin (high impact) and mechanically recycled (low impact) materials will be significant. In this particular case, current EU regulation requires a design and storage conditions, which would ensure a migration lower than 0.017 μg·kg−1 in the beverage for a residual post-consumer contamination of recycled PET conservatively set at 3 mg·kg−1 (Barthélémy et al., 2014). The chief difficulty at this stage is the absence of standardized templates and non-commercial licenses for food packaging design. Bringing the proper amount of plastics to specific food usages, supply chains, and consumer experiences could be accessible while minimizing health and environmental impacts.

Author Contributions

YZ carried out the experiments and ran the simulations. YZ and OV wrote the computational codes and wrote the manuscript. BG reviewed and revised the manuscript. All authors designed the experiments, analyzed and discussed the results.


The study was funded partially by the French National Association of Research and Technology under the grant (N° 2015/1059).

Conflict of Interest Statement

YZ and BG were employed by the group Pernod-Ricard.

The remaining author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


The authors thank Christian Rambaud, Chrystel Bertho, and Laurence Champion of the Pernod Ricard Research Center (CRPR) for their technical assistances.


Abbès, B., Zaki, O., and Safa, L. (2010). Experimental and numerical study of the aging effects of sorption conditions on the mechanical behaviour of polypropylene bottles under columnar crush conditions. Poly. Test. 29, 902–909. doi: 10.1016/j.polymertesting.2010.06.002

CrossRef Full Text | Google Scholar

Accorsi, R., Versari, L., and Manzini, R. (2015). Glass vs. plastic: life cycle assessment of extra-virgin olive oil bottles across global supply chains. Sustainability 7, 2818–2840. doi: 10.3390/su7032818

CrossRef Full Text | Google Scholar

Amaran, S., Sahinidis, N. V., Sharda, B., and Bury, S. J. (2014). Simulation optimization: a review of algorithms and applications. 4OR 12, 301–333. doi: 10.1007/s10288-014-0275-2

CrossRef Full Text | Google Scholar

Barthélémy, E., Spyropoulos, D., Milana, M. R., Pfaff, K., Gontard, N., Lampi, E., et al. (2014). Safety evaluation of mechanical recycling processes used to produce polyethylene terephthalate (PET) intended for food contact applications. Food Addit. Contamin. Part A 31, 490–497. doi: 10.1080/19440049.2013.871755

PubMed Abstract | CrossRef Full Text | Google Scholar

Bettin, H., and Spieweck, F. (1990). A revised formula for the calculation of alcoholometric tables. PTB - Mitteilungen. Forschen und Prüfen 100, 457–463.

Google Scholar

Borwein, J. M., and Lewis, A. S. (2006). Convex Analysis and Nonlinear Optimization: Theory and Examples. New York, NY: Springer Science & Business Media.

Google Scholar

Brezet, H., and Van Hemel, C. (1997). Ecodesign : A Promising Approach to Sustainable Production and Consumption. Paris: United Nations Environment Programme, Industry and Environment, Cleaner Production; Delft: Rathenau Institute; University of Technology.

Google Scholar

Buchert, T., Halstenberg, F. A., Bonvoisin, J., Lindow, K., and Stark, R. (2017). Target-driven selection and scheduling of methods for sustainable product development. J. Clean. Product. 161, 403–421. doi: 10.1016/j.jclepro.2017.05.067

CrossRef Full Text | Google Scholar

Burgess, S. K., Mikkilineni, D. S., Yu, D. B., Kim, D. J., Mubarak, C. R., Kriegel, R. M., et al. (2014a). Water sorption in poly(ethylene furanoate) compared to poly(ethylene terephthalate). Part 2: Kinetic sorption. Polymer 55, 6870–6882. doi: 10.1016/j.polymer.2014.10.065

CrossRef Full Text | Google Scholar

Burgess, S. K., Mikkilineni, D. S., Yu, D. B., Kim, D. J., Mubarak, C. R., Kriegel, R. M., et al. (2014b). Water sorption in poly(ethylene furanoate) compared to poly(ethylene terephthalate). Part 1: Equilibrium sorption. Polymer 55, 6861–6869. doi: 10.1016/j.polymer.2014.10.047

CrossRef Full Text | Google Scholar

CFR (2006). “§§30.66-Gauging manual,” in Title 27-Alcohol, Tobacco Products and Firearms. (U.S. Government Publishing Office).

Chandra, P., and Koros, W. J. (2009). Sorption of lower alcohols in poly(ethylene terephthalate). Polymer 50, 4241–4249. doi: 10.1016/j.polymer.2009.06.066

CrossRef Full Text | Google Scholar

Chanson, E. (2015). Traceability and computerization of alcoholometric tables. OIML Bull. LVI, 3, 1–44.

Google Scholar

Daver, F., Demirel, B., Sutanto, J., and Pang, C. W. (2012). An integrative simulation approach to weight reduction in poly(ethylene terephthalate) bottles. J. Appl. Poly. Sci. 126, E433–E439. doi: 10.1002/app.36780

CrossRef Full Text | Google Scholar

Ddbst (2018). Saturated Vapor Pressure Calculation by Antoine Equation [Online]. Ddbst GmbH. Available online at: (accessed November 18, 2018).

Demirel, B., and Daver, F. (2009). The effects on the properties of PET bottles of changes to bottle-base geometry. J. Appl. Poly. Sci. 114, 3811–3818. doi: 10.1002/app.30990

CrossRef Full Text | Google Scholar

Demirel, B., and Daver, F. (2012). Effects of preform deformation behavior on the properties of the poly(ethylene terephthalate) bottles. J. Appl. Poly. Sci. 126, 1300–1306. doi: 10.1002/app.36925

CrossRef Full Text | Google Scholar

Dombre, C., Rigou, P., and Chalier, P. (2015). The use of active PET to package rosé wine: changes of aromatic profile by chemical evolution and by transfers. Food Res. Int. 74, 63–71. doi: 10.1016/j.foodres.2015.04.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Dubelley, F., Planes, E., Bas, C., Pons, E., Yrieix, B., and Flandin, L. (2017). Water vapor sorption properties of polyethylene terephthalate over a wide range of humidity and temperature. J. Phys. Chem. B 121, 1953–1962. doi: 10.1021/acs.jpcb.6b11700

PubMed Abstract | CrossRef Full Text | Google Scholar

Ducruet, V., Vitrac, O., Saillard, P., Guichard, E., Feigenbaum, A., and Fournier, N. (2007). Sorption of aroma compounds in PET and PVC during the storage of a strawberry syrup. Food Addit. Contam. 24, 1306–1317. doi: 10.1080/02652030701361283

PubMed Abstract | CrossRef Full Text | Google Scholar

EC (2004a). Directive 2004/12/EC of the european parliament and of the council of 11 February 2004 amending Directive 94/62/EC on packaging and packaging waste. Off. J. Eur. Union L 47, 26–31.

EC (2004b). Regulation (EC) no 1935/2004 of the european parliament and of the council of 27 October 2004 on materials and articles intended to come into contact with food and repealing Directives 80/590/EEC and 89/109/EEC. Off. J. Eur. Union L 338, 4–17.

EC (2005). Directive 2005/20/EC of the european parliament and of the council of 9 March 2005 amending Directive 94/62/EC on packaging and packaging waste. Off. J. Eur. Union L 70, 17–18.

EC (2006). Regulation (EC) No 2023/2006 of 22 December 2006 on good manufacturing practice for materials and articles intended to come into contact with food. Off. J. Eur. Union L384, 75–78.

EC (2008a). Regulation (EC) No 282/2008 of 27 March 2008 on recycled plastic materials and articles intended to come into contact with foods and amending Regulation (EC) No 2023/2006. Off. J. Eur. Union L86, 9–18.

EC (2008b). Regulation (EC) No 110/2008 of the european parliament and of the council of 15 January 2008 on the definition, description, presentation, labelling and the protection of geographical indications of spirit drinks and repealing Council Regulation (EEC) No 1576/89. Off. J. Eur. Union L 39, 16–54.

EC (2011). Regulation (EU) No 1169/2011 of the european parliament and of the council. Off. J. Eur. Union L 304:69.

EC (2015). Directive (EU) 2015/720 of the European Parliament and of the Council of 29 april 2015 amending Directive 94/62/EC as regards reducing the consumption of lightweight plastic carrier bags. Off. J. Eur. Union L115, 11–15.

EC (2018). Proposal for A Directive of the European Parliament and of the Council on the Reduction of the Impact of Certain Plastic Products on the Environment. COM/2018/340 final-2018/0172 (COD) (Brussels: European Commission).

EEC (1976). Council directive of 20 January 1976 on the approximation of the laws of the Member States relating to the making-up by weight or by volume of certain prepackaged products. Off. J. Eur. Union Commun. L 46, 1–11.

EEC (1994). European parliament and council directive 94/62/EC of 20 December 1994 on packaging and packaging waste. Off. J. Eur. Union Commun. L 365, 10–23.

Goff, J. A. (1957). Saturation pressure of water on the new Kelvin temperature scale. Trans. Am. Soc. Heat. Ventilat. Eng. 63, 347–354.

Google Scholar

Goff, J. A., and Gratch, S. (1946). Low pressure properties of water from - 160 to 212°F. Trans. Am. Soc. Heat. Ventilat. Eng. 52, 95–122.

Google Scholar

González-García, S., Sanye-Mengual, E., Llorach-Masana, P., Feijoo, G., Gabarrell, X., Rieradevall, J., et al. (2016). “Sustainable design of packaging materials,” in Environmental Footprints of Packaging, ed S. S. Muthu (Singapore: Springer Science+Business Media), 67–74.

Google Scholar

Gustavsson, J., Cederberg, C., Sonesson, U., Van Otterdijk, R., and Meybeck, A. (2011). Global Food Losses and Food Waste – Extent, Causes and Prevention. Rome: Food And Agriculture Organization Of The United Nations (FAO).

Google Scholar

Hamaide, T., Deterre, R., and Feller, J. F. (2014). Environmental Impact of Polymers. London, UK: ISTE Ltd.

Google Scholar

Holdway, R., Walker, D., and Hilton, M. (2002). Eco-design and successful packaging. Des. Manage. J. 13, 45–53. doi: 10.1111/j.1948-7169.2002.tb00330.x

CrossRef Full Text | Google Scholar

Hu, Q. C., Sha, W. J., Li, Y. H., and Wang, Y. S. (2012). Structural optimization and lightweight design of PET bottle based on ABAQUS. Adv. Mater. Res. 346, 558–563. doi: 10.4028/

CrossRef Full Text | Google Scholar

Huang, H.-H., Chen, L.-W., Lu, W.-H., Lin, W.-C., and Chen, Y.-C. (2018). Design and simulation analysis of lightweight HDPE milk bottle. Poly. Poly. Compos. 26, 91–98. doi: 10.1177/096739111802600111

CrossRef Full Text | Google Scholar

Ilevbare, I. M., Probert, D., and Phaal, R. (2013). A review of TRIZ, and its benefits and challenges in practice. Technovation 33, 30–37. doi: 10.1016/j.technovation.2012.11.003

CrossRef Full Text | Google Scholar

ISO (2002). Technical Report ISO/TR 14062. Environmental Management, Integrating Environmental Aspects Into Product Design and Development. Geneva: ISO copiright office.

ISO (2006). International Standard ISO 14040:2006 - Environmental Management-Life Cycle Assessment-Principles and Framework. Geneva: ISO copiright office.

Jambeck, J. R., Geyer, R., Wilcox, C., Siegler, T. R., Perryman, M., Andrady, A., et al. (2015). Plastic waste inputs from land into the ocean. Science 347, 768–771. doi: 10.1126/science.1260352

PubMed Abstract | CrossRef Full Text | Google Scholar

Karlsson, R., and Luttropp, C. (2006). EcoDesign: what's happening? An overview of the subject area of EcoDesign and of the papers in this special issue. J. Clea. Product. 14, 1291–1298. doi: 10.1016/j.jclepro.2005.11.010

CrossRef Full Text | Google Scholar

Keoleian, G. A, and Menerey, D. (1993). Life Cycle Design Guidance Manual: Environmental Requirements and the Product System–Report EPA/600/R-92/226. Washington, DC: U.S. Environmental Protection Agency - Office and Research Development.

Google Scholar

Kutsche, I., Gildehaus, G., Schuller, D., and Schumpe, A. (1984). Oxygen solubilities in aqueous alcohol solutions. J. Chem. Eng. Data 29, 286–287. doi: 10.1021/je00037a018

CrossRef Full Text | Google Scholar

Launay, A., Thominette, F., and Verdu, J. (1999). Water sorption in amorphous poly(ethylene terephthalate). J. Appl. Poly. Sci. 73, 1131–1137.

Google Scholar

Manzardo, A., Loss, A., Mazzi, A., and Scipioni, A. (2016). “Organization life-cycle assessment (olca): Methodological issues and case studies in the beverage-packaging sector,” in Environmental Footprints of Packaging, ed S. S. Muthu (Singapore: Springer Science+Business Media), 67–74. doi: 10.1007/978-981-287-913-4_3

CrossRef Full Text | Google Scholar

Masood, S. H., and Keshavamurthy, V. (2005). Development of collapsible PET water fountain bottles. J. Mater. Process. Technol. 162–163, 83–89. doi: 10.1016/j.jmatprotec.2005.02.176

CrossRef Full Text | Google Scholar

Nguyen, P.-M., Goujon, A., Sauvegrain, P., and Vitrac, O. (2013). A computer-aided methodology to design safe food packaging and related systems. Aiche J. 59, 1183–1212. doi: 10.1002/aic.14056

CrossRef Full Text | Google Scholar

OIML (1975). International Alcoholometric Tables. Paris: International Organization of Legal Metrology.

Osborne, N. S., Mckelvy, E. C., and Bearce, H. W. (1913). Density and thermal expansion of ethyl alcohol and of its mixtures with water. Bull. Bureau Standards 9, 327–474. doi: 10.6028/bulletin.215

CrossRef Full Text | Google Scholar

Piringer, O. G., and Baner, A. L. (2008). Plastic Packaging Interactions With Food and Pharmaceuticals. Weinheim: Wiley-VCH Verlag.

Google Scholar

Plastics Europe (2018). Annual Review 2017–2018. Brussels: PlasticsEurope AISBL.

Press, W. H. (1992). Numerical Recipes in C: The Art of Scientific Computing, 2nd Edn. Cambridge, UK: Cambridge University Press.

Google Scholar

Pyzdek, T. (2003). The Six Sigma Handbook: The Complete Guide for Greenbelts, Blackbelts, and Managers at all Levels, Revised and Expanded Edition. New York, NY: McGraw-Hill.

Google Scholar

Rangaiah, G. P., and Bonilla-Petriciolet, A. (2013). Multi-Objective Optimization in Chemical Engineering: Developments and Applications. Chichester, UK: John Wiley & Sons.

Google Scholar

Rezaei, J., Papakonstantinou, A., Tavasszy, L., Pesch, U., and Kana, A. (2019). Sustainable product-package design in a food supply chain: a multi-criteria life cycle approach. Pack. Techn. Sci. 32, 85–101. doi: 10.1002/pts.2418

CrossRef Full Text | Google Scholar

Rijk, R., and Veraart, R. (eds.). (2010). Global Legislation for Food Packaging Materials. Weinheim: WILEY-VCH Verlag GmbH & Co. KGaA.

Robertson, G. L. (2016). Food Packaging: Principles and Practice, 3rd Edn. Boca-Raton, FL: CRC Press. doi: 10.1201/b21347

CrossRef Full Text | Google Scholar

Silberbauer, A., and Schmid, M. (2017). Packaging concepts for ready-to-eat food: recent progress. J. Pack. Techn. Res. 1, 113–126. doi: 10.1007/s41783-017-0019-9

CrossRef Full Text | Google Scholar

Strongin, R. G., and Sergeyev, Y. D. (2000). Global Optimization with Non-Convex Constraints: Sequential and Parallel Algorithms. Dordrecht: Springer Science & Business Media.

Google Scholar

Swaroop, N., and Gordon, G. A. (1980). Transport in oriented bottles made from polyester resin. Poly. Eng. Sci. 20, 78–81. doi: 10.1002/pen.760200113

CrossRef Full Text | Google Scholar

Toniolo, S., Mazzi, A., Niero, M., Zuliani, F., and Scipioni, A. (2013). Comparative LCA to evaluate how much recycling is environmentally favourable for food packaging. Res. Conserv. Recyc. 77, 61–68. doi: 10.1016/j.resconrec.2013.06.003

CrossRef Full Text | Google Scholar

Truesdale, G. A., and Downing, A. L. (1954). Solubility of oxygen in water. Nature 173, 1236–1236. doi: 10.1038/1731236a0

CrossRef Full Text | Google Scholar

UNEP (2009). Design for Sustainability: A Step-by-Step Approach. Paris: United Nations Environment Programme, DTIE, Sustainable Consumption and production Branch.

Vallet, F., Eynard, B., Millet, D., Mahut, S. G., Tyl, B., and Bertoluci, G. (2013). Using eco-design tools: an overview of experts' practices. Design Stud. 34, 345–377. doi: 10.1016/j.destud.2012.10.001

CrossRef Full Text | Google Scholar

Vitrac, O. (2018). FMECAengine: FMECA Software Developed in the Framework of the Project SafeFoodPack Design [Online]. Available online at: (accessed November 22, 2018).

Vitrac, O., and Hayert, M. (2005). Risk assessment of migration from packaging materials into foodstuffs. AIChE J. 51, 1080–1095. doi: 10.1002/aic.10462

CrossRef Full Text | Google Scholar

Wikström, F., Williams, H., Verghese, K., and Clune, S. (2014). The influence of packaging attributes on consumer behaviour in food-packaging life cycle assessment studies - a neglected topic. J. Clean. Product. 73, 100–108. doi: 10.1016/j.jclepro.2013.10.042

CrossRef Full Text | Google Scholar

Wilson, N. L. W., Rickard, B. J., Saputo, R., and Ho, S.-T. (2017). Food waste: The role of date labels, package size, and product category. Food Qual. Prefer. 55, 35–44. doi: 10.1016/j.foodqual.2016.08.004

CrossRef Full Text | Google Scholar

Wittig, R., Lohmann, J., and Gmehling, J. (2003). Vapor–liquid equilibria by unifac group contribution. 6. revision and extension. Indus. Eng. Chem. Res. 42, 183–188. doi: 10.1021/ie020506l

CrossRef Full Text | Google Scholar

Wunderlich, B. (2005). Thermal Analysis of Polymeric Materials. Berlin: Springer-Verlag.

Google Scholar

Yu, X., Zhang, H., Shu, H., Zhao, W., Yan, T., Liu, Y., et al. (2017). A robust eco-design approach based on new sensitivity coefficients by considering the uncertainty of LCI. J. Adv. Manufact. Syst. 16, 185–203. doi: 10.1142/S0219686717500123

CrossRef Full Text | Google Scholar

Keywords: ecodesign, packaging, mass transfer, alcoholic beverages, shelf-life, polyethylene terephthalate, thermodynamics

Citation: Zhu Y, Guillemat B and Vitrac O (2019) Rational Design of Packaging: Toward Safer and Ecodesigned Food Packaging Systems. Front. Chem. 7:349. doi: 10.3389/fchem.2019.00349

Received: 01 December 2018; Accepted: 29 April 2019;
Published: 15 May 2019.

Edited by:

Thomas Karbowiak, Agrosup Dijon, France

Reviewed by:

Juan Andrés Cárcel, Universitat Politècnica de València, Spain
Masoud Ghaani, University of Milan, Italy
Sara Limbo, University of Milan, Italy

Copyright © 2019 Zhu, Guillemat and Vitrac. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Olivier Vitrac,