Next Article in Journal
Fermented Sargassum fusiforme Mitigates Ulcerative Colitis in Mice by Regulating the Intestinal Barrier, Oxidative Stress, and the NF-κB Pathway
Next Article in Special Issue
Unraveling the Mechanism of StWRKY6 in Potato (Solanum tuberosum)’s Cadmium Tolerance for Ensuring Food Safety
Previous Article in Journal
Effect of Olive Pomace Extract Application and Packaging Material on the Preservation of Fresh-Cut Royal Gala Apples
Previous Article in Special Issue
Analysis of Dynamics and Diversity of Microbial Community during Production of Germinated Brown Rice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ultrasonic Wave Mode-Based Application for Contactless Density Measurement of Highly Aerated Batters

Chair of Brewing and Beverage Technology, TUM School of Life Sciences, Technical University of Munich, 85354 Freising, Germany
*
Author to whom correspondence should be addressed.
Foods 2023, 12(9), 1927; https://doi.org/10.3390/foods12091927
Submission received: 31 March 2023 / Revised: 5 May 2023 / Accepted: 7 May 2023 / Published: 8 May 2023
(This article belongs to the Special Issue Recent Advances in the Food Safety and Quality Management Techniques)

Abstract

:
An ultrasonic wave mode-based method for density measurement in highly foamed batters was developed. Therefore, a non-contact ultrasonic sensor system was designed to generate signals for batch-wise processes. An ultrasonic sensor, containing a piezoelectric ceramic at the fundamental longitudinal frequency of 2 MHz, was used to take impedance measurements in pulse-echo mode. The ultrasonic signals were processed and analysed wave-mode wise, using a feature-driven approach. The measurements were carried out for different mixing times within a container, with the attached ultrasonic sensor. Within the biscuit batter, the change to the ultrasonic signals caused by density changes during the batter-mixing process was monitored (R2 = 0.96). The density range detected by the sensor ranges between 500 g/L and 1000 g/L. The ultrasonic sensor system developed also shows a reasonable level of accuracy for the measurements of biscuit batter variations (R2 > 0.94). The main benefit of this novel technique, which comprises multiple wave modes for signal features and combines these features with the relevant process parameters, leads to a more robust system as regards to multiple interference factors.

1. Introduction

The quality of baked products depends on many factors. To ensure consistently high quality, the production process must be capable of adapting to the fluctuating qualities of the raw materials [1]. Even the short-term storage of flour causes changes to its baking characteristics [2]. Quality-related attributes, such as the specific volume of cake or the porous structure of crumbs, depend on both the cake batter recipe and the general baking conditions [3,4,5].
For highly foamed baked goods, such as biscuits, the air content in the batter affects the texture and quality of the final product significantly. The structure of the final product, however, depends heavily on the process used to mix the cake batter, whereby this constitutes one of the major processes in any bakery system. During mixing, air is introduced during a liquid phase, which consists of egg, sugar, wheat flour and other ingredients, to form the batter. As there is no additional raising agent, the final structure of sponge-cake depends on the air trapped within the continuous phase [6,7] and is, therefore, related to its density [8]. During mixing, the density changes from 1000 g/L (poorly aerated) to under 600 g/L (highly aerated). The mixing time thus depends on the density. This is because of the amount of air bubbles introduced into the batter during the process increases and, ultimately, ends up affecting the cake volume. Therefore, the density of the batter must be measured [9]. To monitor critical process steps such as mixing, real-time data are required: such data afford a deeper insight into the current state of the process. Thus, online sensors, which are able to take measurements directly within the process as well as inline sensors measuring bypass or required samples, are suitable for this purpose [10].
Within the food industry, there is a high interest in online systems measuring the density and the textural properties of food products. Currently, the usual procedure for density measurement uses a gravimetric offline measurement of the batter. It does this by analysing samples of fixed volumes. Because of this analogue, time-consuming step, there is a need to develop a non-invasive, online method for determining the density in baking processes.
Another common method used to study the mixing process is based on measuring work input and torque, using a strain gauge. The mixing process is stopped at a predefined total energy input, which is calculated beforehand by integrating the power over time [11,12]. Using this method, however, only non-specific physical modifications can be analysed, while the method itself depends heavily on the type of vessel and the mixture employed. Other direct rheological measurements were performed for different dough systems, such as yeasted bread doughs [13], sour doughs and cake batters. The mixing process can also be investigated by means of X-ray CT (computed tomography) measurements. The detection of bubble size and oxygen dispersion within a dough mix has already been shown [14]. This method has a high level of accuracy but is not suitable for inline measurement, due to the complex evaluation process it necessitates.
Ultrasonic-based sensors are perfectly suited to measurement systems operated during food processing, due to their hygienic, non-invasive and non-destructive behaviour [15,16,17]. The use of ultrasonic-based measurements for solid dough compounds is already well-established [18]. It is already possible to measure wheat quality [19] and the precise composition of wheat-flour–water doughs [20] using ultrasonic methods. The influence of the properties of bread dough during mixing [20,21,22] and fermentation [23,24] as detected by ultrasonic signals has been shown. Recently, it was also shown that it is possible to monitor the mixing process for liquid model systems using ultrasound [25]. It is even possible to detect air bubbles in dough with ultrasonic methods [26,27]. The impact of air within a more liquid-like product seems to be more challenging, however. There are only a few reports of ultrasonic-based investigations of batters [28]. Because of the higher air content in batters compared to dough systems, ultrasonic signals are dampened almost completely within the medium and are therefore not suitable for transmission. Hence, only ultrasonic measurement methods that do not depend on acoustic transparency can be used.
Ultrasonic waves for measurement systems are generated mostly by piezo-ceramic transducers, which convert electrical energy into mechanical energy and vice versa. The measurement method can be divided up into three basic components: reflection, transmission and emission.
Reflection measurements, using either one (pulse-echo mode) or multiple transducers (pitch-catch mode), send and receive the reflected wave at a given interface. Transmission measurements are realized with a sending and receiving transducer, where the ultrasonic signal is transmitted through the material, whereas for emission measurements, the sound signals recorded are emitted by the process itself. Ultrasonic waves for measurement applications are mostly low power and high frequency, in order to prevent changes to the media they are transmitted through [15]. In the reflection mode, impedance, as measured at the surface of the batter [9], correlates with density.
Other studies have also shown the feasibility of ultrasonic techniques for density measurement in batter [9,29,30,31] but only for poorly aerated batters (~1000 g/L). For highly aerated batters, however, the high reflection and distortion rates of ultrasonic waves in gas bubbles do not allow transmission measurements. Therefore, only reflected, scattered waves can be used for analysis purposes. Furthermore, the ultrasonic waves can be influenced by reflection and scattering at each interface used. The resulting ultrasonic signals interact at all of the interfaces individually and are then subjected to mode conversion [32]. Mode conversion occurs when an ultrasonic wave faces an interface between materials of different impedances and the incident angle to the interface is abnormal. Some of the sound energy introduced is therefore reflected at the interface between the buffer rod and the surface of the batter, due to the acoustic impedance. The acoustic impedance Z can be described with Formula (1).
Z = p u = ρ c
p : sound pressure, p : sound particle velocity, 𝜌: density, c : ultrasonic velocity, Z : impedance.
Due to the different impedances of the two different media (buffer rod—batter), reflections occur at their mutual interface. The reflection coefficient R can be described with Formula (2).
R = p d p e = Z 2 Z 1 Z 2 + Z 1
Z 1 : acoustic impedance of medium 1, Z 2 : acoustic impedance of medium 2, p e : emitted sound energy, p d : detected sound energy, R : reflection coefficient.
The ratio between the emitted sound energy pe and detected sound energy pd depends on the acoustic properties of both interfaces, as represented by the acoustic impedance Z1 of the medium M1 and Z2 of the medium M2. Due to the differing impedances of the phases (solid/liquid/air), a mode conversion also occurs between longitudinal, shear and surface–acoustic waves [33,34]. At the changing interfaces between the batter structure and the air bubbles within the batter, the original ultrasonic waves scatter and meet the air bubbles at different angles. Due to the characteristic of batter being a non-Newtonian fluid (i.e., neither solid nor liquid), the ultrasonic waves are also mode converted within the batter (longitudinal to shear waves) [35].
As regards the interpretation of ultrasonic signals, the use of special analytical strategies plays a decisive role. There are two main kinds of approaches to predicting the desired process parameters by means of ultrasonic signal data: model-based approaches and machine-learning approaches.
Model-based approaches rely on expert knowledge of the processes and the underlying physical relations. Gaining expert knowledge of the whole system is a very complicated undertaking, while, ultimately, due to their inherent complexity, both the process and the production elements here are affected by a large number of factors [36].
On the other hand, the machine-learning approaches do not require expert knowledge but rather, a considerable amount of sensor data, in order to create a decent model. For modelling complex processes, using only expert knowledge-based approaches is very difficult and, therefore, not possible for many process steps. The downside of the machine-learning approach using deep learning algorithms is that a large amount of data is required in order to come up with a proper prediction. The ultrasonic features of the approach used in this study are based on the physical behaviour of the biscuit masses. Thus, expert knowledge can be integrated into the machine-learning approach via ultrasonic features.
This work shows an ultrasonic-based, non-invasive low-cost method for determining the density of highly aerated batter during a batch process. The method therefore uses a feature analysis, which is based on impedance measurement regarding mode conversion, to correlate ultrasonic signals with batter density. This also applies to different recipes (e.g., standard, gluten free).

