Chem. Senses 29: 235-245,
2004
© Oxford University Press 2004
Modeling the Kinetics of Flavour Release during Drinking
Firmenich S.A., 7 rue de la Bergère, 1217 Meyrin 2 Geneva, Switzerland
Correspondence to be sent to: Dr Valéry Normand, Firmenich S.A., 7 rue de la Bergère, 1217 Meyrin 2 Geneva, Switzerland. e-mail: valery.normand{at}firmenich.com
| Abstract |
|---|
|
|
|---|
Novel mathematical models for flavour release during drinking are described, based on the physiology of breathing and swallowing. Surprisingly, we conclude that most flavour molecules arriving in the nose are extracted from liquid left in the throat, after swallowing. The models are fit to real time flavour release data obtained using APCImass spectrometry. Before modelling, raw data are corrected for the effects of varying airflow rate, using the signal from acetone in exhaled air. A simple equilibrium batch extraction model correctly describes flavour release during the first breaths after swallowing a flavoured liquid. It shows that for eight volatiles, whose in vitro airwater partition coefficients vary by a factor of 500, the apparent in vivo airsaliva partition coefficients vary only by a factor of five. To interpret the kinetics of flavour release longer after swallowing, diffusion of flavour into the throat lining is included. This is done using a three-layer model for mass transfer in the throat. An analytical solution of this model gives good fits to typical data. These models de-couple the physiological and physico-chemical aspects of flavour release, clarifying the effect of behaviour on in-vivo flavour release.
Key words: : flavour release, mass-transfer, physiology
| Introduction |
|---|
|
|
|---|
The goal of cognitive science is to understand how perception is related to sensory stimuli. For the olfactive component of flavour (the aroma of food during eating) this approach was impossible until recently because the stimulus, the concentration of flavour molecules on the nasal epithelium, was unknown. The situation has changed now that the concentration of flavour molecules in air exhaled through the nose can be measured in real time, using mass spectrometry (a recent review is Taylor and Linforth, 2000
The aim of this study is to better understand the physiological and physico-chemical origins of the variation in the kinetics of flavour release. We are particularly interested in the influence of food composition (fat content, for instance). Improved modelling should improve our understanding of how flavour composition, anatomy and physiology influence flavour release.
Existing mathematical models for in vivo (or in mouth) flavour release are adaptations of models for in vitro release (for a review, see Taylor, 2002
). However until now they have not correctly incorporated the physiology of eating and drinking. Invariably, flavour-rich air in the mouth is assumed to pass directly into the nose (e.g. Harrison, 1998
, 2000). This assumption is almost always incorrect, as there is a gas-tight seal between mouth and nose for most of the time (Buettner et al., 2001
). During eating this seal is not perfect, so flavour can reach the nose intermittently during chewing. During drinking, which we study here, the seal is perfect, except when swallowing.
The paper is structured as follows. First we describe the physiology of swallowing. Then we discuss the raw data from flavour release and describe some essential improvements to the data treatment. Finally, we describe two different mathematical models for flavour release, comparing each with experimental data to show its strengths and weaknesses.
| Physiology of swallowing |
|---|
|
|
|---|
The throat is used for both breathing and feeding, so each swallow involves a risk. This is especially true during drinking, as even a small opening at the back of the mouth will allow liquid food to flow into the lungs, with a risk of infection or choking. Recent studies of the physiology of swallowing (Dodds et al., 1990a
|
At the upper end this tube separates to form two cavities: the mouth for feeding and the nasal cavity for breathing. At the lower end the tube also divides into two: the oesophagus, for feeding, and the trachea, for breathing. The mechanisms represented in Figure 1b,c are well described elsewhere (Bosma, 1980
| Materials and methods |
|---|
|
|
|---|
Flavour release in expired air was measured using atmospheric pressure chemical ionizationmass spectrometry (APCIMS; Micromass, Manchester, UK) (Taylor et al., 2000
Panellists were asked to swallow according to the following protocols.
- Protocol 1: freely.
- Protocol 2: (i) hold mouth shut without swallowing for 2 min; (ii) inject 8 ml of flavoured solution into closed mouth; (iii) hold solution in mouth for 30s; (iv) swallow all the solution; and (v) breath normally through the nose, keeping mouth shut. The subjects were given no other instructions on how to breathe.
- Protocol 3: identical to protocol 2 except that the subjects were asked to breathe at a given frequency.
Volatiles in the exhaled breath were recorded using the nose space technique for several minutes. The flow rate of inhaled and exhaled air was recorded in real time using a physiological air flow meter (World Precision Instruments, Stevenage, UK). It was placed in-line with the APCIMS interface, so that the flow rate and composition of expired air were measured simultaneously.
| Results |
|---|
|
|
|---|
Effect of behaviour on the kinetics of flavour release
Real time measurements show that the kinetics of flavour release is extremely subject-dependent. This point is illustrated in Figure 2, which shows the flavour concentration in exhaled air for three subjects drinking the same menthol solution with no artificial restrictions on their behaviour (protocol 1).
|
These patterns of flavour release show huge differences, due to the different ways in which people breathe and swallow. It is very likely that the perceptions resulting from these different stimuli will also be different. Electromyography measurements are needed to directly determine the panellists behaviour whilst measuring flavour release (Hodgson et al., 2003
In fact, Land speculated 10 years ago (Land, 1994
) that the only source of aroma compounds in the nose during eating was from the small swallow breath immediately after swallowing. Buettner et al. (2001
) agreed with this point of view, on the basis of their study using videofluoroscopy and real-time MRI. Figure 2 shows that this is not the case. Significant amounts of aroma compound are exhaled from the nose in every breath, even 1 min after swallowing.
Improved data treatment
A major problem in interpreting raw nose-space data is that the signal is proportional to the amount of flavour per unit time, which is the product of the air flow rate and the flavour concentration. Breathing faster will have exactly the same effect on the signal as increasing the flavour concentration, so an unusually large peak could be due to either higher flavour concentration or faster exhalation. Unambiguous interpretation of the raw data is impossible. The simplest way to remove this ambiguity is to measure the air flow rate in the nose. However, physiological air flow meters are expensive and integration of the flow rate and flavour release data is tedious. A more practical solution would be to use an internal standard. The internal standard should be easy to measure in the exhaled air and its concentration must not change during an experiment. We find that acetone is suitable. It is produced by fat metabolism in the liver and then passes into the blood. It partitions from the blood stream into the lungs and so is present in exhaled air. It is easily measured by the nose-space technique.
We made a series of measurements of acetone peak area and exhaled volume per breath. The measurements used one subject and lasted >1 h. Plotting peak area against breath volume gave a good straight line (r2 = 0.978) passing through the origin (data not shown). This result shows that over the time scale of our experiments, the acetone concentration in exhaled air is sufficiently constant. Therefore, it was measured simultaneously with the flavour molecules and used to allow for changing air flow rate in the following way.
The average flavour concentration in the ith exhalation is:

Here V is the volume of expired air and j is the ratio between the amount of flavour molecules and the signal intensity (I). The integral is over one breath peak, so rewriting in terms of the peak area, Si, gives:

Using this expression for acetone and rearranging gives the volume exhaled:
![]()
Substituting (2) into (1b) gives:

J is constant for each flavour molecule, subject and experiment. Equation (3) shows that the average concentration of flavour in each breath peak is directly proportional to the ratio between the areas of the flavour and acetone peaks. However, the absolute flavour concentration cannot be determined, because the acetone concentration is unknown. Nevertheless, the shape of the kinetics of flavour release can be compared. To eliminate the effect of air flow rate, all the peak areas are re-scaled relative to the area of the largest acetone peak. Figure 3 shows a typical result of this transformation for menthol (Figure 3a,c) and acetone (Figure 3b,d).
|
Figure 3a shows that the subject swallowed twice during the experiment (peaks 2 and 11). Figure 3c shows that the raw data for the peak areas (squares) do not decrease regularly. However, the corrected peak areas (circles) are much more regular. This regularity is hidden in the raw data. For the rest of this article, flavour peak areas for each experiment are normalised relative to the largest acetone peak in the experiment.
Before describing the modelling, we want to re-emphasize the importance of using this procedure. In the raw data, peak area for each breath is proportional to the product of exhaled volume and flavour concentration. Data in this form cannot be interpreted in terms of flavour concentration. After normalization, the peak area is proportional to the average concentration (or amount) of flavour released. Using the internal standard removes the ambiguity in interpretation of the data. For instance, in Figure 3c, the raw data after each swallow show a rise in the peak area (sixth breath after the first swallow and fourth breath after the second). The corrected data show that these increases are artefacts, caused by variation in breath volume.
A large set of experimental data (following protocol 2) was obtained from nine subjects drinking aqueous solutions of eight flavour molecules with a wide range of volatilities (see Table 1). The raw data were first treated as described above to average over breath peaks and correct for the effects of air flow rate. Figure 4 shows a typical result, with the area of the breath peaks on a logarithmic scale.
|
|
In general, we find that after the first breath, the points fall onto two straight lines, with the first steeper than the second. The first point always falls above the initial steeper line. Following Land (1994
Model I: equilibrium batch extraction
Surprisingly, the simplest possible model for mass transfer in the throat gives a good description of our data. We assume that the air from each breath has time to equilibrate with the liquid remaining in the throat after swallowing. In other words, we assume that flavour release is due to equilibrium batch extraction. Time is absent, as the system is assumed to be at equilibrium. The same model has been used in vitro to measure the partition coefficients of volatile organic compounds (Mackay et al., 1979
; Nielsen et al., 1994
; Chaintreau et al., 1995
). We explain this method in detail, as it is applied in exactly the same way to our in vivo data. A syringe is partially filled with a known volume, VWater, of flavoured aqueous solution containing NWater moles of flavour compound. The rest of the syringe is filled with a volume VAir of clean air. Once airwater equilibrium has been reached, the air is pushed out of the syringe and the number of flavour molecules in it, NAir, is determined. The syringe is then refilled with clean air and the whole process is repeated n times. The amount of flavour in the nth batch of expelled air can be calculated using:

KAW is the airwater partition coefficient of the flavour molecule. These equations can be used to calculate NAir as a function of the initial number of flavour molecules, (NWater)0, using:

Taking logarithms this becomes:

This equation has the form y = px + q, with p and q constants, so the logarithm of the number of molecules present in a batch of air is proportional to the batch number. Plotting the left-hand side of this equation against the batch number gives a straight line whose slope, p, depends on the partition coefficient, KAW, but is independent of the initial amount of flavour. Rearranging the expression for the slope, we obtain:

The key advantage of the method is that partition coefficients can be determined without absolute concentration measurements. Only the ratio of the volume of liquid to the volume of air is needed to calculate the partition coefficient from the slope. By analogy, flavour release data from in vivo experiments can be used to obtain an apparent partition coefficient, without needing the absolute flavour concentration in the breath.
We assume that the volume ratio of saliva in the throat to air in the upper respiratory tract is 0.008 (i.e. 100 ml of air for 0.8 ml of saliva in the throat). We justify this value later. This ratio is used in equation (7) to calculate an apparent in vivo airsaliva partition coefficient (KAS) from the experimental data. This method of analysis was applied to all the data measured above. Table 1 summarizes the results and compares them with literature values for the volatility (KAW), interpolated to 25°C (J.-Y. de Saint Laumer, private communication). The molecules are listed in order of increasing volatility.
Note that the coefficients of variation tend to be larger for high volatilities. These data are plotted in Figure 5. A solid line indicates the result if the in vivo and in vitro volatilities were equal.
|
The most striking observation is that the apparent in vivo volatilities all fall within an extremely narrow range: they only vary by a factor of five, whereas the in vitro volatilities vary by a factor of >500. Choosing a different air/liquid volume ratio would shift the in vivo values vertically, but would not change the slope. Despite the greatly reduced range, KAS is still roughly proportional to KAW. In fact, it is almost exactly proportional to (KAW)1/5. This formula is extremely useful for predicting the apparent in vivo volatilities of flavour molecules.
When mass transfer between saliva and air is fast, the flavour molecule will be in apparent equilibrium between the two phases and the slope p (and hence KAS) will not depend on the breathing frequency. In this case model I will be applicable. On the other hand, if mass transfer is slow, the system will not be in apparent equilibrium and the slope (and hence KAS) will depend on the breathing frequency.
To test the effect of breathing frequency, two subjects were instructed to breathe with fixed frequency (protocol 3). Experimentally we find that the assumption of equilibrium in vivo is reasonable for menthol, ethanol, octanone, ethyl acetate and citral. For anethole and iso-amyl acetate the apparent partition coefficient depends on the breathing frequency, so not even apparent equilibrium is reached. Figure 6 shows the experimentally measured effects of breathing frequency for anethole, which shows breathing frequency dependence, and menthol, which does not.
|
For breathing frequencies of 2 and 10/min, menthol gives comparable results, whereas the release of anethole depends strongly on breathing rate.
This first model describes the initial slope that begins with the second breath after swallowing. Also, it is very useful for determining the apparent in vivo partition coefficient of a flavour molecule, which controls its release after swallowing. However, Figure 4 showed that this simple model only applies for a few breaths after swallowing a liquid. To explain the slower flavour release that occurs longer after swallowing a more elaborate model is needed.
Model II: non-equilibrium and absorption
Model I only applies to the first few breaths after swallowing. Figure 4 showed that after this, the slope decreases continuously and then reaches another regime with a different, constant slope. This can only mean that, after a few breaths, flavour molecules lining the throat move into an environment where they are released more slowly. This environment can only be the throat lining, or mucosa as suggested elsewhere (Buettner et al., 2002
). If the flavour molecules were insoluble in the mucosa, then model I would always apply and the logarithm of the flavour peak area would be linear with the breath number throughout the experiment.
Model II includes a layer of mucosa underneath the layer of saliva, so we now have three layers: air, saliva and mucosa. However, just adding another layer to model I is insufficient for two slopes to be observed. Two conditions must be fulfilled: (i) the affinity of flavour molecules must be higher for the mucosa than for air, that is KSM > KAS and (ii) the rate of mass transfer from saliva to air must be much faster than from saliva to mucosa. These two conditions are necessary and sufficient to observe two distinct slopes, with the initial slope being steeper.
Under these conditions, at short times flavour will be released rapidly from the saliva layer into the air. It will also be transferred slowly into the mucosa, but at this point, the transfer is so slow that its effect is not apparent. As time goes on, the loss of flavour to the mucosa becomes large enough to affect the transfer from saliva to the air, which starts to slow, so the slope decreases. Next, the concentration in the mucosa increases to a point where it is in equilibrium with the decreasing concentration in the saliva. This is its maximum concentration. After the next breath, the concentration in the saliva falls below that required to transfer from saliva to mucosa and the process reverses. Finally, the second linear regime is reached where flavour release is limited by transfer out of the mucosa.
Model II abandons the assumption of equilibriumtime is included. Mass transfer between saliva/air and saliva/mucosa are treated as continuous in time and coupled.
Following Hills and Harrison (1995
), we choose to model mass transfer using the two-film model. This assumes that both phases are perfectly mixed, so flavour concentrations are equal throughout each phase, except at a very thin region next to each side of the interface. The theoretical expression for the rate of concentration flow through the interface per unit area per unit time (the flux) is given by:
![]()
Here, k is the overall mass transfer coefficient. This equation shows that mass transfer is driven by the difference between: (i) the actual concentration in one phase and (ii) the concentration that it would have at equilibrium with the actual concentration in the other phase. Note that this expression uses concentration differences in the air. A slightly different expression, but the same flux, would be found if the concentrations in saliva were chosen (see, for instance, Cussler, 1997
). To emphasize this point, k should really be called the overall air-side mass transfer coefficient.
We now write the mass balance for each environment, leading to a system of three differential equations:

This system of first-order differential equations can be solved analytically. The solution of each is the sum of three exponentials:

where
i,
i,
i and µi depend on: (i) the mass transfer coefficients (kSA and kSM), (ii) the phase volumes (VA, VS and VM), (iii) the interfacial areas (SSA and SSM) and (iv) the partition coefficients (KAS and KMS). Qi also depends on the initial quantities [nA(0), nS(0) and nM(0)]. Of these 12 parameters, only the mass transfer coefficients and partition coefficients are not fixed by the anatomy, which decreases the number of free parameters to four. In addition, model I gives a good approximation to the apparent KAS, so only three parameters are free.
We now make some reasonable approximations for the variables in the model. The total volume of the upper respiratory tract of a 75 kg human is
100 cm3 (Ménache et al., 1997
). Then 8 cm3 of flavoured solution was put into the mouth. We assume that 10% of the solution mixed with saliva covers the throat lining after swallowing. Therefore, 0.8 cm3 covers the inner surface of a cylinder delimited by the soft palate at the top and the epiglottis at the bottom. The length of the cylinder is estimated as 3.5 cm. The diameter of the throat is taken as 3 cm. These assumptions give the thickness of the saliva coating the throat as 0.25 mm and the contact surfaces for mucosa/saliva and saliva/air as 33 cm2. At this point, the only remaining adjustable parameters are the two mass transfer coefficients and the apparent saliva-mucosa partition coefficient. Following the careful study of release from water by Marin et al. (2000
), we set the mass transfer coefficient for airsaliva (kSA) to 3 cm/s. Finally, the presence of two slopes, with the first steeper than the second, implies that flavour in saliva is preferentially partitioned into the mucosa rather than air, i.e. KMS > KAS. If this were not true, the second slope would not exist.
This model is continuous in time. However, it retains the idea of batch extraction, as all the air in each breath is assumed to arrive instantly in the throat at the beginning of each exhalation and to disappear instantly at the end of each exhalation. Time and the concentration of flavour molecules in the air are both reset to zero at the start of each exhalation, when clean air arrives in the throat. At the start of each breath, the numbers of molecules in saliva and mucosa are assumed to be equal to those at the end of the previous breath.
We can now calculate the concentrations of flavour in air, saliva and mucosa. Figure 7 shows some typical results.
|
At time zero, all the flavour is in the saliva. At short times it is transferred to both mucosa and air. At this stage, the amount in the mucosa increases, because it can accumulate from breath to breath. After
1 min, in this case, the slope for the mucosa changes sign and those for saliva and air become parallel. This is the point at which the flavour concentration in the saliva falls below the level required for transfer from saliva to mucosa. Transfer changes direction to go from the mucosa to the saliva. Transfer to the air is now limited by the rate of flavour transfer from the mucosa, so flavour release falls to a new lower rate. A non-linear least-squares method was used to fit model II to experimental data. Model I showed that the overall shape of a flavour molecules release depends on its volatility, but not on its initial concentration. However, the amplitude of the signal is proportional to the initial concentration, so an arbitrary vertical shift in log space is needed to adjust the model to the experimental data. Figure 8 shows the good fit of model II to typical experimental data.
|
The best fit is obtained with the following parameter values: kSA = 0.3 cm/s, kSM = 1.7 103 cm/s and KMS = 192. As implied by the existence of two slopes in Figure 4, the mass transfer coefficient for salivaair is
200 times larger than that for saliva-mucosa. In addition, the large value of KMS shows that flavour partitioning is in favour of the mucosa. Effects of varying parameter values in model II
Figure 9 shows how varying each of the three key parameters and the breathing frequency affects the kinetics of flavour release predicted by model II. Arrows show how the shape of the curve changes as the parameter increases.
|
Figure 9a shows that small changes in the salivaair mass transfer coefficient have a large effect on both the initial and final slopes. The initial slope decreases as this mass transfer coefficient increases, whereas the final slope increases. However, both slopes tend towards limiting values as the rate of transfer increases.
In contrast, Figure 9b shows that changing the saliva-mucosa partition coefficient has a smaller effect on the kinetics of flavour release. The final slope is more affected than the initial slope.
Figure 9c shows that varying the salivamucosa mass transfer coefficient has a significant effect on the point of inflection separating the two slopes. The lower the salivamucosa mass transfer coefficient, the later the point of inflexion occurs. Intuitively, slowing the rate of transfer increases the time needed to reach saliva-mucosa equilibrium.
Figure 9d shows the effects of changing the breathing frequency. Under our chosen conditions, model II agrees with the observation illustrated in Figure 6b: changes in breathing frequency have no effect on the initial slope under normal human breathing frequencies. Only unrealistically low breathing frequencies influence the initial flavour release. At normal breathing frequencies, the flavour does not have time to penetrate sufficiently into the mucosa.
| Discussion |
|---|
|
|
|---|
The models formulated here for flavour release during drinking combine realistic physiological assumptions with the simplest possible models for the liquid/air mass transfer. This is a deliberate choice: we want to gain insight from approximate models with a few parameters before progressing to more detailed models. The realistic physiology distinguishes our work from that of Harrison (1998
The models for mass transfer are clearly unrealistic: model I assumes equilibrium. Both models assume instant, perfect mixing in all phases. Mass transfer is lumped over each exhalation. They are far from the sophistication of Keyhani et al.s (1997) model for odorant mass transfer in the nose, during sniffing. The simplicity of our models means that, for instance, we have nothing to say about the change in flavour concentration during a breath. Nor can we address the question of how changes in the concentration of aroma molecules over the nasal epithelium might affect perception, as discussed by Keyhani et al. (1997
). Nevertheless, we think that they give insight into the basic processes by which aroma molecules are transferred to exhaled air during drinking. The models make non-trivial, testable predictions that we have shown are compatible with a sizeable data set.
Unexpectedly, model I shows that apparent airsaliva partition coefficients fall within a very narrow range. The decreasing slope at long times, described by model II, can only be interpreted by assuming that the airmucus partition coefficient is significantly lower than that for airsaliva. Keyhani et al. (1997
) assume that airmucus partition coefficients are equal to those for airwater, stating that there are no experimental measurements for airmucus partition. Our results are incompatible with this assumption.
| Summary and conclusions |
|---|
|
|
|---|
The kinetics of flavour release during drinking vary hugely between subjects. These differences are mainly due to variations in the pattern of swallowing and breathing. Interpretation of nose space data in terms of flavour concentrations can only be done after correcting for the effect of air flow rate. Exhaled acetone can be used to do this. We have shown that during drinking, only a small amount of aroma-rich air passes direct from mouth to nose, just after swallowing. The rest of the time, the aroma signal in the nose is extracted from liquid in the throat by exhaled air. Flavour release after a swallow during drinking falls into three regimes: (i) the first breath, which is abnormally large, contains air from the mouth; (ii) the next few breaths where release can be assumed to be from a liquid film coating the throat; and (iii) the rest of the breaths for which interaction with the mucosa must be considered.
Two simple models for mass transfer in the throat correctly encapsulate the trends observed experimentally when the data are treated breath by breath. More sophisticated models will be needed to understand the fine structure of flavour release during each breath.
Our overall aim was to obtain physico-chemical parameters from real time flavour release data that are independent of the panellists behaviour. We feel that we have succeeded in this task, perhaps surprisingly well, considering the simplicity of the models.
| Acknowledgements |
|---|
Our thanks to Tony Blake, Bénedicte Ribot and Jérôme Barra for their help and advice.
| References |
|---|
|
|
|---|
Bosma, J.F. (1980) Physiology of the mouth, pharynx, and oesophagus. Physiology of the mouth. In Paparella, M.M. and Shumrick, D.A. (eds), Otolaryngology. W.B. Saunders, Philadelphia, PA, pp. 319331.
Buettner, A., Beer, A., Hannig, C. and Settles, M. (2001) Observation of the swallowing process by application of videofluoroscopy and real time magnetic resonance imagingconsequences for retronasal aroma stimulation. Chem. Senses, 26, 12111219.
Buettner, A., Beer, A., Hannig, C., Settles, M. and Schieberle, P. (2002) Physiological and analytical studies on flavor perception dynamics as induced by the eating and swallowing process. Food Qual. Pref., 13, 497504.[CrossRef]
Cussler, E.L. (1997) Diffusion Mass Transfer in Fluid Systems. Cambridge University Press, Cambridge.
Chaintreau, A., Grade, A. and Munoz-Box, R. (1995) Determination of partition coefficients and quantification of headspace volatile compounds. Anal. Chem., 67, 33003304.[CrossRef]
Dodds, W.J., Stewart, E.T. and Logemann, J.A. (1990a) Physiology and radiology of the normal oral and pharyngeal phases of swallowing. Am. J. Roentgenol., 154, 953963.
Dodds, W.J., Logemann, J.A. and Stewart, E.T. (1990b) Radiologic assessment of abnormal oral and pharyngeal phases of swallowing. Am. J. Roentgenol., 154, 965974.
Firmin, H., Reilly, S. and Fourcin, A. (1997) Non-invasive monitoring of reflexive swallowing. Speech, Hearing Lang., 10, 171184.
Harrison, M. (1998) Effect of breathing and saliva flow on flavor release from liquid foods. J. Agric. Food Chem., 46, 27272735.[CrossRef]
Harrison, M. (2000) Mathematical Models of Release and Transport of Flavors from Foods in the Mouth to the Olfactory Epithelium, ACS Symposium Series 763 (Flavour Release), pp. 179191.
Hills, B.P. and Harrison, M. (1995) Two-film theory of flavour release from solids. Int. J. Food Sci. Technol., 30, 425436.
Hodgson, M., Linforth, R.S.T. and Taylor, A.J. (2003) Simultaneous real time measurements of mastication, swallowing, nasal airflow and aroma release. J. Agric. Food Chem., 51, 50525057.[CrossRef][Web of Science][Medline]
Keyhani, K., Scherer, P.W. and Mozell, M.W. (1997) A numerical model of nasal odorant transport for the analysis of human olfaction. J. Theor. Biol., 186, 279301.[CrossRef][Web of Science][Medline]
Land, D.G. (1994) Perspectives on the Effects of Interactions on Flavor Perception: An Overview. ACS Symposium Series 633, pp. 211.
Linforth, R.S.T. and Taylor, A.J. (1998) Apparatus and methods for the analysis of trace constituents of gases. European patent EP 0819 937 A2.
Linforth, R.S.T. and Taylor, A.J. (2000) Persistence of volatile compounds in the breath after their consumption in aqueous solutions. J. Agric. Food Chem., 48, 54195423.[CrossRef][Web of Science][Medline]
Mackay, D., Shiu, W.Y. and Sutherland, R.P. (1979) Determination of airwater Henrys law constants for hydrophobic pollutants. Environ. Sci. Technol., 13, 333337.[CrossRef]
Marin, M., Baek, I. and Taylor, A.J. (2000) Flavor Release as a Unit Operation: A Mass Transfer Approach Based on a Dynamic Headspace Dilution Method. ACS Symposium Series 763, pp. 153165.
Ménache, M.G., Hanna, L.M., Gross, E.A., Lou, S.-R., Zinreich, S.J., Leopold, D.A., Jarabek, A.M. and Miller, F.J. (1997) Upper respiratory tract surface areas and volumes of laboratory animals and humans: considerations for dosimetry models. J. Toxicol. Environ. Health, 50, 475506.[CrossRef][Web of Science][Medline]
Nielsen, F., Olsen, E. and Fredenslund, A. (1994) Henrys law constants and infinite dilution activity coefficients for volatile organic compounds in water by a validated batch air stripping method. Environ. Sci. Technol., 28, 21332138.
Taylor, A.J. (2002) Release and transport of flavours in vivo: physicochemical, physiological, and perceptual considerations. Compr. Rev. Food Sci. Food Safety, 1, 4557.
Taylor, A.J. and Linforth, R.S.T. (2000) Techniques for Measuring Volatile Release In-vivo during Consumption of Food. ACS Symposium Series 763, pp. 821.
Taylor, A.J., Linforth, R.S.T., Harvey, B.A. and Blake, A. (2000) Atmospheric pressure chemical ionisation mass spectrometry for in vivo analysis of volatile flavour release. Food Chem., 71, 327338.[CrossRef]
Accepted January 28, 2004
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
R. M.A.J. Ruijschop, M. J.M. Burgering, M. A. Jacobs, and A. E.M. Boelrijk Retro-Nasal Aroma Release Depends on Both Subject and Product Differences: A Link to Food Intake Regulation? Chem Senses, June 1, 2009; 34(5): 395 - 403. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. C. Trelea, S. Atlan, I. Deleris, A. Saint-Eve, M. Marin, and I. Souchon Mechanistic Mathematical Model for In Vivo Aroma Release during Eating of Semiliquid Foods Chem Senses, February 1, 2008; 33(2): 181 - 192. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. C. Thomas, T. C. Wright, A. Formusiak, J. P. Cant, and V. R. Osborne Use of Flavored Drinking Water in Calves and Lactating Dairy Cattle J Dairy Sci, August 1, 2007; 90(8): 3831 - 3837. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||










