Coffee extraction kinetics in a well mixed system

The extraction of coffee solubles from roasted and ground coffee is a complex operation, the understanding of which is key to the brewing of high quality coffee. This complexity stems from the fact that brewing of coffee is achieved through a wide variety of techniques each of which depends on a large number of process variables. In this paper, we consider a recent, experimentally validated model of coffee extraction, which describes extraction from a coffee bed using a double porosity model. The model incorporates dissolution and transport of coffee in the coffee bed. The model was shown to accurately describe extraction of coffee solubles from grains in two situations: extraction from a dilute suspension of coffee grains and extraction from a packed coffee bed. The full model equations can only be solved numerically. In this work we consider asymptotic solutions, based on the dominant mechanisms, in the case of coffee extraction from a dilute suspension of coffee grains. Extraction in this well mixed system, can be described by a set of ordinary differential equations. This allows analysis of the extraction kinetics from the coffee grains independent of transport processes associated with flow through packed coffee beds. Coffee extraction for an individual grain is controlled by two processes: a rapid dissolution of coffee from the grain surfaces in conjunction with a much slower diffusion of coffee through the tortuous intragranular pore network to the grain surfaces. Utilising a small parameter resulting from the ratio of these two timescales, we construct asymptotic solutions using the method of matched asymptotic expansions. The asymptotic solutions are compared with numerical solutions and data from coffee extraction experiments. The asymptotic solutions depend on a small number of dimensionless parameters, so the solutions facilitate quick investigation of the influence of various process parameters on the coffee extraction curves.