2. Materials and Methods

2.1. Materials

Wheat flour (Type 550) was supplied by Rosenmehl, Rosenmühle GmbH (Ergolding, Germany). Pasteurized whole egg, white sugar and wheat starch were purchased from the local market. The emulsifier (Colco) was supplied by Aromatic (Aromatic Marketing GmbH—Berlin, Germany).

2.2. Batter Formulation

A standard procedure based on industrial production was used for the batter formulation. A set of four batters was prepared with different sugar, flour, wheat starch and emulator contents (see Table 1). The standard mixture was varied (with/without starch/emulsifier) in order to obtain both different batter properties (standard and gluten free) and thermally fixed cakes with varied physical characteristics as a result of varying aeration. All of the ingredients, which were tempered to room temperature, were pre-mixed at a low speed for 120 s. The foaming process was started afterwards by high-speed mixing for various durations, using a KitchenAid Professional mixer 5KPM5 (KitchenAid, Benton Harbor, MI, USA). The temperature measured after the completion of the foaming process was 22 ± 2 °C.

2.3. Density Measurement

The batter was characterized after the completion of the foaming process. The measurements were performed at a weight of 1 L (CSM Deutschland GmbH, Bingen, Germany), with 100 cm3 at a temperature of 22 ± 2 °C. The cup was filled with batter and levelled. The weight was determined using a laboratory balance (Kern, accuracy to ±0.5 mg). To determine the density, the average of three samples was taken and calculated, using the ratio of the mass of batter and the volume of the litre weight. The density can be described using Formula (3).
ρ = m l i t e r   w e i g h t + b a t t e r m l i t e r   w e i g h t 100   c m 3 10
m l i t e r   w e i g h t : mass of empty litre weight [g]; m l i t r e   w e i g h t + b a t t e r : filled litre weight with batter [g]; ρ : density of batter [g/L]

2.4. Ultrasonic Measurement System

