The behaviour of foams at rest, but particularly during fluid mechanical transport is not sufficiently investigated yet. The present article deals with protein foams as they have a great importance in food production. In the first part, the foaming process of a highly viscous liquid due to gaseous materials dispersed under pressure in the liquid and mass transport of volatile components dissolved in the liquid is considered. The aim is to calculate the foam volume and the concentration of the dissolved, volatile components as a function of the material and process parameters. In the second part, material equations for bubble suspensions with gas volume fractions φ ≤ 0.6 and small bubble deformations (i.e. NCa 1) are presented. The basics form two constitutive laws which are used for describing a steady shear flow. If the rates of work of the two models are compared, material equations for the shear viscosity and the normal stress differences can be derived.
Foaming processThe model presented in [4] and developed further in this article describes the foaming process of a highly viscous liquid due to gaseous materials dispersed under pressure in the liquid and mass transport of volatile components dissolved in the liquid. The model is divided into two model phases, which describe spherical foam and polyhedron foam, respectively. The increase of the foam volume is calculated by solving problem-specific conservation equations for mass, momentum and energy. Alike, concentration equations for the dissolved, volatile components, which are diffusing into the bubbles during the foaming process, are solved. The equations are subjected to a Lagrangian coordinate transformation and transformed into dimensionless equations by means of a set of dimensionless π-numbers. The equation of motion for the spherical foam is given bywhere R B is the bubble radius, R S is the radius of the liquid volume surrounding the bubble, P g is the internal pressure of the bubble, P ∞ is the ambient pressure, T Y Y and T θθ denote the normal stresses, T is the time and Y is the transformed spatial coordinate. The integral on the right side of equation (1) characterises the viscous elongation stresses around the bubble. The rheological behaviour of the liquid is described by the constitutive equation for Newtonian liquids or the generalised contravariant Maxwell-Oldroyd model for viscoelastic liquids. In the latter case, the extra stress tensor T results from the sum of the partial stresses T k based on a discrete spectrum that covers relaxation times λ R,k and viscosities η k , the sum of which is equal to the zero shear viscosity η 0 of the liquid:The function a takes account of viscosity changes due to changes in the temperature Θ and the concentration ξ i of the volatile components dissolved in the liquid. A similar differential equation is derived for the other relevant normal stress component T θθ of the extra stress tensor T. The equations for the partial stresses and the equation (1) are combined, resulting in a large system of non-linear ...