Introduction
Coffee is a popular beverage prepared from the beans (seeds) of the coffee plant. The consumption of coffee, which has been increasing in recent years, is on a global scale. In  coffee production reached . million tonnes worldwide, which represents more than  billion cups []. Typically coffee preparation involves three main stages. First the raw coffee beans are roasted. Following this the roasted beans are ground or milled to facilitate a faster extraction during the final brewing stage. The brewing stage involves the leaching of coffee solubles from the roasted and ground coffee grains with hot water. Generally the extract is filtered to remove undissolved solids from the final coffee beverage before consumption. There are a wide variety of techniques used in the brewing stage. These brewing methods fall into three broad categories: decoction methods, infusion methods and pressure methods. A large number of these techniques is described in [, ]. Despite the wide variety of methods available to brew coffee, each method relies on solid-liquid extraction of coffee solubles from coffee grains with hot water as the central operation. Even in its simplest manifestation the brewing of coffee is a complex operation which is dependent a large number of process variables. Important parameters include the brew ratio (dry coffee mass to water volume used), grind size and distribution, brewing time, water temperature, agitation, water quality and uniformity of extraction [, ].
A fundamental goal for the manufacturers of coffee brewing machines and end users of their products is to be able to consistently brew the highest quality coffee possible. This presents a number of challenges. First of all coffee quality is difficult to define and to some degree a matter of taste. In addition to this, even if the composition of the extracted coffee solution can be identified for different tastes, the coffee machine should be able to change the brewing parameters to move between these ideal compositions. To achieve this goal, both a model relating the final composition of the brewed coffee to the process parameters, and a measure of the quality of the brewed coffee, are required. The complex chemistry of coffee makes it difficult to find correlations between the individual chemical constituents of the extracted coffee solubles and the quality of the final beverage. In the absence of such a description, a set of standards, summarised in the coffee brewing control chart [], is often used as a simple measure of coffee quality. This chart gives target ranges for the brew strength and extraction yield of the coffee based on preferences recorded in organised taste tests. Brew strength is defined as the ratio of the mass of dissolved coffee in the beverage to volume. Extraction yield is the percentage of dry coffee grind mass that has been extracted as solubles into the water. Brew strength and extraction yield are related by the brew ratio. Given that the most widely used measure of coffee quality considers coffee as a single component, it seems logical to model extraction based on a single coffee constituent []. Clearly modelling of the coffee concentration (brew strength) of the final beverage for a particular brewing apparatus, with the process parameters as inputs, could be a valuable tool in order to choose the optimal set of parameters to achieve a targeted coffee quality on the coffee brewing control chart. To achieve this goal for any brewing apparatus, an accurate physical description of the extraction kinetics of coffee solubles from the coffee grains is a key requirement.
The physics of coffee extraction from various coffee brewing systems has been the subject of some investigation and modelling over the years. Industrial scale coffee extraction to produce instant coffee has been studied with the aim of optimising the design of these systems. Early work [, ] looked at models of coffee extraction in large packed columns, called diffusion batteries, with the aim of extracting highly concentrated solution. In these systems, extraction was performed by forcing hot water through these columns to extract solubles. In order to achieve the maximum possible concentration in the extract, extraction was performed at high temperatures and columns were connected consecutively in a series. Much of this work is summarised in []. At a smaller scale, domestic and catering brewing of coffee has also been subject of some investigation. Experiments on the operation and efficiency of the stove-top of Moka pot are detailed in [, ]. Fasano et al. have developed some very general multiscale models of coffee extraction in different sit-uations, with a particular focus on the espresso coffee machine [-]. The influence of some brewing parameters on coffee extraction was investigated by Voilley and Simatos []. A number of different extraction experiments was conducted on a well-mixed system of coffee grounds and water. The response of brew strength to variations in process parameters such as brewing time, granule size, brew ratio and water temperature was considered. A simple model was used to describe the variations of brew strength during the experiments. The model assumed the coffee grains were spherical and suspended in a homogeneous system. Extraction was then modelled as diffusion of a single component from a sphere with a diameter equal to the mean grain diameter. The model was fitted to the data using the diffusion coefficient as a fitting parameter and was found to provide reasonable agreement with the experimental extraction curve. It was noted, however, that the initial extraction proceeded at a faster rate in the experiment compared to the model extraction curve.
More recently, the drip filter brewing system has been the subject of increased attention. In this brewing technique hot water is poured over coffee grounds contained in a filter which is often conical in shape. The hot water flows down through the coffee bed under the influence of gravity, leaching soluble content. The only impediment to flow is the resistance provided by the porous bed and the filter. A schematic of this process is shown in Figure . The drip filter brewing system was analysed by a group of applied mathematicians at the ESGI  study group with industry at the University of Limerick []. The problem was brought to the study group by coffee researchers at Philips Research in Eindhoven. Some of the aspects considered included the evolution of the shape of the coffee bed during brewing, correlations observed between the final bed shape and the quality of the brewed coffee and the use of a single jet or multiple jets (shower head) to deliver water to the coffee bed. Following this, a paper by one of the study group participants [], analysed gravity driven flow in a conical filter, similar to that used in a drip filter machine, and constructed a mathematical model describing the flow. In a recent paper, Moroney et al. [] presented a new multiscale model of coffee extraction from a coffee bed. The model is derived from first principles, by forming balance equations for coffee solubles and water in the different phases in the bed, at different length scales. The parameters in the macroscale model are related to the microscale parameters in the coffee bed by a volume averaging procedure. Transport of water or coffee solubles across the phase boundaries on the microscale appear as source terms on the macroscale. The general model is then specialised based on observations in coffee extraction experiments. In particular it was noted that extraction seemed to proceed in two stages; an initial rapid extraction over a short period, followed by a significantly slower extraction during the rest of the brewing period. This phenomenon was explained by assuming that the fast initial extraction was due to a reduced mass transfer resistance in fine particles and the broken cells on the surfaces of grains following grinding. The slower extraction was assumed to occur due to a much larger mass transfer resistance in the intact coffee cells in the kernels of larger coffee grains. The model was specialised to describe coffee extraction in two experimental situations: extraction from a dilute suspension of coffee grains in water and extraction from a packed coffee bed. Numerical simulations of the model equations showed that the model could quantitatively reproduce the extraction profiles of these experiments for one fine coffee grind and one coarse coffee grind considered.
In this paper the coffee extraction model from [] is analysed in the case of the extraction from a dilute suspension of coffee grains. The analysis in the case of extraction from a packed coffee bed is the subject of a separate paper []. Analysing extraction from a well-mixed, dilute suspension of coffee grains is of interest, as it allows us to investigate the extraction kinetics from the coffee grains, independent of the complications introduced by the advection and mechanical dispersion occurring in pressurised flow through a tortuous packed coffee bed. The general model from [] will be introduced and specialised to describe extraction from a dilute suspension of coffee grains. The dominant mechanisms during the extraction are then identified by non-dimensionalising the equations. Approximate solutions are formed based on the dominant processes during different stages of the extraction. The approximate solutions are compared with numerical solutions of the system and also with experimental data presented in [].