The ultrasonic-based measuring system developed consists of a microcontroller for signal generation, as well as a signal recording element. The data obtained from the ultrasonic transducer was stored and evaluated, using Virtual Expert 4.0 software (Gimbio GmbH, Freising, Germany). The self-developed round-shaped ultrasonic transducers with a diameter of 12 mm, which are based on a piezoceramic (PIC 255, PI Ceramic GmbH, Lederhose, Germany), used a median frequency of 2 MHz and were encapsulated in resin and tungsten-based backing material. This was to avoid ringing and to guarantee an adequate signal resolution. The piezo transducer was excited by a rectangular 100 V excitation of 250 ns, with three consecutive pulses, in order to obtain stable signals.
Because of the heavily whipped batters, it was difficult to obtain the transmission measurements. Therefore, the measurement principle was based on the pulse–echo method. Thus, with the pulse–echo measurement method used, the exiting transducer and the receiving transducer were the same. Directly after sending the electrical impulse to the transducer, the set-up switched from send to record. The ultrasonic signals returning were recorded for 60 µs. To avoid the product contact and to guarantee measurements in the far field, the ultrasonic signals were introduced into the medium by an ultrasonic transducer via a polymethyl methacrylate (PMMA) buffer rod (see Figure 1A). Because the pressure wave originated from several points along the piezo transducer disc, the waves interfered in the region directly in front of the sensor (near field). Therefore, the pressure field was not constant up to a certain distance (far field). The natural focus (transition between the near field and far field) can be calculated using Formula (4) [37].
N D 2 4 λ = A π λ
N: natural focus length (mm), D: ultrasonic transducer diameter (mm), A: sensor area (mm2), λ: wavelength of ultrasonic wave propagation in the medium (mm).
Figure 1. Schematic representation of the ultrasonic-based measuring system. (A) The ultrasonic waves are propagating through the buffer rod and are reflected at the surface between buffer-rod/batter and batter/bubbles. (B) Sectional view of the measurement set-up.
Figure 1. Schematic representation of the ultrasonic-based measuring system. (A) The ultrasonic waves are propagating through the buffer rod and are reflected at the surface between buffer-rod/batter and batter/bubbles. (B) Sectional view of the measurement set-up.
Foods 12 01927 g001
The wavelength can be calculated using Formula (5) [37].
λ = c f
c: wave velocity of ultrasonic wave propagation in the medium (m/s), f: ultrasonic frequency (Hz).
For PMMA at room temperature, the ultrasonic wave velocity for 2 MHz is approximately 2740 m/s [38]. The far field calculated for a 12 mm piezo disc sensor for PMMA was thus 26.28 mm. The buffer rod therefore had a length of 30 mm, in order to ensure a measurement in the ultrasonic far-field. The 2 MHz ultrasonic sensor was attached at the bottom of the measurement chamber, which could be filled easily from above (see Figure 1B).
The acoustic properties of batter are subject to change. This is due to batter whipping whenever air is introduced and to the consequent effect on the reflection coefficient R (see Formula (2)). The signal energy, which penetrates the batter, is reflected and scattered at gas bubbles and other surfaces of the batter, as well as mode converted. Therefore, considering the reflection coefficient alone is not enough to determine the density of highly aerated batters.

2.5. Signal Pre-Processing

To eliminate noise interference, e.g., due to pumps or motors, the ultrasonic signals were analysed and filtered by an 8th order Butterworth band pass filter with band threshes at 1.75 and 2.25 MHz (see Figure 2A,B). The signals were collected using Virtual Expert 4.0 software at a sampling rate of 50 MHz for 60 µs. For analysis, the signal-processing toolbox within MATLAB (Version 2018a) was used. Due to the pulse–echo measurement, the signal could be divided into four sections in the time domain (see Figure 2B):
  • Sensor ringing (0–10 µs)
  • Echo of impedance signal (PMMA/batter) (10–30 µs)
  • First echo of reflections and scattering effects in batter (30–50 µs)
  • Second echo of reflections and scattering effects in batter (50–60 µs)
Figure 2. Schematic representation of the ultrasonic signal for evaluation by means of feature analysis. (A) Ultrasonic raw signal, (B) filtered ultrasonic signal divided into the features areas, (C) rank of importance of the features and selection of highest importance weighting features (above red line).
Figure 2. Schematic representation of the ultrasonic signal for evaluation by means of feature analysis. (A) Ultrasonic raw signal, (B) filtered ultrasonic signal divided into the features areas, (C) rank of importance of the features and selection of highest importance weighting features (above red line).
Foods 12 01927 g002
The sensor ringing in the first section was independent from the actual measurement and thus was not considered. The remaining time sections (Section 1, Section 2 and Section 3) were analysed subsequently via ultrasonic features. Since frequency shifts can occur in ultrasonic-based mode conversion, as well as with scattering and reflection, the ultrasonic signal was also analysed in the frequency domain. Therefore, the signal was converted by means of a discrete Fourier transformation [39] (Formula (6)).
X k = n = 0 N 1 x n e 2 π i N k n , k = 0 , , N 1
The Fourier transformation was calculated for the combined sections of the time domain signal and used subsequently for the frequency domain features.

2.6. Ultrasonic Feature Creation

Suitable signal parameters describing the density change were generated from the pre-processed signal data. Physically related parameters, such as acoustic impedance or signal energy, as well as non-specific parameters from the time and frequency domains, were considered to predict the density changes. Because of the major temperature dependencies of ultrasonic propagation, the process temperature was monitored and included in the analysis. The stirring time as a physically based variable was also included as a feature. The following features were adapted from the field “speech, music and environmental sounds” [40] and used for the analysis.

2.6.1. Time Domain Features

