We propose a variational framework for the resolution of a
non-hydrostatic Saint-Venant type model with bottom topography. This model is a shallow
water type approximation of the incompressible Euler system with free
surface and slightly differs from the Green-Nagdhi model,
see [13] for more details about the model derivation.
The numerical approximation relies on a prediction-correction type scheme initially introduced by Chorin-Temam [17] to treat the incompressibility in the Navier-Stokes equations. The hyperbolic part of
the system is approximated using a kinetic finite volume solver
and the correction step implies to solve a mixed problem where the velocity and the pressure are defined in compatible finite element spaces.
The resolution of the incompressibility constraint leads to an elliptic problem involving the non-hydrostatic part of the pressure. This step uses a variational formulation of a shallow water version of the incompressibility condition.
Several numerical experiments are performed to confirm the relevance of our approach.
Citation: Nora Aïssiouene, Marie-Odile Bristeau, Edwige Godlewski, Jacques Sainte-Marie. A combined finite volume - finite element scheme for a dispersive shallow water system[J]. Networks and Heterogeneous Media, 2016, 11(1): 1-27. doi: 10.3934/nhm.2016.11.1
Related Papers:
[1]
Xavier Blanc, Claude Le Bris .
Improving on computation of homogenized coefficients in the
periodic and quasi-periodic settings. Networks and Heterogeneous Media, 2010, 5(1): 1-29.
doi: 10.3934/nhm.2010.5.1
[2]
Martin Heida, Benedikt Jahnel, Anh Duc Vu .
Regularized homogenization on irregularly perforated domains. Networks and Heterogeneous Media, 2025, 20(1): 165-212.
doi: 10.3934/nhm.2025010
[3]
Junlong Chen, Yanbin Tang .
Homogenization of nonlinear nonlocal diffusion equation with periodic and stationary structure. Networks and Heterogeneous Media, 2023, 18(3): 1118-1177.
doi: 10.3934/nhm.2023049
[4]
Fabio Camilli, Claudio Marchi .
On the convergence rate in multiscale homogenization of fully nonlinear elliptic problems. Networks and Heterogeneous Media, 2011, 6(1): 61-75.
doi: 10.3934/nhm.2011.6.61
[5]
Patrick Henning .
Convergence of MsFEM approximations for elliptic, non-periodic homogenization problems. Networks and Heterogeneous Media, 2012, 7(3): 503-524.
doi: 10.3934/nhm.2012.7.503
[6]
Ben Schweizer, Marco Veneroni .
The needle problem approach to non-periodic homogenization. Networks and Heterogeneous Media, 2011, 6(4): 755-781.
doi: 10.3934/nhm.2011.6.755
[7]
Grigor Nika, Adrian Muntean .
Hypertemperature effects in heterogeneous media and thermal flux at small-length scales. Networks and Heterogeneous Media, 2023, 18(3): 1207-1225.
doi: 10.3934/nhm.2023052
[8]
Mogtaba Mohammed, Mamadou Sango .
Homogenization of nonlinear hyperbolic stochastic partial differential equations with nonlinear damping and forcing. Networks and Heterogeneous Media, 2019, 14(2): 341-369.
doi: 10.3934/nhm.2019014
[9]
Rémi Goudey .
A periodic homogenization problem with defects rare at infinity. Networks and Heterogeneous Media, 2022, 17(4): 547-592.
doi: 10.3934/nhm.2022014
[10]
Tasnim Fatima, Ekeoma Ijioma, Toshiyuki Ogawa, Adrian Muntean .
Homogenization and dimension reduction of filtration combustion in heterogeneous thin layers. Networks and Heterogeneous Media, 2014, 9(4): 709-737.
doi: 10.3934/nhm.2014.9.709
Abstract
We propose a variational framework for the resolution of a
non-hydrostatic Saint-Venant type model with bottom topography. This model is a shallow
water type approximation of the incompressible Euler system with free
surface and slightly differs from the Green-Nagdhi model,
see [13] for more details about the model derivation.
The numerical approximation relies on a prediction-correction type scheme initially introduced by Chorin-Temam [17] to treat the incompressibility in the Navier-Stokes equations. The hyperbolic part of
the system is approximated using a kinetic finite volume solver
and the correction step implies to solve a mixed problem where the velocity and the pressure are defined in compatible finite element spaces.
The resolution of the incompressibility constraint leads to an elliptic problem involving the non-hydrostatic part of the pressure. This step uses a variational formulation of a shallow water version of the incompressibility condition.
Several numerical experiments are performed to confirm the relevance of our approach.
1.
Introduction
Tequila is an alcoholic beverage obtained from Agave tequilana Weber blue variety; its production is protected under a well-established region in Mexico by the Official Mexican Norm for Tequila [1]. The production process of Tequila features five stages: cooking, grinding, fermentation, distillation and aging or resting [2]. Usually, research about Tequila production is focused in each stage of the process, and also in how to increase the Tequila's yield and quality. Therefore, fast analytical techniques have been developed to ensure the authenticity and quality of Tequila to satisfy the Mexican regulation [3]. Commercial regulations enforce that the maximal permissible concentration of higher alcohols must be 500 mg per 100 ml of anhydrous alcohol. In addition, Tequila aroma is strongly related to these higher alcohols and volatile compounds. Among the 200 volatile compounds, only a few are essential in the quality and bouquet of Tequila [4]. Therefore, the milestone stage for the process of Tequila production relies in the fermentation stage; in this stage the sugars (principally fructose) are metabolized into ethanol, carbon dioxide (CO2), and also into a wide variety of secondary products such as the volatile compounds which are very important in the bouquet of the liquor of Tequila [5]. The composition of Tequila can be affected by a wide range of factors from which the fermented most and the yeast metabolism are the most important. For the production of higher alcohols, the composition of Tequila liquor can be affected by the nitrogen source (amount and type) and how the yeast will metabolize it [5,6,7].
Saccharomyces cerevisiae is a versatile microorganism commonly used to ferment agave juice in the production of Tequila and also employed in many different industrial applications. This microorganism has been investigated in metabolic studies which elucidate its gene functions, integration, and metabolism [8]. Moreover, key research areas are focused in the fermentation process by using novel strains such as Saccharomyces and non-Saccharomyces [9,10].
Genomic Scale Models (GSM) are used to estimate the metabolic fluxes distribution; these models contain a wide collection of stoichiometrically-branded biochemical reactions related to enzymes in the cell/tissue [11]. Applications of GSM go from topological networks analyses, phenotype behaviors until the simulation of metabolism under certain metabolic engineering strategies [8]. Thus, over 100 GSMs have been developed for a wide variety of microorganisms such as bacteria, eukaryotes and archea i.e., Haemophilus influenzae, Escherichia coli, Saccharomyces cerevisiae, Helicobacter pylori, Staphyloccocus aureus, Bacillus subtilis, Homo sapiens, Pseudomona aeruginosa and for the genus Synechocystis, among others [12,13,14]. By taking advantage of the capabilities described in the literature, it is possible to test novel hypotheses on metabolic functions for a set of organisms-of-interest [14].
In the last decade, GSMs for Saccharomyces cerevisiae have been reported, such as the models iFF708, iND750, iLL672, iIN800, iMM904, iTO977 and Yeast 6 [13,15,16,17,18,19,20]. These models have been developed gradually by adding open reading frames (ORF) reactions and compartments. The first version of the model iFF708 includes 708 ORF with 1175 biochemical reactions and 3 compartments. In contrast, the last version, Yeast 6, contains 900 ORF with 1888 biochemical reactions and 15 compartments [8]. As mentioned before, there are already models which their number in biochemical reactions is around one thousand, however, using these models in practical applications represents a challenge due to their complexity. Moreover, computational algorithms for FBA have been commonly used in different strain models [12,21]. The metabolic phenotype of the microorganism is usually analyzed based on the flux distributions in the metabolic network, for example, for S. cerevisiae, it was found that the global cellular functions under aerobic and anaerobic culture were consistent with the experimental data [22]. Pereira et. al. (2016) selected four GSMs to analyze S. cerevisiae; iFF708; IMM904; iTO977 and Yeast 6, however, a comparison among these models was performed in order to find the best approximation with regard to experimental fluxes in vivo resulting that iFF708 showed the best approximation in terms of the metabolic flux distribution for the central metabolism. To analyze the GSM, FBA is used to estimate the flux distribution in the metabolic network. This tool has demonstrated to be a key alternative to foresee the metabolic-networks capabilities especially when there is a scarcity of kinetic parameters [23,24]. However, for practical applications it is helpful to have more simple models that can help in the development of a robust control system for alcoholic fermentation. To our understand, there is no information available with regard to the metabolic flux estimation for S. cerevisiae in the flux distribution of higher alcohols in Tequila production which, traditionally, it is performed using a batch process. Usually, studies in this area focus on the fermentative capabilities of the microorganism and the yields of the synthesis of volatiles compounds [3,5,7,9,10,25,26].
Therefore, in this research a GSM is reconstructed and validated comparing FBA results with experimental data for the continuous fermentation process of Tequila [5,22] using Saccharomyces cerevisiae. This model is tested under aerobic and anaerobic conditions considering the central metabolism which features the catabolic and anabolic biochemical reactions involved in the synthesis of higher alcohols.
2.
Materials and methods: mathematical developments
2.1. Flux balance analysis, FBA
FBA is a methodology used to analyze GSM applying thermodynamic restrictions as well as network and transport capabilities [21,27]. The distribution of metabolic fluxes is found by optimizing an objective function which in the present research is the microorganism growth (μ) [1/s]. In addition, an optimal flux distribution is calculated using a linear programming technique under steady state conditions. According to the the optimal flux distributions under the set of considerations or restrictions, it is possible to generate a quantitative hypothesis in silico which can be tested experimentally [21].
2.1.1. Mathematical representation
According to the biochemical information, the mathematical representation of the GSM is obtained from a mass balance for each metabolite in the biological system and it is depicted as a matrix as shown in Eq (2.1)
dxdt=S⋅v,
(2.1)
where, x [mmol/g D.W.] is the concentration vector of all metabolites considered in the metabolic network, S is the stoichiometric matrix and v [mmol/g D.W. h−1] is the internal flux vector and the whole exchange fluxes considered between the organism and the culture medium. By considering steady-state, the above equation becomes
S⋅v=0.
(2.2)
The system is solved by applying a linear programming technique that maximizes the microorganism growth function which has been experimentally corroborated [21]. The growth function z(xm) is defined by the microorganism biomass as shown by Eq (2.3)
z(xm)=Growth∑m=1dm⋅xm→biomass,
(2.3)
where dm represents each metabolite proportion and xm is the metabolite in the biomass composition.
Applying thermodynamic restrictions as well as the whole interchange fluxes related to the transport capabilities as shown in Eqs (2.4) and (2.5), the flux distribution is obtained
vj≥0,
(2.4)
α≤bj≤βj,
(2.5)
where α and β represent the thermodynamic limits, low and superior, respectively, for each biochemical reaction of the GSM. Furthermore, bj is the estimated value by the FBA or a real value which could be specified as a restriction in the optimization problem.
2.1.2. Reconstruction of the genomic scale model
The reconstructed stoichiometric model took advantage of other models reported in the literature to analyze the metabolism of S. cerevisiae [8,22,28,29,30,31], and it is reported in Appendix A. The model considers the central metabolism for glycolysis (reactions 1-8), the pentose phosphate pathway (reactions 17-22), and the Krebs cycle (reactions 23-31), as well as the fermentative paths for acetate, glycerol and ethanol production (reactions 9-12, 15 and 16), under anaerobic conditions. For the metabolic production of acetate and ethanol, other two key reactions 13 and 14 are included. Therefore, the aim of our GSM model is that it can be used in practical applications as the Tequila production. Usually, there are simplified models which can be aerobic or anaerobic even though they do not include the metabolic and anabolic pathways related with the synthesis of characteristic volatile compounds like the higher alcohols in Tequila production. Therefore, this last characteristic as well as the small number of biochemical reactions to the complex process-production of Tequila are the main features with respect to other models that justify why the GSM is reconstructed. Additionally, our GSM model includes the catabolic pathway which is also called Erlych's pathway and the anabolic pathway that is result of the metabolization of amino acids like Val, Leu, Ileu, Thre and Phe.
2.1.3. Growth objective function
The objective function is defined as the biomass composition and it is considered in the biochemical reaction 65 (flux) for the macromolecules composition such as proteins, carbohydrates, lipids, deoxyribonucleic acid (DNA), and ribonucleic acid (RNA) using the cited values by Nissen et al. (1997) [30]. Under continuous culture fermentation, the cellular biomass composition varies according to the rate dilution (D) [h−1]. The main variation in biomass composition is found in proteins and RNA that increases linearly as a function of D and carbohydrates (see Table S1). In this research, the biomass composition at D = 0.1 h−1 was taken from the literature [30]. The model encompasses both the aerobic and anaerobic conditions where reactions 12, 14, 16, 26 and 33 take place in aerobic conditions while 13, 15 and 27 happen in anaerobic conditions. In order to synthesize different macromolecules present in S. cerevisiae, the polymerization energy reported by Stephanopolous (1998) [29] was considered for reaction 65.
2.1.4. Biochemical protein synthesis
Protein synthesis considers 20 amino acids (reaction 63), 15 of them feature only one amino acid reaction (reactions 36-49) [29]. In particular, the amino acid syntheses; leucine (LEU), valine (VAL), threonine (THR), isoleucine (ILEU) and phenylalanine (PHE), are considered in two reactions, reactions 81 and 91 [28]. The first one uses an amino acid as a precursor to the formation of α-keto acid and the second one from the α-keto acid produces the corresponding amino acid. Those amino acids are strongly related to the synthesis of higher alcohols, particularly, they are a nitrogen source which is presented in the complex mediums of agave juice used for Tequila production.
2.1.5. Higher-alcohol synthesis
The synthesis of higher alcohols is strongly related to amino acids through α-keto acids which are precursors of n-propanol, isoamyl and amyl alcohol, 2-phenyl-ethanol and isobutanol. Higher alcohols production reactions 70-79 can be portrayed by two main pathways as:
1) Catabolic pathway, also known as the Ehrlich pathway, considers the transamination and deamination of the amino acids in the medium culture in reactions 81-92.
2) Anabolic pathway allows the α-keto acid synthesis from the present sugar in the medium. Especially, from the α-keto acids, the higher alcohols are obtained with a similar mechanism as the catabolic pathway in reactions 71-80.
2.1.6. DNA and RNA synthesis
For the syntheses of RNA (reaction 65) and DNA (reaction 66) macromolecules, the formation of molecules for the nucleotides biosynthesis is considered. For RNA, adenosine-monophosphate (AMP), guanosine-monophosphate (GMP), cytidine-monophosphate (CMP) and uridine-monophosphate (UMP) are considered in reactions 50-53 while for DNA, deoxyadenosine monophosphate (dAMP), deoxyguanosine monophosphate (dGMP), deoxycytidine monophosphate (dCMP), and deoxyuracil monophosphate (dUMP) are contemplated in reactions 54-57 [29].
2.1.7. Carbohydrates synthesis
Carbohydrate composition for S. cerevisiae includes glycogen, trehalose, mannose, and other carbohydrates as reported for reaction 64.
2.1.8. Lipids synthesis
The lipids for S. cerevisiae are mainly made of phospholipids, sterols and, triacylglycerols which their composition is 3, 54 and 20%, respectively. The main blocks for lipid synthesis consider fatty acids where the palmitic, oleic and linoleic fatty acids are 75% of the total [29]. With regard to the phospholipid composition for S. cerevisiae, it is well known that phosphatidylcholine, phosphatidylethanolamine, and phosphatidylinositol represent more than 90%. For sterols composition in S. cerevisiae, ergosterol represents the most abundant sterol with a porcentage of 90%. Therefore, the lipid synthesis in reaction 67 is presented [29].
2.1.9. Oxidative phosphorylation
For aerobic conditions, the oxidative phosphorylation is included (reactions 32 to 34). Furthermore, in other models that use S. cerevisiae, the stoichiometric oxidative ratio for phosphorylation is P/O = 1.0.
2.1.10. Cellular maintenance
The cellular maintenance in reaction 69 considers the Adenosine Triphosphate (ATP) consumption as for cell-pair functions. The value of S. cerevisiae cellular maintenance is 0.7 mmol ATP/ h g Dry Weight (DW) as it is reported in [29].
2.1.11. Model description
We propose a stoichiometric model that encompasses 117 biochemical reactions and 94 metabolites. From the 117 biochemical reactions, 93 are focused in internal fluxes and 24 to either transport or exchange fluxes between the external environment and the cell. The exchange fluxes analyzed are: glucose, fructose, ammonium, GLN, LEU, ILEU, VAL, THR, PHE, oxygen, sulphate, carbon dioxide, propanol, isobutanol, isoamylic alcohol, amylic alcohol, phenylethanol, ethyl acetate, ethanol, glycerol, acetate, acetaldehyde, succinate and biomass.
2.1.12. Compartmentalization
The compartmentalization was performed by the compounds OAA, ACCOA, NADH and NADPH which participate in reactions in the mitochondria as well as in the cytosol [30].
3.
Results
3.1. Stoichiometric model validation
The proposed stoichiometric model is validated under both conditions, fermentation (anaerobic) and growth (aerobic), which are evaluated with other models reported in the literature [15,30]. For a simulation in silico, the cellular maintenance (m) is considered to be m=0.7 mmol ATP/h g Dry Weight as well as a molecular weight (MW) for the biomass defined as 28 g/C-mol [30].
3.1.1. Stoichimetric model validation under anaerobic conditions
To validate the model under anaerobic conditions, the flux for the consumption of glucose at D = 0.1 h−1 was calculated using the Pirt maintenance model [32]. The glucose flux was 5.917 mmol/g DW h and it was established as an experimental reaction according to Eq (2.5). For oxygen, the flux was zero for anaerobic conditions. Moreover, ammonium was considered as a unique nitrogen source without any restriction (see Table S2). The results showed that large fermentation products (ethanol, biomass and CO2) presented a minimal deviation of 2.41, 0.9 and 2.6%, respectively, while glicerol had a deviation of 21% with respect to the values predicted in silico by Nissen et al. (1997), Table 1. The algorithm estimates only optimal solutions therefore non production of acetate is reported [22].
Table 1.
Theoretical yields obtained with the proposed fermentation model.
3.1.2. Stoichiometric model validation under aerobic conditions
For aerobic conditions, our model was tested against iFF708 which was conformed with 1035 reactions and had been validated experimentally under three conditions: microaerobic fermentation, oxide-fermentative growth, and aerobic growth featuring different oxygen consumptions [15].
1) Microaerobic fermentation
For microaerobic fermentation, the performed restrictions by the simulation were: glucose consumption flux = 14 mmol/g DW h and, oxygen consumption flux = 1 mmol/g DW h. However, ammonium was considered only as a nitrogen source without restrictions for either transport or consumption (see Table S3). The shown analysis of microaerobic fermentation in Table 2 predicted a growth specific velocity of μ=0.322 mmol/g DW h. For this value, deviations of 2.42 and 3.87% with respect to the values in silico and experimental were calculated, respectively. Moreover, the value of ethanol flux = 21.01 mmol/g DW h showed deviations of 1.31 and 4.63% with respect to the in silico and experimental values [15].
Table 2.
Comparison of predictions in silico in microaerobic fermentation.
The oxide-fermentative growth was modelled by considering the following fluxes restrictions: glucose consumption of 12 mmol/g DW h, maximal oxygen consumption of 9 mmol/g DW h and ammonium source with a flux of 1 mmol/g DW h (see Table S4). For this simulation, a growth specific velocity of μ = 0.51 h−1 was obtained which represented a minor value of 3.77% compared with the value in silico, μ = 0.53 h−1, this value is shown in Table 3 and it was similar to the experimental value, both reported by Duarte (2004) [15]. With regard to the ethanol flux (14.03 mmol/g DW h), this represented a variation of 17.11% with respect to the value in silico and 21.10% with respect to the experimental value. Nevertheless, for acetate, our proposed model was not able to predict its production as it also had been reported for larger metabolic models [22]. Furthermore, this behavior could be attributed to the fact that the optimization was focused on maximizing the biomass production.
Table 3.
Comparison of predictions in silico under oxido-fermentative growth.
To analyze this physiological state, a set of restrictions in the fluxes were performed: the flux of glucose consumption was 2.5 mmol/g DW h, the maximal flux consumption for oxygen was 8 mmol/g DW h and ammonium was the nitrogen source without any restriction (see Table S5). Table 4 shows the numerical result of μ = 0.20 h−1 which is 9.09% lower than the value calculated in silico, and it is similar to the experimental value reported by Duarte (2004) [15]. It is observed that our model and the one reported by Duarte (2004) [15] do not predict the production of ethanol and acetate, fluxes equal to zero, this behavior could be attributed to the fact that the objective function is the maximization of growth (μ) in aerobic conditions.
Table 4.
Comparison of predictions in silico in glucose-limited aerobic growth.
In particular, our model showed good results once it was tested with all the above conditions, therefore, it was useful to analyze the effects of dilution rate for the Tequila fermentation process using S. cerevisiae in agave juice.
3.2. Simulation of agave-juice fermentation under continuous culture for tequila production: Dilution rate effect
Fluxes for fructose consumption, ethanol production and the specific growth rate were calculated (see Table S6), with base on the experimental data reported by Moran (2011) [5]. Additionally, the fluxes of volatile compounds were also estimated (see Table S7). These estimated fluxes were specified in our model for performing simulations in silico and for estimating the metabolic fluxes distribution (see Figures 1-3). The optimization process was performed under anaerobic conditions, restricting the oxygen consumption flux in the optimization problem and considering the ammonium phosphate consumption as the unique source of nitrogen. The higher-alcohols production fluxes were specified in the optimization problem as experimental restrictions to obtain a physiological state of the microorganism. The simulation showed that the model predicted an ethanol flux of 5.8 mmol/g DW h (see Table 5), which represented a deviation of -13% with respect to the calculated value from the experimental data of 6.654 mmol/g DW h. The specific growth velocity (μ) was established as the fermentor's continuous dilution rate of μ=D=0.04 h−1. Moreover, the considered cellular maintenance was m=0.7 mmol ATP/g DW h as the value reported for S. cerevisiae by Stephanopolous (1998) [29].
Figure 1.
Central metabolism flux distribution and higher alcohol fermentations of agave juice under continuous culture with Saccharomyces cerevisiae D = 0.04 h−1.
Figure 2.
Central metabolism flux distribution and higher alcohol fermentations of agave juice under continuous culture with Saccharomyces cerevisiae D = 0.08 h−1.
Figure 3.
Central metabolism flux distribution and higher alcohol fermentations of agave juice under continuous culture with Saccharomyces cerevisiae D = 0.12 h−1.
The performed in silico optimizations through FBA were analyzed for the dilution rates of D = 0.04 h−1, D = 0.08 h−1 and D = 0.12 h−1 obtaining a good adjustment regarding ethanol flux (see Table S8). However, the adjustment was reduced as the dilution rate was increased. In addition, for the dilution rate of D=0.16 h−1, and a sugar flux of 12.569 mmol glucose/g DW h, it was not possible to obtain a flux map distribution.
4.
Discussion
Observing the flux distribution map under continuous culture, it is possible to identify the activation of different metabolism paths for Saccharomyces cerevisiae such as glycolysis, pentoses phosphate, fermentative pathways, Krebs cycle, oxidative, phosphorylation, and the consume fluxes for the production of metabolites (see Figures 1-3). Particularly, when the dilution rate increases within the range D = 0.04 to 0.08 h−1 there is an increment in the fluxes of glycerol, acetaldehyde and ammonium as well as in the higher alcohols production (see Figures 1 and 2).
This phenomenon occurs due to a direct relationship between ammonium and the production of higher alcohols. In some fermentation processes, the consumption of ammonium in high levels is correlated with high glycerol, acetate and acetaldehyde syntheses resulting in a reduction of the higher alcohols synthesis [32]. If the dilution rate increases, the microorganism requires higher ammonium concentration but due to the fact that it is limited, the pathway regulation reduces the anabolic pathways of amino acids biosynthesis and then the higher alcohol synthesis increases due to the anabolic pathway from amino acids (LEU, ILEU, VAL, THR and PHE) present in the agave juice [7]. The three physiological states for Saccharomyces cerevisiae are shown in the Figures 1-3 according to the experimental restrictions used in the FBA tool.
Table 6.
Estimation of ethanol flux based on experimental data for different dilution rates.
● a) Map one represents the metabolic flux distribution for (D = 0.04 h−1) (Figure 1) (see Table S6). This physiological state had the highest fermentative capacity of the four states explored by Moran et. al. [5]. Moreover, the highest ethanol concentration was 43.92 g/L and the highest biomass concentration (5.83 g/L) has attained the minimum-residual sugar concentration in the fermentation media (3.94 g/L) by consuming the whole sugar. Additionally, the alcohol productivity reached the second place (1.76 g/L h) in comparison with the other three dilution rates.
● b) Map two represents the metabolic flux distribution for (D = 0.08 h−1) (Figure 2) (see Table S6). The physiological state features one of the two states that had the largest alcohol productivity (2.37 g/L h) similar to the dilution rate D = 0.12 h−1, in comparison with the other two dilution rates. This physiological state had the maximum sugar consumption (5.08 g/L h), the second highest ethanol concentration (29.63 g/L), a residual sugar concentration of 35.34 g/L and a biomass concentration of 3.38 g/L.
● c) Map three represents the metabolic flux distribution for (D = 0.12 h−1) (Figure 3) (see Table S6). It was one of two states that had the highest alcohol productivity (2.37 g/L h), similar to the dilution rate D = 0.08 h−1. Furthermore, it had the second highest value in the sugar consumption rate (4.69 g/L h) and the third highest value in ethanol concentration (19.76 g/L), the second highest value of the residual sugar in the media (59.75 g/L) and also including the third highest value in the biomass concentration (3.04 g/L).
The range of dilution rate (D = 0.08 to 0.12 h−1) showed a low fermentative capacity (biomass concentration) which could be due to a low content of nutrients in agave juice [25] (see Table S6). While the dilution rate increased over the range of D = 0.04 to 0.12 h−1, the concentrations for higher alcohols (n-propanol, iso-butanol, isobutyl alcohol, amyl and isoamyl alcohol and 2-Phenyl-ethanol) also increased.
5.
Conclusions
A metabolic model to analyze Saccharomyces cerevisiae in a continuous fermentation process to obtain agave juice for Tequila production was proposed and implemented. The metabolic model encompassed 94 metabolites and 117 reactions in contrast with far more complex models featuring up to 1035 biochemical reactions such as iFF708. From the 117 reactions, 93 of those were linked to biochemical internal reactions and 24 to transport fluxes between the medium and the microorganism. The developed model was validated under anaerobic and aerobic conditions and it predicted the flux values of the principal metabolites associated with Tequila fermentation. The model allowed us to obtain an estimate of the metabolic flux-map distribution of the strain Saccharomyces cerevisiae and the physiological states of the yeast during the continuous agave-juice fermentation for Tequila production. The estimated fluxes by the model were similar to experimental results reported in literature in which it was possible to visualize the central metabolism and the synthesis pathways of the higher alcohols. The importance of the present genomic model is that it is a simple model which allows the estimation of the metabolic fluxes and determines the cellular physiology state of Saccharomyces cerevisiae.
Acknowledgments
To Consejo Nacional de Ciencia y Tecnología (CONACYT) for the fellowship granted (fellowship number: 174828 and registry number: 165004) and to CIATEJ and UdeG for the facilities and support during the development of the present work.
Conflicts of interest
The authors declare no conflicts of interest.
Appendix
Appendix A. Reaction's list considered in the stoichiometric model.
Glucolysis.
1. GLC + ATP → G6P + ADP
2. FRUC + ATP → F6P + ADP
3. G6P ↔ F6P
4. F6P + ATP → DHAP + G3P + ADP
5. DHAP → G3P
6. G3P + NAD + ADP + Pi → PG3 + ATP + NADH
7. PG3 ↔ PEP
8. PEP + ADP → PYR + ATP
Fermentative pathways.
9. DHAP + NADH → GL3P + NAD
10. GL3P → GL + Pi
11. PYR → ACD + CO2
12. ACD + NADH ↔ ETH + NAD
13. ACD + NADHm ↔ ETH + NADm
14. ACD + NAD → AC + NADH
15. ACD + NADP → AC + NADPH
16. COAm + NADm + PYR → ACCOAm + CO2 + NADHm
17. AC + ATP + COA → ACCOA + ADP + Pi
Pentose phosphate pathway.
18. G6P + 2 NADP → CO2 + 2 NADPH + RL5P
19. RL5P ↔ X5P
20. RL5P ↔ R5P
21. R5P + X5P ↔ G3P + S7P
22. E4P + X5P ↔ F6P + G3P
23. G3P + S7P ↔ E4P + F6P
Citric acid cycle.
24. ATP + OAA → ADP + OAAm + Pi
25. ATP + CO2 + PYR → ADP + OAA + Pi
26. ACCOAm + OAAm → COAm + ICIT
27. ACCOA + OAAm → COA + ICIT
28. ICIT + NADm ↔ AKG + CO2 + NADHm
29. ICIT + NADPm ↔ AKG + CO2 + NADPHm
30. ADP + AKG + NADm + Pi ↔ ATP + CO2 + NADHm + SUCC
B. Alvarez-Samaniego and D. Lannes, Large time existence for 3D water-waves and asymptotics, Invent. Math., 171 (2008), 485-541. doi: 10.1007/s00222-007-0088-4
[3]
B. Alvarez-Samaniego and D. Lannes, A Nash-Moser theorem for singular evolution equations. Application to the Serre and Green-Naghdi equations, Indiana Univ. Math. J., 57 (2008), 97-131. doi: 10.1512/iumj.2008.57.3200
[4]
E. Audusse, F. Bouchut, M.-O. Bristeau, R. Klein and B. Perthame, A fast and stable well-balanced scheme with hydrostatic reconstruction for Shallow Water flows, SIAM J. Sci. Comput., 25 (2004), 2050-2065. doi: 10.1137/S1064827503431090
J.-L. Bona, T.-B. Benjamin and J.-J. Mahony, Model equations for long waves in nonlinear dispersive systems, Philos. Trans. Royal Soc. London Series A, 272 (1972), 47-78. doi: 10.1098/rsta.1972.0032
[7]
P. Bonneton, E. Barthelemy, F. Chazel, R. Cienfuegos, D. Lannes, F. Marche and M. Tissier, Recent advances in Serre-Green Naghdi modelling for wave transformation, breaking and runup processes, European Journal of Mechanics - B/Fluids, 30 (2011), 589-597, URL http://www.sciencedirect.com/science/article/pii/S0997754611000185, Special Issue: Nearshore Hydrodynamics. doi: 10.1016/j.euromechflu.2011.02.005
[8]
F. Bouchut, An introduction to finite volume methods for hyperbolic conservation laws, ESAIM Proc., 15 (2005), 1-17.
[9]
F. Bouchut, Nonlinear Stability of Finite Volume Methods for Hyperbolic Conservation Laws and Well-Balanced Schemes for Sources, Birkhäuser, 2004. doi: 10.1007/b93802
[10]
F. Brezzi, On the existence, uniqueness and approximation of saddle-point problems arising from Lagrangian multipliers, Rev. Française Automat. Informat. Recherche Opérationnelle Sér. Rouge, 8 (1974), 129-151.
[11]
M.-O. Bristeau and B. Coussin, Boundary Conditions for the Shallow Water Equations Solved by Kinetic Schemes, Rapport de recherche RR-4282, INRIA, 2001, URL http://hal.inria.fr/inria-00072305, Projet M3N.
[12]
M.-O. Bristeau, N. Goutal and J. Sainte-Marie, Numerical simulations of a non-hydrostatic Shallow Water model, Computers & Fluids, 47 (2011), 51-64. doi: 10.1016/j.compfluid.2011.02.013
M.-O. Bristeau and J. Sainte-Marie, Derivation of a non-hydrostatic shallow water model; Comparison with Saint-Venant and Boussinesq systems, Discrete Contin. Dyn. Syst. Ser. B, 10 (2008), 733-759. doi: 10.3934/dcdsb.2008.10.733
[15]
R. Camassa, D. Holm and J. Hyman, A new integrable shallow water equation, Adv. Appl. Math., 31 (1994), 1-33. doi: 10.1016/S0065-2156(08)70254-0
[16]
F. Chazel, D. Lannes and F. Marche, Numerical simulation of strongly nonlinear and dispersive waves using a Green-Naghdi model, J. Sci. Comput., 48 (2011), 105-116. doi: 10.1007/s10915-010-9395-9
[17]
A. J. Chorin, Numerical solution of the Navier-Stokes equations, Math. Comp., 22 (1968), 745-762. doi: 10.1090/S0025-5718-1968-0242392-2
[18]
M.-W. Dingemans, Wave Propagation Over Uneven Bottoms, Advanced Series on Ocean Engineering - World Scientific, 1997. doi: 10.1142/9789812796042
[19]
M. Dingemans, Comparison of Computations with Boussinesq-like Models and Laboratory Measurements, Technical Report H1684-12, AST G8M Coastal Morphodybamics Research Programme, 1994.
E. Godlewski and P.-A. Raviart, Numerical Approximations of Hyperbolic Systems of Conservation Laws, Applied Mathematical Sciences, vol. 118, Springer, New York, 1996. doi: 10.1007/978-1-4612-0713-9
[24]
A. Green and P. Naghdi, A derivation of equations for wave propagation in water of variable depth, J. Fluid Mech., 78 (1976), 237-246. doi: 10.1017/S0022112076002425
[25]
First Int. Conf. on Comput. Methods in Flow Analysis, Okayama University, Japan.
[26]
J.-L. Guermond, Some implementations of projection methods for Navier-Stokes equations, ESAIM: Mathematical Modelling and Numerical Analysis, 30 (1996), 637-667, URL http://eudml.org/doc/193818.
D. Lannes and P. Bonneton, Derivation of asymptotic two-dimensional time-dependent equations for surface water wave propagation, Physics of Fluids, 21 (2009), 016601. doi: 10.1063/1.3053183
[30]
O. Le Métayer, S. Gavrilyuk and S. Hank, A numerical scheme for the Green-Naghdi model, J. Comput. Phys., 229 (2010), 2034-2045. doi: 10.1016/j.jcp.2009.11.021
[31]
R.-J. LeVeque, Finite Volume Methods for Hyperbolic Problems, Cambridge University Press, 2002. doi: 10.1017/CBO9780511791253
[32]
O. Nwogu, Alternative form of Boussinesq equations for nearshore wave propagation, Journal of Waterway, Port, Coastal and Ocean Engineering, ASCE, 119 (1993), 618-638. doi: 10.1061/(ASCE)0733-950X(1993)119:6(618)
[33]
D. Peregrine, Long waves on a beach, J. Fluid Mech., 27 (1967), 815-827. doi: 10.1017/S0022112067002605
[34]
O. Pironneau, Méthodes Des Éléments Finis Pour Les Fluides., Masson, 1988.
[35]
R. Rannacher, On Chorin's projection method for the incompressible Navier-Stokes equations, in The Navier-Stokes Equations II - Theory and Numerical Methods (eds. G. Heywood John, K. Masuda, R. Rautmann and A. Solonnikov Vsevolod), vol. 1530 of Lecture Notes in Mathematics, Springer Berlin Heidelberg, 1992, 167-183, URL http://dx.doi.org/10.1007/BFb0090341. doi: 10.1007/BFb0090341
[36]
11th AIAA Computational Fluid Dynamic Conference, Orlando, FL, USA.
[37]
J. Shen, On error estimates of the penalty method for unsteady Navier-Stokes equations, SIAM J. Numer. Anal., 32 (1995), 386-403. doi: 10.1137/0732016
[38]
W. C. Thacker, Some exact solutions to the non-linear shallow water wave equations, J. Fluid Mech., 107 (1981), 499-508. doi: 10.1017/S0022112081001882
Assyr Abdulle, Doghonay Arjmand, Edoardo Paganoni,
Exponential decay of the resonance error in numerical homogenization via parabolic and elliptic cell problems,
2019,
357,
1631073X,
545,
10.1016/j.crma.2019.05.011
3.
François Bignonnet, Karam Sab, Luc Dormieux, Sébastien Brisard, Antoine Bisson,
Macroscopically consistent non-local modeling of heterogeneous media,
2014,
278,
00457825,
218,
10.1016/j.cma.2014.05.014
4.
Doghonay Arjmand, Olof Runborg,
A time dependent approach for removing the cell boundary error in elliptic homogenization problems,
2016,
314,
00219991,
206,
10.1016/j.jcp.2016.03.009
5.
Eric Cancès, Virginie Ehrlacher, Frédéric Legoll, Benjamin Stamm, Shuyang Xiang,
An Embedded Corrector Problem for Homogenization. Part I: Theory,
2020,
18,
1540-3459,
1179,
10.1137/18M120035X
J.-C. Mourrat,
Efficient Methods for the Estimation of Homogenized Coefficients,
2019,
19,
1615-3375,
435,
10.1007/s10208-018-9389-9
8.
Ivo Babuška, Mohammad Motamed,
A fuzzy-stochastic multiscale model for fiber composites,
2016,
302,
00457825,
109,
10.1016/j.cma.2015.12.016
9.
Antti Hannukainen, Jean-Christophe Mourrat, Harmen T. Stoppels,
Computing homogenized coefficientsviamultiscale representation and hierarchical hybrid grids,
2021,
55,
0764-583X,
S149,
10.1051/m2an/2020024
10.
Julian Fischer,
The Choice of Representative Volumes in the Approximation of Effective Properties of Random Materials,
2019,
234,
0003-9527,
635,
10.1007/s00205-019-01400-w
11.
Assyr Abdulle, Doghonay Arjmand, Edoardo Paganoni,
A parabolic local problem with exponential decay of the resonance error for numerical homogenization,
2021,
31,
0218-2025,
2733,
10.1142/S0218202521500603
12.
Matti Schneider, Marc Josien, Felix Otto,
Representative volume elements for matrix-inclusion composites - a computational study on the effects of an improper treatment of particles intersecting the boundary and the benefits of periodizing the ensemble,
2022,
158,
00225096,
104652,
10.1016/j.jmps.2021.104652
13.
Antoine Gloria, Zakaria Habibi,
Reduction in the Resonance Error in Numerical Homogenization II: Correctors and Extrapolation,
2016,
16,
1615-3375,
217,
10.1007/s10208-015-9246-z
14.
Ronan Costaouec, Claude Le Bris, Frédéric Legoll,
Variance reduction in stochastic homogenization: proof of concept, using antithetic variables,
2010,
50,
1575-9822,
9,
10.1007/BF03322539
15.
ANTOINE GLORIA,
REDUCTION OF THE RESONANCE ERROR — PART 1: APPROXIMATION OF HOMOGENIZED COEFFICIENTS,
2011,
21,
0218-2025,
1601,
10.1142/S0218202511005507
16.
Assyr Abdulle, Doghonay Arjmand, Edoardo Paganoni,
An Elliptic Local Problem with Exponential Decay of the Resonance Error for Numerical Homogenization,
2023,
21,
1540-3459,
513,
10.1137/21M1452123
Sean P. Carney, Milica Dussinger, Björn Engquist,
On the Nature of the Boundary Resonance Error in Numerical Homogenization and Its Reduction,
2024,
22,
1540-3459,
811,
10.1137/23M1594492
19.
Nicolas Clozeau, Marc Josien, Felix Otto, Qiang Xu,
Bias in the Representative Volume Element method: Periodize the Ensemble Instead of Its Realizations,
2024,
24,
1615-3375,
1305,
10.1007/s10208-023-09613-y
Nora Aïssiouene, Marie-Odile Bristeau, Edwige Godlewski, Jacques Sainte-Marie. A combined finite volume - finite element scheme for a dispersive shallow water system[J]. Networks and Heterogeneous Media, 2016, 11(1): 1-27. doi: 10.3934/nhm.2016.11.1
Nora Aïssiouene, Marie-Odile Bristeau, Edwige Godlewski, Jacques Sainte-Marie. A combined finite volume - finite element scheme for a dispersive shallow water system[J]. Networks and Heterogeneous Media, 2016, 11(1): 1-27. doi: 10.3934/nhm.2016.11.1