General coffee extraction model
The general coffee extraction model outlined here is described in detail in []. The full model will not be analysed in this paper, but it is useful to summarise the model before specialising to the case of interest here. The coffee bed is represented using a double porosity model. This means that the coffee bed consists of grains which are themselves porous. Temperature variations during brewing are considered to be small, so isothermal conditions are assumed, and the variation in any temperature dependent parameters is considered negligible. At the macroscale the model consists of conservation equations for coffee and liquid in three phases. The phase consisting of the pores between the coffee grains (intergranular pores) is called the h-phase. The pores inside the grains (intragranular pores) make up the v-phase. Finally the solid coffee cellular matrix in the grains is called the s-phase. The main transfers of coffee solubles and fluid between phases in the coffee bed are illustrated in Figure . The coffee concentrations (mass per unit volume) in each of the phases are denoted by c * h , c * v and c s . The solid coffee matrix density, c s , is assumed constant. Variables in the dimensional model which will appear in scaled form in the dimensionless model are denoted by an asterisk. Note that this includes some intrinsically dimensionless variables which are normalised in the dimensionless model. The porosity or volume The diagram shows the transfers of water and coffee which are described by the coffee extraction model presented in [5]. Boundary conditions must be prescribed to describe any transfers of fluid or solubles in or out of the bed where necessary. fraction of the intergranular pores is denoted by φ h and is assumed constant. The volume fraction of the grains (φ h ) is divided into two domains. The intragranular pores have a volume fraction (of the total grain volume) of φ * v , while the solid coffee matrix has a grain volume fraction of φ * s = φ * v . The actual solid coffee volume fraction is further divided into three parts as illustrated in Figure . The volume fraction of coffee which is insoluble under the conditions in the coffee bed is denoted by φ s,i and may depend on water temperature, coffee grind distribution and other variables. Experiments in [] show that for water at °C the extractable mass in coffee grains can vary from % for coarse grinds to % for fine grinds of the same coffee. The soluble coffee grain volume fraction is split into two parts according to its position in the coffee grains. Fines, which are broken cell fragments produced during grinding of the coffee beans, can account for a significant volume of the coffee grind distribution, particularly for finer grinds. Following grinding the surfaces of the coffee grains consist mainly of broken coffee cells. Coffee in each of these regions in the coffee bed is expected to have a significantly lower mass transfer resistance relative to that in intact cells in the kernels of larger grains. The volume fraction of this coffee is denoted φ * s,s . The intact cells in the kernels of larger grains have a much higher mass transfer resistance. The volume fraction of this coffee is denoted φ * s,b . Thus the volume fraction of soluble coffee in the coffee grains is It is also useful to track the fraction of the original amount of coffee (in the dry coffee grains) present in the grains surfaces and in the grain kernels at a given time during brewing. The volume fractions of coffee in the dry coffee grains, grain surfaces and grain kernels are given by φ cd , φ s,sd and φ s,bd respectively. The fractions of the original amount of coffee left on the grain surfaces and in the grain kernels are denoted by ψ * s and ψ * v . In extractions involving flow through a bed of grains the surface area of the coffee grain distribution is a key parameter influencing the flow. The surface area is also key to extraction of coffee during the brewing process. The specific surface area of the entire grain distribution, which can influence flow and extraction of coffee from grain surfaces, is described using the Sauter mean diameter k sv . This is the diameter of a spherical particle which has the same surface area to volume ratio as the entire grind size distribution []. The specific surface area of larger grains (diameter larger than  μm), which influences extraction from the kernels of these larger grains, is described using the Sauter mean diameter k sv . The effective diffusion coefficient of coffee in water is given by D h in the h-phase and D v in the v-phase. The effective diffusion distance between the v-phase and the h-phase is given by l l . The average coffee cell radius in the grains is given by m. It is assumed that there is some coffee concentration c sat , which is the concentration in the liquid phase that would be in equilibrium with the concentration in the solid. Here we take it to be the maximum solubility of coffee in the liquid. Mechanical dispersion can be an important factor in flows in porous media []. The dispersion tensor is denoted byD b . Other important parameters necessary to complete the description include the fluid density ρ, the dynamic viscosity of the fluid μ and the shape factor κ from the Kozeny-Carman equations. The system of equations in [] included a term for transfer of fluid (and coffee solubles as a result) between the h-phase and v-phase due to differences in pressures. This correction is difficult to model accurately without further experimental insight and is thought to occur much faster than other bed processes. For this reason these terms are neglected in the presentation of the equations here. Neglecting these terms, the general mathematical model becomes The constants r s and r v are the reciprocals of φ s,sd and φ s,bd respectively. Boundary conditions need to be prescribed depending on the geometry of the coffee bed considered. Initial conditions need to be determined or inferred from experiment once the bed is saturated with water, following the initial addition of water to the dry coffee bed. Alternatively, modelling of the unsaturated flow during the initial infiltration of water could be used. A complete description of equations ()-() is presented in []. The exact form of the coefficients in terms of the process parameters is arrived at through a volume averaging procedure. This process is detailed in the appendices of []. For the purposes of this paper, it is sufficient to point out the meaning of each of the terms in the equations, in the context of the transport process taking place in the system. Equation () describes how the concentration of coffee in the interstitial fluid (intergranular pores) is changing during brewing. The first term on the right-hand side of the equation describes the advection of coffee in the fluid flow based on Darcy's Law and the Kozeny-Carman equations. The second and third terms describe diffusion and mechanical dispersion of the coffee solubles relative to the flow. The fourth term describes transport of coffee into the interstitial fluid from the fluid in the intact cells in the grain kernels. The final term describes the rapid dissolution of coffee from the fines and broken cell surfaces directly into the interstitial fluid. Equation () gives the equation for flow of the fluid in the intergranular pores in a coffee bed according to Darcy's Law. Equation () describes the changes in the quantity of coffee solubles in the intragranular pores during brewing. Note here that along with changes in the coffee solubles concentration c * v , the intragranular porosity φ * v may also change as coffee dissolves from the cell walls within the grains. On the right-hand side of the equation, the quantity of coffee in the intragranular pores is decreased by slow diffusion of coffee out of the grains (first term), but increased by coffee dissolving from the cell walls (second term). The diffusion length scale has been adjusted for grain tortuosity τ v , using the functional Equations () and () keep track of the fraction of the initial amount (in dry coffee) of coffee present on the grain surfaces and in the grain kernels at a given time respectively. The terms on the right-hand side of these equations model the dissolution of solid coffee in these two regions. As this coffee dissolves the intragranular porosity grows from its initial value. This increase in intragranular porosity is described by (). In the following section this dimensional system of equations will be specialised to model experimental results for extraction of coffee from a dilute suspension of coffee grains. The model will be non-dimensionalised and analysed to develop approximate solutions.