MaxAmp shows the highest amplitude value of the signal (Formula (7)).
maxAmp:
m a x A m p = m a x x n
The temporal crest factor (tcf) describes the ratio between the maximum amplitude and the average amplitude and represents the singularity of the signal (Formula (8)).
tcf (Temporal crest factor):
t c f = m a x x n 1 N n = 1 N x n
Temporal centroid is the time at which half the energy of the signal is reached. It is not necessary for this time to be at the zero point (Formula (9)).
tempCentroid:
C t = n = 1 N n · x n 2 n = 1 N x n 2
Temporal energy shows the energy of the whole signal, described by the area under the signal curve by summing up the absolute values of the signal (Formula (10)).
tempEnergy:
t E = n = 1 N x n 2
Temporal skewness is a measure of the asymmetry of signal amplitudes, based on the normal distribution (Formula (11)).
tempSkew:
m 3 = n = 1 N n C t 3 · x n n = 1 N x n t s p = n = r 1 r 2 n C t 2 · x n n = 1 N x n T s = m 3 t s p 3 / 2
Temporal smoothness shows the variation of the signal amplitude related to the direct neighbours of a signal point (Formula (12)).
tempSmooth:
t m p S m = 20 · n = 2 N 1 l o g x n l o g x n 1 + l o g x n + l o g x n + 1 3
The zero crossing rate is the number of points at which the amplitude crosses the zero point (Formula (13)).
zeroCrossing:
z C r o s = 1 2 · N n = 2 N s i g n x n s i g n x n 1 s i g n = 1 ; x n > 0 s i g n = 0 ; x n = 0 s i g n = 1 ; x n < 0

2.6.2. Spectral Domain Features

Spectral attenuation shows the energy of the entire frequency spectrum, as described by the area under the spectrum curve by summing the values of the signals and the subsequent logarithmic conversion (Formula (14)).
attenuation:
a t t = 10 · log 10 m = 1 N X ( m ) ²
The bandwidth shows the bandwidth around the main frequency that is above a certain threshold value (Formula (15)).
bandwidth:
b w = Δ f = f 2 f 1
The spectral centroid is the frequency at which half the energy of the signal is reached (Formula (16)).
centroid:
C s = n = 1 N n · x n 2 n = 1 N x n 2
The spectral crest factor describes the ratio between the maximum frequency and the average frequency and represents the singularity of the frequencies (Formula (17)).
crestFactor:
C r e s t f = m a x x n 1 N n = 1 N x n
The spectral energy shows the energy of the whole frequencies, as described by the area under the frequency curve by summing up the values of the Fourier-transformed signal (Formula (18)).
energy:
E n e r g y = n = 1 N x n 2
The spectral entropy shows the total loss of frequency energy and calculates the loss by summing the energy of the frequencies point by point (Formula (19)).
entropy:
E n t r o p y = n = 1 N x n n = 1 N x n 2 l n x n n = 1 N x n 2
The spectral kurtosis is a measure of the flatness (slope) and the singular variation of the frequency spectrum around the frequency maximum (Formula (20)).
kurtosis:
K u r t = m = 1 N m C t 4 · x m m = 1 N x m 2 m = 1 N m C t 2 · x m m = 1 N x m 2 2 < 3 f l a t t e r   d i s t r i b u t i o n = 3 n o r m a l   d i s t r i b u t i o n > 3 p e a k e r   d i s t r i b u t i o n
The rollOff shows the frequency at which 90% of the energy of the frequencies has already been received (Formula (21)).
rollOff:
r o l l O f f = m = 1 m r f X m = 0.9 · m = 1 1024 X m
The spectral skewness is a measure of the asymmetry of the frequency pattern, based on the normal distribution (Formula (22)).
Skewness  (Ts):
m 3 = n = 1 N n C t 3 · x n n = 1 N x n t s p = n = r 1 r 2 n C t 2 · x n n = 1 N x n T s = m 3 t s p 3 / 2
The spectral slope shows the decrease in the frequency energy over the frequencies by linear regression (Formula (23)).
slope:
S l o p e = N · m = 1 N m · X m m = 1 N m · m = 1 N X m m = 1 N X m · m = 1 N m 2 m = 1 N m 2
The spectral smoothness shows the variation of the frequency power with respect to the direct neighbours of the main frequency (Formula (24)).
smooth:
s m t = 20 · n = 2 N 1 l o g x n l o g x n 1 + l o g x n + l o g x n + 1 3
The spectral spread shows the variance of the main frequency energy around the mean value of the total frequency energy (Formula (25)).
spread:
s s p = m = 1 N m C s 2 · X m m = 1 N X m
The spectral summFS shows the energy of a frequency band defined by two boundaries, described by the area under the signal curve by summing its values (Formula (26)).
summFS:
s u m m F s = n = 1 N x n 2

2.7. Feature Selection and Regression Method

In a first step, all ultrasonic features presented previously were generated for each ultrasonic signal (see Figure 2A). To enhance the informative value of each signal, the features were taken for different time-areas (impedance, Echo I, Echo II). The time-areas therefore can be divided into echo parts and an impedance part (see Figure 2B). The spectral features were computed for the whole ultrasonic signal. Additionally, the stirring time and the temperature measured were used as features too. To avoid overfitting, the importance ranking of the predictors was applied using the nearest-neighbour method [41] based on the ten nearest neighbours. For further analysis, only the ten features (above the red line) with the highest importance weighting were kept (see Figure 2C).
To regress the features generated to the density of the batter, a Gaussian process regression (GPR) was performed [42]. The Gaussian process regression is a non-parametric kernel-based probabilistic model, which permits modelling data with varying multiple scales. The rational GPR model used can be described with Formula (27).
k x i , x j θ = σ f 2 ( 1 + r 2 2 α σ i 2 ) where : r = x i x j T ( x i x j )
θ : maximum of posteriori estimates, σ f : signal standard deviation, α: non-negative parameter of the covariance.
To select and train the GPR model, 70% of the measured data were used; 30% of the data were used for validation.

3. Results and Discussion

3.1. Density Measurement for Standard Batter without Variation

A first experiment was performed, in order to test the ultrasonic sensor with differently aerated standard biscuit batter. All tests were performed in three-fold determination. The time profile of density, as a function of whipping time for standard biscuit batter (see Table 1), is shown in Figure 3A. During mixing at high speed, the density of the standard batter drops quickly and remains steady after about 100 s. Longer aeration (more than 180 s) showed a negligible change on the extended aeration and density loss. Due to the foam-stabilizing effect of emulsifiers, only the foam structure changes [43] after 180 s of whipping time, whereas the density remains constant. The changes in the ultrasonic signal between 0 s and 180 s is shown in Figure 3B. The signal deflection in the signal from the ultrasonic sensor in the range between 0 and 10 µs is due to ringing of the sensor; and it is, therefore, it is not conducive to the measurement of material properties. Changes in the signal can mainly be detected because of the impedance signal (reflection PMMA/batter, 10–30 µs), although barely for the first echo (30–50 µs). The second echo (>50 µs) can only be detected with poorly aerated batters. Due to the strong attenuation and scattering effects of aerated batter for ultrasonic waves, the sound waves probably penetrate the biscuit batter only marginally and, ideally, are reflected back in such a way that these signals reach the ultrasonic receiver again (see Figure 1A). Therefore, the features generated from the echo areas play a minor role (see Figure 2C). At low densities, due to the increase in air bubbles in the batter, the ultrasonic signal is more scattered and therefore dampened. The prevailing longitudinal waves of the ultrasonic transducer are thus reflected at each interface (batter/air) and converted into other waveforms (surface waves, shear waves, etc.). By filtering the ultrasonic signals at the longitudinal frequency (1.80–2.25 MHz), the signal areas of the ultrasonic signals converted disappear from the signals shown. This is due to other frequency ranges and so they can only be detected within a certain frequency range. Therefore, the frequency-based features are also accorded greater importance (see Figure 2C). Utilizing the resulting features for the GPR model (see Formula (27)), a correlation of the density can be realized (see Figure 4). Since the measurement mainly takes place at the boundary surface of the buffer rod via the ultrasonic sensor, sticking batter residues or unfavourable bubble adhesion onto this surface can cause measurement outliers (see Figure 4). Except for these rare erroneous measurements, the density of standard biscuit batter can be measured in total up to a NRMSE (Normalized Root Mean Squared Error) of 6.32% (see Table 2).

3.2. Density Measurement for Batter Variations

A second experiment was performed in order to test the ultrasonic sensor with differently aerated standard and gluten-free biscuit batter. Furthermore, variations without an emulsifier were tested, as this has a major impact on the ability of biscuit batters to retain air bubbles. The time profile of density as a function of whipping time for all of the batter variations tested (see Table 1) is shown in Figure 5. The experiment series without an emulsifier could hardly be foamed within the period examined of 180 s. Only with significantly longer whipping times (>10 min), was it possible for the densities of the samples containing an emulsifier to be achieved approximately. The experiments also showed that the required density of approximately 600 g/L for highly aerated batters can be achieved, for both the standard batter and the gluten-free batter, by foaming with emulsifier. The whipping time required to reach this density hardly differs between the standard and gluten-free batter. It turns out that the gluten-free variant can produce lower overall densities. However, this only applies with the use of emulsifiers. Without emulsifiers, the standard biscuit batter whips up faster and can reach lower densities in the same time. This is because gluten itself can serve as an emulsifier, which would stabilize the structure in a standard biscuit. With regard to the features used, this means that the external feature of “whipping time” plays a less significant role when viewed overall. Due to the mode conversion of the ultrasonic signal at every boundary phase, the frequency spectrum of the ultrasonic signal showed changes in the density of the batter. Therefore, seven out of the ten features selected are spectral features (see Figure 2C). When using the same features from the previous experiment for modelling purposes, better results can be achieved overall with the extended dataset (see Figure 6). Outliers can also occur in these series of measurements, due to sticking batter residues or unfavourable bubble adhesion. However, these outliers occur preferentially in gluten-free samples for highly aerated batters due to the lack of stabilization by the gluten networks. The density prediction of the model for standard batter with emulsifier (R2 = 0.99, NRMSE = 1.31%) and without emulsifier (R2 = 0.99, NRMSE = 1.55%) is therefore slightly better than for gluten-free batter with (R2 = 0.96, NRMSE = 8.88%) and without (R2 = 0.94, NRMSE = 3.89%). A summary is presented in Table 3. In this work, only two varying biscuit masses (gluten containing, gluten free) without oil and fats were investigated. It was shown that the method also works for this type of variation in the foam system. However, for the determination of the density of other masses (oil containing, etc.) with varying rheological properties, the ultrasonic-based measurement system must be further adapted, since the surface properties change considerably. While in earlier publications the ultrasound-based measurement of density in grain-based foams was only presented as a potential (R2 < 0.7) [9,31], the method developed here can be used to make a significantly better prediction (R2 > 0.9) of the density, even for highly foamed batters.

4. Conclusions

In this paper, the design of an ultrasonic sensor for the mode-conversion measurement via ultrasonic features of highly aerated batters has been described. The ultrasonic sensor system developed was realized as a non-invasive system and, therefore, does not affect the process. The variation of bubbles in the batter due to alternating mixing times changes the density of the batter. In this way, by measuring the mode conversion with respect to the acoustic features of the batter, its density can be measured down to >500 g/L. Although the results already yield a reasonably prediction for highly aerated batters, (NRMSE = 6.32%), including for variated batter formulations (NRMSE = 4.15%), the method can still be extended and optimized.
Future research will extend the pulse–echo approach with multiple sensors for surface acoustic wave measurements in pitch–catch mode. This is in order to improve the accuracy and robustness and to increase the data from the signals due to the ultrasonic mode conversion. In this study, the final ultrasonic model was adjusted to the entire measurement range, as well as to the batter variations. To enhance this model, a restriction to one variation and a specific density range would be beneficial. In addition, for the machine-learning approach used, the rational GPR model was selected, by way of example, in order to provide the proof of concept. By comparing different machine-learning models (support vector regression (SVR), partial least-square regression (PLSR), neural network etc.) and by optimizing them [44], the approach shown in this study could be improved further. Another research direction would be to use the methodology developed in this study to measure further structural parameters (e.g., void fraction) of cereal-based foams [43] or to extend the method for further dough matrices. This would offer [43] a deeper insight into inline measurement systems.

Author Contributions