Coffee extraction model for coffee grains in a fixed volume of water
One of the experiments described in [] involves mixing  g of coffee grind with . L of hot water and measuring the concentration of the extracted species as a function of time. In this situation the coffee extraction model can be simplified. Firstly it is assumed that the solution in the liquid outside the grains (h-phase) is well mixed, since only the average concentration is measured anyway. As a result, the solution in the h-phase has uniform concentration. This means that the spatial derivatives in the model disappear and we have the following system of ordinary differential equations: The last term now appearing on the right-hand side of equation () is a volume correction term to allow for changes in the intragranular porosity. Additional parameters α * , β * and γ * are introduced here to fit the model to experiment. These parameters are necessary due to the difficulty in obtaining accurate values for some of the model parameters.
For example the diffusion coefficient used for coffee solubles will be for that of caffeine in water. However caffeine diffusivity may differ from the average coffee diffusivity in water and the extraction of many coffee components from the grains cellular matrix is hindered to a greater degree than that of caffeine []. Parameters relating to coffee grain properties, such as the cell radius and external and particularly internal specific surface areas are similarly difficult to measure. Thus fitting parameters are required. Clearly to solve this system we need to prescribe initial conditions for each of the dependent variables. This model does not include the addition of the water to the dry coffee grains or the infiltration of water into the intragranular pores. Rather it is assumed that this infiltration occurs quickly and the model starts with the grains saturated with water. This raises the question as to what initial conditions we should use for the model. How much of the coffee has dissolved during the addition of water to the coffee and the infiltration of water into the coffee grains? In the absence of experimental data or separate modelling of this initial filling stage we will need to prescribe appropriate initial conditions. The simplest choice would be to assume that no coffee dissolves during this filling stage so that all initial concentrations are zero but this may not be the most useful.
To reduce the steps we have to take when non-dimensionalising the problem it is useful to make some simplifying steps at this stage based on our knowledge of the problem. First of all the sum of the intragranular pore volume fraction, and the soluble coffee volume fraction is conserved. This is described by equation (). Integrating this equation with respect to time we obtain with K a constant. We can determine this constant by denoting the final intragranular porosity when all soluble coffee has been dissolved as φ ∞ v . The value of φ ∞ v can be estimated from data on the percentage of soluble coffee for a particular coffee grind. Thus we have Using this equation and () we can eliminate φ * v from the model system of equations. We can further simplify the system by noting that the change in intragranular porosity during extraction is small relative to the final (or initial) value of intragranular porosity. More Thus we make the approximation that φ * v ≈ φ ∞ v in the equations. Initial conditions should be adjusted to account for this. Also the volume correction term now drops out of equation ().
Finally we consider the fitting parameters α * , β * and γ * . Experimental data for c * h obtained in [] reveals a rapid initial extraction followed by a much slower extraction. We can use this data to fit α * and β * . As dissolution of coffee within the grains and diffusion of coffee from the grains occur consecutively, fitting the parameter α * rather than γ * corresponds to the assumption that diffusion of coffee from the grains is the rate controlling step in this series. The fitting parameter γ * for the rate of dissolution of coffee from the cell walls within the grains cannot be determined independently from α * from experimental data available in []. This is because the intragranular concentration c * v is not measured experimentally. A reasonable guess however may be that it has a value similar to the surface dissolution fitting parameter β * . In reality, due to a larger surface area per unit volume within the coffee grains, it is likely γ * > β * . The timescale on which surface dissolution occurs is much shorter than that of diffusion of coffee from the grain kernels. Also based on the assumption γ * > β * , dissolution from cell walls within the grains occurs on a similar timescale to surface dissolution. Thus diffusion through grain cellular structure is the rate controlling step in the extraction of solubles from the intact cells in the grain kernels. This suggests that we can approximate these processes by assuming that all the coffee in the grain kernels dissolves very quickly initially into the intragranular pores and extraction is just controlled from there. This of course means that the concentration in the intragranular pores and the intragranular porosity will be incorrect at least for a short period initially but in any case we have no data for comparison. The approximation has the advantage of giving us easy to calculate initial conditions. To ensure conservation of coffee during extraction we need to be careful when prescribing initial conditions. Finally we assume that initially none of the coffee on the grain surfaces has dissolved. Based on this discussion we choose the following dimensional initial conditions: Since ψ * v () = , we have reduced the system to the following three ordinary differential equations:

Non-dimensionalisation
There are three main timescales evident in the equations: The timescale on which the coffee dissolves within the grains, the timescale on which coffee dissolves from the grain surfaces and finally the timescale on which coffee diffuses from the intragranular pores to the surrounding fluid. We have already made the assumption that the diffusion of coffee from the intragranular pores occurs on a much longer timescale than dissolution of coffee within the grains. On the basis that these processes occur consecutively, only diffusion of coffee from the intragranular pores is considered as the rate limiting step. Thus there are two remaining timescales to consider. We denote the bulk diffusion timescale by t d and The scale for ψ * s can be chosen either from the initial condition ψ * s () = ψ * s =  or by balancing the source term for the transfer of coffee from the grain surfaces with the other two terms in equation (). Either of these scales allow us to solve the problem so we adopt the former which is more tractable. Thus the scales are Before the dimensionless equations are presented we tidy up the presentation by introducing some dimensionless parameters. First the ratio of the surface dissolution timescale to the timescale of diffusion of coffee from the grain kernels is denoted by For parameters of both fine and coarse grinds the parameter , while b  and b  are both O(). We represent the non-dimensional initial concentration of coffee in the intragranular pores by the parameter γ  . Thus our equations and initial conditions on the diffusion timescale are On the long timescale, this suggests that at leading order s relaxes to a steady state zero value, while C v decays exponentially with time. It is clear from the form of the equations that we have a singular perturbation and so initial conditions may not be applied to the outer solutions. For this system we have a number of ways to solve the problem. To begin with we consider the solution of the problem in a phase plane.