Conceptualization, M.M. and D.G.; methodology, M.M.; validation, M.M.; formal analysis, M.M. and D.G.; writing—original draft preparation, M.M.; writing—review and editing, M.M., D.G. and T.B.; supervision, D.G. and T.B. All authors have read and agreed to the published version of the manuscript.

Funding

This IGF Project “Optimization of Distributions Concerning Dimensions of Pores in Cereal-Based Foams Used for Production of Fine Pastry” of the FEI was supported via AiF (grant number AiF 18238 N) within the programme for promoting the Industrial Collective Research (IGF) of the Federal Ministry of Economic Affairs and Climate Action (BMWK), based on a resolution of the German Parliament. This research is part of the “REIF” project, supported via DLR within the innovation competition “Artificial Intelligence as a Driver for Economically Relevant Ecosystems” and is funded by the German Ministry of Economics and Climate Action (BMWK), based on a resolution of the German Parliament under the grant number 01MK20009Q. The Chair of Brewing and Beverage Technology thanks its industrial partners for supporting this research. Without the industrial data and the concomitant efforts expended by our partners, this research would not have been possible.

Data Availability Statement

The datasets generated for this study are available on request to the corresponding author, subject to demonstration of a legitimate interest.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Færgestad, E.M.; Magnus, E.M.; Sahlström, S.; Næs, T. Influence of Flour Quality and Baking Process on Hearth Bread Characteristics Made Using Gentle Mixing. J. Cereal Sci. 1999, 30, 61–70. [Google Scholar] [CrossRef]
  2. Hruskova, M.; Machova, D. Changes of wheat flour properties during short term storage. Czech J. Food Sci.—UZPI 2002, 20, 125–130. [Google Scholar] [CrossRef]
  3. Psimouli, V.; Oreopoulou, V. The effect of alternative sweeteners on batter rheology and cake properties. J. Sci. Food Agric. 2012, 92, 99–105. [Google Scholar] [CrossRef]
  4. Alifakı, Y.Ö.; Şakıyan, Ö. Dielectric properties, optimum formulation and microwave baking conditions of chickpea cakes. J. Food Sci. Technol. 2017, 54, 944–953. [Google Scholar] [CrossRef] [PubMed]
  5. Yildiz, E.; Guner, S.; Sumnu, G.; Sahin, S.; Oztop, M.H. Monitoring the Effects of Ingredients and Baking Methods on Quality of Gluten-Free Cakes by Time-Domain (TD) NMR Relaxometry. Food Bioprocess Technol. 2018, 11, 1923–1933. [Google Scholar] [CrossRef]
  6. Conforti, F.D. Cake Manufacture. In Bakery Products Science and Technology; Wiley: Hoboken, NJ, USA, 2014; pp. 563–584. [Google Scholar]
  7. Wilderjans, E.; Luyts, A.; Brijs, K.; Delcour, J.A. Ingredient functionality in batter type cake making. Trends Food Sci. Technol. 2013, 30, 6–15. [Google Scholar] [CrossRef]
  8. Campbell, G. Creation and characterisation of aerated food products. Trends Food Sci. Technol. 1999, 10, 283–296. [Google Scholar] [CrossRef]
  9. Salazar, J.; Turó, A.; Chávez, J.A.; García, M.J. Ultrasonic inspection of batters for on-line process monitoring. Ultrasonics 2004, 42, 155–159. [Google Scholar] [CrossRef]
  10. Bowler, A.L.; Bakalis, S.; Watson, N.J. A review of in-line and on-line measurement techniques to monitor industrial mixing processes. Chem. Eng. Res. Des. 2020, 153, 463–495. [Google Scholar] [CrossRef]
  11. Layton, R.A.; Murray, W.R.; Garbini, J.L. The feasibility of controlling power for efficient batch mixing. Mechatronics 1996, 6, 21–34. [Google Scholar] [CrossRef]
  12. Layton, R.A.; Murray, W.R.; Garbini, J.L. The Control of Power for Efficient Batch Mixing. Propellants Explos. Pyrotech. 1997, 22, 269–278. [Google Scholar] [CrossRef]
  13. Kilborn, R.H.; Preston, K.R. A dough height tracker and its potential application to the study of dough characteristics. Cereal Chem. 1981, 58, 198–201. [Google Scholar]
  14. Trinh, L.; Lowe, T.; Campbell, G.M.; Withers, P.J.; Martin, P.J. Bread dough aeration dynamics during pressure step-change mixing: Studies by X-ray tomography, dough density and population balance modelling. Chem. Eng. Sci. 2013, 101, 470–477. [Google Scholar] [CrossRef]
  15. Awad, T.S.; Moharram, H.A.; Shaltout, O.E.; Asker, D.; Youssef, M.M. Applications of ultrasound in analysis, processing and quality control of food: A review. Food Res. Int. 2012, 48, 410–427. [Google Scholar] [CrossRef]
  16. Chandrapala, J.; Oliver, C.; Kentish, S.; Ashokkumar, M. Ultrasonics in food processing—Food quality assurance and food safety. Trends Food Sci. Technol. 2012, 26, 88–98. [Google Scholar] [CrossRef]
  17. Khairi, M.T.M.; Ibrahim, S.; Yunus, M.A.M.; Faramarzi, M. Contact and non-contact ultrasonic measurement in the food industry: A review. Meas. Sci. Technol. 2015, 27, 12001. [Google Scholar] [CrossRef]
  18. Sarkar, T.; Salauddin, M.; Kirtonia, K.; Pati, S.; Rebezov, M.; Khayrullin, M.; Panasenko, S.; Tretyak, L.; Temerbayeva, M.; Kapustina, N.; et al. A Review on the Commonly Used Methods for Analysis of Physical Properties of Food Materials. Appl. Sci. 2022, 12, 2004. [Google Scholar] [CrossRef]
  19. Salimi Khorshidi, A.; Thandapilly, S.J.; Ames, N. Application of low-intensity ultrasound as a rapid, cost-effective tool to wheat screening: A systematic frequency selection. J. Cereal Sci. 2018, 82, 190–197. [Google Scholar] [CrossRef]
  20. Létang, C.; Piau, M.; Verdier, C.; Lefebvre, L. Characterization of wheat-flour–water doughs: A new method using ultrasound. Ultrasonics 2001, 39, 133–141. [Google Scholar] [CrossRef] [PubMed]
  21. Elmehdi, H.M.; Page, J.H.; Scanlon, M.G. Ultrasonic Investigation of the Effect of Mixing Under Reduced Pressure on the Mechanical Properties of Bread Dough. Cereal Chem. J. 2004, 81, 504–510. [Google Scholar] [CrossRef]
  22. Ross, K.A.; Pyrak-Nolte, L.J.; Campanella, O.H. The use of ultrasound and shear oscillatory tests to characterize the effect of mixing time on the rheological properties of dough. Food Res. Int. 2004, 37, 567–577. [Google Scholar] [CrossRef]
  23. Elmehdi, H.M.; Page, J.H.; Scanlon, M.G. Monitoring Dough Fermentation Using Acoustic Waves. Food Bioprod. Process. 2003, 81, 217–223. [Google Scholar] [CrossRef]
  24. Skaf, A.; Nassar, G.; Lefebvre, F.; Nongaillard, B. A new acoustic technique to monitor bread dough during the fermentation phase. J. Food Eng. 2009, 93, 365–378. [Google Scholar] [CrossRef]
  25. Bowler, A.L.; Bakalis, S.; Watson, N.J. Monitoring Mixing Processes Using Ultrasonic Sensors and Machine Learning. Sensors 2020, 20, 1813. [Google Scholar] [CrossRef]
  26. Scanlon, M.G.; Page, J.H. Probing the Properties of Dough with Low-Intensity Ultrasound. Cereal Chem. J. 2015, 92, 121–133. [Google Scholar] [CrossRef]
  27. Strybulevych, A.; Leroy, V.; Shum, A.L.; Koksel, H.F.; Scanlon, M.G.; Page, J.H. Use of an ultrasonic reflectance technique to examine bubble size changes in dough. In IOP Conference Series: Materials Science and Engineering, Proceedings of the International Symposium on Ultrasound in the Control of Industrial Processes (UCIP 2012), Madrid, Spain, 18–20 April 2012; Volume 42, p. 12037. [CrossRef]
  28. Allais, I.; Edoura-Gaena, R.-B.; Dufour, É. Characterisation of lady finger batters and biscuits by fluorescence spectroscopy—Relation with density, color and texture. J. Food Eng. 2006, 77, 896–909. [Google Scholar] [CrossRef]
  29. Fox, P.D.; Smith, P.P.; Sahi, S.S. Buffer rod design for measurement of specific gravity in the processing of industrial food batters. In 2002 IEEE Ultrasonics Symposium, Proceedings of the An International Symposium, Forum Hotel, Munich, Germany, 8–11 October 2002; Yuhas, D.E., Schneider, S.C., Eds.; IEEE: Piscataway, NJ, USA, 2002; pp. 679–682. ISBN 0-7803-7582-3. [Google Scholar]
  30. Fox, P.; Smith, P.P.; Sahi, S. Ultrasound measurements to monitor the specific gravity of food batters. J. Food Eng. 2004, 65, 317–324. [Google Scholar] [CrossRef]
  31. Gómez, M.; Oliete, B.; García-Álvarez, J.; Ronda, F.; Salazar, J. Characterization of cake batters by ultrasound measurements. J. Food Eng. 2008, 89, 408–413. [Google Scholar] [CrossRef]
  32. Colombi, A.; Ageeva, V.; Smith, R.J.; Clare, A.; Patel, R.; Clark, M.; Colquitt, D.; Roux, P.; Guenneau, S.; Craster, R.V. Enhanced sensing and conversion of ultrasonic Rayleigh waves by elastic metasurfaces. Sci. Rep. 2017, 7, 6750. [Google Scholar] [CrossRef]
  33. Mayer, W.G. Energy partition of ultrasonic waves at flat boundaries. Ultrasonics 1965, 3, 62–68. [Google Scholar] [CrossRef]
  34. Chaplain, G.J.; de Ponti, J.M.; Colombi, A.; Fuentes-Dominguez, R.; Dryburg, P.; Pieris, D.; Smith, R.J.; Clare, A.; Clark, M.; Craster, R.V. Tailored elastic surface to body wave Umklapp conversion. Nat. Commun. 2020, 11, 3267. [Google Scholar] [CrossRef] [PubMed]
  35. Clement, G.T.; White, P.J.; Hynynen, K. Enhanced ultrasound transmission through the human skull using shear mode conversion. J. Acoust. Soc. Am. 2004, 115, 1356–1364. [Google Scholar] [CrossRef] [PubMed]
  36. Wall, K. Complexity of chemical products, plants, processes and control systems. Chem. Eng. Res. Des. 2009, 87, 1430–1437. [Google Scholar] [CrossRef]
  37. Gao, W.; Liu, W.; Hu, Y.; Wang, J. Study of Ultrasonic Near-Field Region in Ultrasonic Liquid-Level Monitoring System. Micromachines 2020, 11, 763. [Google Scholar] [CrossRef]
  38. Amirkhani, M.; Taschin, A.; Cucini, R.; Bartolini, P.; Leporini, D.; Torre, R. Polymer thermal and acoustic properties using heterodyne detected transient grating technique. J. Polym. Sci. Part B Polym. Phys. 2011, 49, 685–690. [Google Scholar] [CrossRef]
  39. Kanatov, I.; Kaplun, D.; Butusov, D.; Gulvanskii, V.; Sinitca, A. One Technique to Enhance the Resolution of Discrete Fourier Transform. Electronics 2019, 8, 330. [Google Scholar] [CrossRef]
  40. Alías, F.; Socoró, J.; Sevillano, X. A Review of Physical and Perceptual Feature Extraction Techniques for Speech, Music and Environmental Sounds. Appl. Sci. 2016, 6, 143. [Google Scholar] [CrossRef]
  41. Robnik-Šikonja, M.; Kononenko, I. Theoretical and Empirical Analysis of ReliefF and RReliefF. Mach. Learn. 2003, 53, 23–69. [Google Scholar] [CrossRef]
  42. Zhang, N.; Xiong, J.; Zhong, J.; Leatham, K. Gaussian Process Regression Method for Classification for High-Dimensional Data with Limited Samples. In Proceedings of the 8th International Conference on Information Science and Technology, ICIST 2018, Granada/Cordoba/Seville, Spain, 30 June–6 July 2018; IEEE: Piscataway, NJ, USA, 2018; pp. 358–363, ISBN 978-1-5386-3782-1. [Google Scholar]
  43. Chesterton, A.; de Abreu, D.P.; Moggridge, G.D.; Sadd, P.A.; Wilson, D.I. Evolution of cake batter bubble structure and rheology during planetary mixing. Food Bioprod. Process. 2013, 91, 192–206. [Google Scholar] [CrossRef]
  44. Bowler, A.L.; Pound, M.P.; Watson, N.J. A review of ultrasonic sensing and machine learning methods to monitor industrial processes. Ultrasonics 2022, 124, 106776. [Google Scholar] [CrossRef]