Phase plane analysis
In the case of extraction of coffee from a suspension of coffee grains, unlike the system for extraction from a packed coffee bed, we have conservation of coffee in the system. Thus we can reduce the number of equations in the system by one by using this property. We note that and integrating with respect to τ gives with C constant. Using the initial conditions we can write Thus s can be eliminated from the system to give with initial conditions given by (). Writing the equations in this form has some advantages. The number of equations is reduced to two, which facilitates the investigation of the solutions for C h and C v on a phase plane. We find the nullclines by setting the right-hand sides of equations () and () equal to zero. The physically relevant region of interest in the problem is Once the nullclines are drawn on the phase plane it is relatively straightforward to plot the solution on the plane. The phase planes and solution trajectories for the fine and coarse grind experiments are shown in Figure . In both cases, we see a rapid transition from the initial point to the C h nullcline before a much slower transition along the slow manifold towards a stable equilibrium point at the intersection of the nullclines. This stable equilibrium is given by It is clearly useful to be able to reduce the number of equations to two, but it is more difficult to understand the physical meaning of the individual terms when the equation is written in this from. Thus when proceeding with the determination of approximate solutions of the problem the system of three differential equations will be used.

Perturbation solutions on the bulk diffusion (outer) timescale
We use the following expansions for the bulk diffusion timescale Substituting these expansions into the equations ()-() collecting terms at each order, it is straightforward to find the solutions on the bulk diffusion timescale. As mentioned above the problem is singularly perturbed. Thus these outer solutions will involve constants which need to be determined by matching solutions obtained from the rescaling of the problem to the initial layer.

.. Leading order equations
The leading order equations are As expected we see that s = . Solving the second equation we find that

.. Order equations
Substituting in the known terms, the order equations are Once again we see immediately that s = . Solving the remaining two equations we find that .. Order  equations Substituting in the known terms, the order  equations are Solving the non-trivial equations we find that

.. Outer solutions
Collecting the outer solutions we have Thus the composite solutions on the inner timescale and concentration scale are

Comparison of asymptotics with numerical simulation
The composite solutions are compared with the numerical solution of equations ()-() and available experimental data in Figure . The comparison with data is shown for two different grind size distributions from []: one fine grind called JK drip filter grind and one relatively coarse grind called Cimbali # grind. Available data for coffee concentration in the intergranular fluid is compared to the approximate and numerical solutions based on the parameters for these two grinds. The asymptotic and numerical solutions to the model for intragranular concentration c v (t) and the fraction of coffee remaining on the grain surfaces ψ s (t) are also plotted. No data is available for comparison in this case.

Conclusion
In this paper the mathematical model of coffee extraction described in [] is analysed. The model equations are simplified and non-dimensionalised to describe extraction of coffee from a dilute suspension of coffee grains. The dimensionless form of the equations represents a much simpler description of the important processes in this situation and depends on a much smaller set of parameters than the original equations. This is a good example of model reduction []. Approximate solutions of the specialised dimensionless equations can be found using perturbation techniques. These solutions depend on dimensionless parameters which are defined in terms of the physical parameters of the system. Thus the approximate solutions are useful to both quickly fit solutions to a particular set of process parameters and to investigate the influence of changing a particular parameter on the system. The model equations used to describe the extraction of coffee in a fixed water volume admit approximate solutions due to the existence of two different extraction mechanisms over different timescales. Extraction of coffee from the surfaces of the coffee grains proceeds much faster than diffusion of coffee through the intragranular pore network to the grain surface. The ratio of these timescales gives us a small dimensionless parameter. We utilise this parameter to form solutions on the inner (surface dissolution) and outer (diffusion from grain kernel) timescales based on the dominant mechanisms in these regimes. These approximate solutions are found to match the numerical solutions well and provide a much simpler equation to fit to data.
Extraction of coffee solubles from individual coffee grains is a key operation in many brewing techniques. The physical description of extraction here is observed to describe extraction well in the case of extraction from a dilute suspension of coffee grains. Thus it may be useful in models of more complicated brewing techniques to describe the grain extraction kinetics component of the procedure. This has already been shown in the case of extraction in flow through a packed coffee bed []. There is a wide variety of topics in the area of coffee extraction and coffee brewing in general which may benefit from further mathematical modelling and investigation. The model in [] may be generalised to describe the unsaturated flow during the filling and draining of a coffee filter. The model can also easily be extended to model the extraction of a number of different coffee constituents rather than just a single entity. This may be useful if the influence on flavour of particular constituents (or groups of constituents which extract at similar rates) can be identified. The modelling of coffee in a drip filter presents a number of challenges. The filter geometry needs to be accounted for, the fluid flow in the filter may be complex and the coffee grains can be transported around the filter in the flow. On the scale of a coffee grain there is scope to investigate the dissolution of coffee within a coffee grain in more detail to investigate some of the assumptions made in this paper.