Figure 3. (A) Time profile of density as a function of whipping time for the batch process of standard biscuit batter from 0 s (red) to 180 s (blue). (B) Ultrasonic signal change between 0 s (red) and 180 s (blue) whipping time.
Figure 3. (A) Time profile of density as a function of whipping time for the batch process of standard biscuit batter from 0 s (red) to 180 s (blue). (B) Ultrasonic signal change between 0 s (red) and 180 s (blue) whipping time.
Foods 12 01927 g003
Figure 4. Relationship between density predicted by ultrasonic signals and the reference (gravimetric determination) for standard I batter; model (dot (●), R² = 0.99, RMSE = 0.03), validation (circle (○), R² = 0.91, RMSE = 54.69).
Figure 4. Relationship between density predicted by ultrasonic signals and the reference (gravimetric determination) for standard I batter; model (dot (●), R² = 0.99, RMSE = 0.03), validation (circle (○), R² = 0.91, RMSE = 54.69).
Foods 12 01927 g004
Figure 5. Time profile of density as a function of whipping time for (A) standard I batter (black square) and standard II (square) without an emulsifier and (B) gluten-free I batter (black dot) and gluten-free II (dot) without an emulsifier.
Figure 5. Time profile of density as a function of whipping time for (A) standard I batter (black square) and standard II (square) without an emulsifier and (B) gluten-free I batter (black dot) and gluten-free II (dot) without an emulsifier.
Foods 12 01927 g005
Figure 6. Relationship between density predicted by ultrasonic signals and the reference (gravimetric determination) for all variations (circle: standard I, R² = 0.99, RMSE = 6.32; square: standard II, R² = 0.99, RMSE = 7.46; diamond: gluten-free I, R² = 0.96, RMSE = 42, triangle: gluten-free II, R² = 0.94, RMSE = 18.66).
Figure 6. Relationship between density predicted by ultrasonic signals and the reference (gravimetric determination) for all variations (circle: standard I, R² = 0.99, RMSE = 6.32; square: standard II, R² = 0.99, RMSE = 7.46; diamond: gluten-free I, R² = 0.96, RMSE = 42, triangle: gluten-free II, R² = 0.94, RMSE = 18.66).
Foods 12 01927 g006
Table 1. The composition of the different sponge cake batter variants.
Table 1. The composition of the different sponge cake batter variants.
IngredientIngredient Consumption g/5 kg
Standard IStandard IIGluten-Free IGluten-Free II
Whole egg1417141714171417
Wheat flour (Type 550)75075000
Wheat starch92597516751725
White sugar1858185818581858
Emulsifier500500
Table 2. Summary of the R2, RMSE (Root Mean Squared Error) and NRMSE (Normalized Root Mean Squared Error) of the model for standard batter.
Table 2. Summary of the R2, RMSE (Root Mean Squared Error) and NRMSE (Normalized Root Mean Squared Error) of the model for standard batter.
R2RMSENRMSE [%]
Model0.990.030.006
Validation0.9154.6911.41
Total0.9630.316.32
Table 3. Summary of the R2, RMSE (Root Mean Squared Error) and NRMSE (Normalized Root Mean Squared Error) of the models for batter variations.
Table 3. Summary of the R2, RMSE (Root Mean Squared Error) and NRMSE (Normalized Root Mean Squared Error) of the models for batter variations.
R2RMSENRMSE [%]
Standard I0.996.321.31
Standard II0.997.461.55
Gluten-free I0.9642.538.88
Gluten-free II0.9418.663.89
Total0.9819.924.15
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Metzenmacher, M.; Geier, D.; Becker, T. Ultrasonic Wave Mode-Based Application for Contactless Density Measurement of Highly Aerated Batters. Foods 2023, 12, 1927. https://doi.org/10.3390/foods12091927

AMA Style

Metzenmacher M, Geier D, Becker T. Ultrasonic Wave Mode-Based Application for Contactless Density Measurement of Highly Aerated Batters. Foods. 2023; 12(9):1927. https://doi.org/10.3390/foods12091927

Chicago/Turabian Style

Metzenmacher, Michael, Dominik Geier, and Thomas Becker. 2023. "Ultrasonic Wave Mode-Based Application for Contactless Density Measurement of Highly Aerated Batters" Foods 12, no. 9: 1927. https://doi.org/10.3390/foods12091927

APA Style

Metzenmacher, M., Geier, D., & Becker, T. (2023). Ultrasonic Wave Mode-Based Application for Contactless Density Measurement of Highly Aerated Batters. Foods, 12(9), 1927. https://doi.org/10.3390/foods12091927

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop