Start main page content

Study Group Problems

Study Group Problem 1 : The Batch Pan Scheduling Problem in a white sugar refinery

Industry:  Sugar

Industry Representative:  Dr Richard Loubser, Sugar Milling Research Institute, University of KwaZulu-Natal, Durban

Moderator: Professor Montaz Ali, University of the Witwatersrand, Johannesburg

Student Moderator:

Problem Statement

In a typical white sugar refinery, the sugar is crystalised in batch “Vacuum Pans”. These are evaporative crystallisers where evaporation and crystallisation takes place simultaneously under vacuum (to be able to operate at a reduced pressure where the boiling temperature is reduced and there is thus less thermal degradation of sucrose). At the end of each batch pan cycle the pan will contain a mixture of crystals and mother liquor with 50% of the sugar (typically) being in crystal form.

The crystallisation is followed by centrifugation to separate the crystal sugar from the mother liquor. There is some water addition at this stage to wash the sugar and dissolve small crystals that can pass into the separated mother liquor (normally called jet). The crystallisation is done in multiple stages with the impurities, predominantly “colour”, remaining in the mother liquor.

The flow scheme for a typical refinery multiple stage crystallisation process is shown in the diagram below:

The pieces of process equipment used in this flow scheme are:

The crystallisation takes place in a batch pan where water is removed from the liquor fed to the pan, using heating steam to drive the evaporation. When each batch is complete, the contents of the batch pan are discharged into a “strike receiver” which will have a capacity of 1.5 times the size of the batch. This means that the “strike receiver” does not have to be completely empty before a batch pan is discharged into it.

The mixture of crystals and mother liquor (called massecuite) will be processed in centrifugals to separate the crystals. Although this is a batch process, it can be considered as continuous for this exercise as the centrifugal batch sizes are small relative to the pan batch sizes.

The impurities concentrate up across each crystallisation cycle. Of particular interest is “colour” – an analytical result that can be treated as if it is concentration (expressed as parts per million on solids).

The colour in the sugar is dependant on the colour remaining in the mother liquor. With 50% of the sucrose being removed from the liquor in each stage of crystallisation the mother liquor colour and the sugar colour doubles for each stage. At the same time the quantity of sugar produced in each stage halves. Typical sugar colours would be :

  • First sugar : 20
  • Second sugar: 40
  • Third sugar: 80
  • Fourth sugar: 160.

These four sugars will be mixed together to give an average final sugar colour of :

               (1 * 20) + (0.5 * 40 ) + (0.25 * 80) + (0.125 * 160) / (1 + 0.5 + 0.25 + 0.125) =  70.7

It is important that the rate of production of sugars is maintained in this ratio so that the average sugar colour (of the mixed sugar) does not deviate from the long term average value.

The length of the batch cycles becomes longer for the later crystallisation stages as a result of the concentrating impurities slowing down the crystallisation rate.

In terms of steam demand, a typical batch pan cycle can be approximated as follows.

Typical values for a batch pan of a standard design that is able to produce 75 tons of massecuite in each batch pan cycle would be as follows:

The following performance data for pans of this design can be assumed:

As a first approximation, assume that the demand for liquor to a pan is spread evenly over the batch cycle.

The task is to develop a scheduling program that allocates pans appropriately to be able to process a specified quantity of fine liquor (say 100 tons/hr with a concentration of 76% - thus 76 tons of dissolved sugar) with the following constraints:

  • Achieve the smoothest possible total steam demand from all pans
  • Add waiting times between batches to achieve the required scheduling
  • Do not start a batch pan unless there is sufficient liquor/jet in the feed tank to complete the batch
  • There must be sufficient space in the strike receiver to accommodate the contents of a batch pan at the end of its cycle
  • Processing of sugar through the centrifugals must match the requirements for correct proportional mixing of sugars of different grades (so as to maintain the required average sugar colour)
  • Determine the minimum number of pans necessary to meet the requirements.

 

References

 

Supporting Material

 

Presentation

MISG 2024 Study Group Problem 1 Presentation - Scheduling

Report-back Presentation

MISG 2024 Study Group Problem 1 Report-back presentation Batch_Pan_Scheduling

Study Group Problem 2: Rogue wave in the Agulhas region

Industry:  Shipping

Industry Representative:  Mr. Morwakoma M Matabane, Research and   Development Scientist, South African Weather Service

Moderator: Dr Thama Duba, School of Computer Science and Applied Mathematics, University of the Witwatersrand, Johannesburg

Student Moderator:

 

Problem Statement

Ocean waves in the Southern African region are caused by the interplay of wind, current, bathymetry, tides, etc. Recently South Africa has been experiencing a high frequency of rogue/freak waves in the Agulhas region [or Southern Eastern Coastlines, etc.] that has resulted in loss of life, livelihood, flooding, and damage to infrastructure. The problem is centered around finding out the contribution of the wind-current, wind-bathymetry, and wind-tide interplay towards making waves extreme/severe in the Southeastern parts of Southern Africa.

References

MISG 2024 Study Group Problem 2 Rogue Waves Reference - Didenkulova et al 2023

MISG 2024 Study Group Problem 2 Rogue Waves Reference - Fedele et al 2016

Supporting Material

MISG 2024 Study Group Problem 2 Rogue Waves Supporting Material - MMMatabane

Presentation

MISG 2024 Study Group Problem 2 Presentation - Rogue wave

Report-back Presentation

MISG 2024 Study Group Problem 2 Report back presentation - Rogue Waves

Study Group Problem 3: Modelling a gas chromatograph

Industry:   Agriculture, food and beverage, pharmaceutical

Industry Representative:    Dr Alba Cabrera-Codony, Laboratory of Chemical and Environmental  Engineering, University of Girona, Spain

Moderator:    Professor Tim Myers, Centre de Recerca Matemática, Barcelona, Spain

Student Moderator:

Problem Statement

Gas chromatography (GC) is a widely used technique for separating and analysing volatile compounds in gas mixtures. GC has a broad range of applications, for example in detecting and quantifying pollutants, pesticides and environmental contaminants in air, water and soil samples. In food and beverage analysis, it can be used to determine the presence and concentrations of flavour compounds or additives. In the pharmaceutical industry, it permits the analysis of drugs, including purity and quantifications of active ingredients. In forensic analysis it plays a crucial role in toxicology and arson investigations. Additionally, in clinical and medical laboratories, GC is used for analysing blood, urine and other biological samples.

Here's a brief overview of how gas chromatography works:

  1. A small amount of the sample containing different compounds is injected into the chromatograph.
  2. A carrier gas, usually helium or nitrogen, carries the vaporized sample through the chromatographic column. The column is a long, thin tube: there are two types of column, i) Packed Columns, which contain a solid support material coated with a liquid stationary phase, and ii) Capillary Columns, which are narrower and contain a thin film of stationary phase on the inner wall.
  3. Adsorption: This causes separation, based on the different interactions of the compounds, with the stationary phase.
  4. Separation: As the sample travels through the column, different compounds interact differently with the stationary phase. The extent of interaction depends on factors such as molecular size, polarity, and volatility. Compounds that interact strongly with the stationary phase will spend more time in the column, leading to a slower movement through the column. Conversely, compounds that interact weakly will move faster.
  5. Detection: As each separated compound exits the column, it passes through a detector. The detector responds to the presence of compounds and generates signals proportional to the concentration. Common detectors include flame ionization detectors (FID), electron capture detectors (ECD), and mass spectrometers (MS).
  6. The signals from the detector are recorded and analysed to produce a chromatogram, which is a graphical representation of the separation of compounds over time that permits the data analysis.

Figure 1: Schematic of a gas chromatography system: a sample is added to the carrier gas flow. This passes through the column and concentrations are measured at the outlet and the data passed to a computer.

Adsorption and desorption are fundamental processes that occur in GC, contributing to the separation of compounds as they move through the chromatographic column. Adsorption refers to the adhesion of molecules to the surface of a solid. In gas chromatography, the stationary phase, which is typically a solid support or a liquid film on the column wall, acts as the adsorbent. The compounds in the vaporized sample interact with the stationary phase and the strength of this interaction depends on the chemical nature of both the stationary phase and the sample compounds.

Compounds that strongly adsorb to the stationary phase spend more time in the column, resulting in longer retention times. The degree of adsorption is influenced by factors such as the temperature, type of stationary phase, and the nature of the sample compounds.

Desorption is the process of releasing adsorbed molecules from the surface of the stationary phase. As the carrier gas (mobile phase) transports the sample through the column, the temperature of the column is typically controlled.

GC analyses use oven temperature programming, where the column temperature is increased gradually during the analysis. Increasing temperature speeds up the desorption of compounds by progressively reducing their affinity for the stationary phase. This allows for efficient separation of different compounds based on their boiling points and affinities for the stationary phase.

The final outcome is the chromatogram, a graphical representation illustrating the compounds analysed by the detector over time. In Figure 2 we present a typical example.

Figure 2: The chromatogram demonstrates the separation of 4 volatile organic compounds using a Flame Ionization Detector (FID): The x-axis represents time. Each peak corresponds to a specific compound and is related its retention time, i.e. the time it takes to travel through the chromatographic column and reach the detector. The y-axis represents the signal obtained with the FID. The chromatogram shows how the signal intensity changes over the course of the analysis, providing information about the concentration of the sample.

References

 

Supporting Material

 

Presentation

 MISG 2024 Study Group Problem 3 Presentation - Gas chromatograph

Report-back Presentation

 

Study Group Problem 4: Axial strain evaluation without the use of strain gauges

Industry:  Mining

Industry Representative:  Dr Halil Yilmaz, Rock Mechanics Laboratory, CC CSIR,  Melville, Johannesburg

Moderator:

Student Moderator:

Problem statement

Uniaxial compressive strain testing with the use of axial and lateral strain gauges is one of the most common tests required by the mining industry. Strain gauges are used for the calculation of the Young’s modulus and the Poisson ratio. Young’s modulus is generated  from the axial stress versus actual strain data and the Poisson ratio is calculated using the lateral strain versus axial strain data. Most universal testing machines are capable of measuring the force and the machine displacement as a crosshead movement in the direction of the applied force.  The machines unfortunately cannot measure the deformation of the specimen in the same direction.

Is it possible to infer the axial deformation taking place on the rock specimen without using a  strain gauge by using  the data generated by the machine.

A strain gauge can be used only once and is expensive. The materials remain elastic throughout the deformation and linear elasticity described by the generalised Hooke’s law applies.

Shown  below are:

  • A universal testing machine with a rock specimen
  • A rock specimen before and after testing with axial and lateral strain gauges
  • Graphs of stress plotted against strain and load plotted against machine crosshead movement

Universal testing machine with rock specimen

A rock specimen before and after testing with axial and lateral (radial) strain gauges.

Example: Stress vs Strain curves

Example: Load vs Machine Crosshead Movement

 

References

 

Supporting Material

 

Presentation

MISG 2024 Study Group Problem 4 Presentation - Axial Strain

Report-back Presentation

MISG 2024 Study Group Problem 4 Report-back Presentation - Axial strain

Study Group Problem 5: Tourist attractions capping visitor numbers

Industry: Tourism Sector

Industry Representative: Dr Lombuso Precious Shabalala, University of South Africa 

Moderator: Dr M Sejeso, Wits University and Dr S Simelane, University of Johannesburg

Problem Statement:

Case study of Manyeleti Nature Reserve, Mariepskop Nature Reserve, Bushbuckridge N

The continuing tourism growth will eventually result in increased visitation to some destinations / tourist attractions. Certain tourist attractions are limiting the number of visitors they welcome each day. The main reasons are to protect sensitive environments and provide a more enjoyable visitor experience by lessening the crowd.  It must be noted that tourism caps have been around for decades and the pandemic travel patterns have encouraged new restrictions to take effect.   It is perceived that visitation caps make for an inherently less flexible travel experience but as well as minimizing crowd and putting less strain on staff (which continue to be in short supply), limiting the number of visitors also helps to preserve and conserve natural resources.

It is vital for natural attractions to sustain the physical or ecological impact of visitors. The issue for managers surrounds the number of visitors that can be accommodated before the experience provided by the attraction is compromised. This challenge can be resolved through determining the attraction's social carrying capacity (SCC) taking social comfort level (SCL) into account and traffic flow management.

However, the challenge of managing tourism sustainably for residents, tourists and day visitors has been recognised.  Hence, there is a radical change in the perceptions of local people to tourism, and in many destinations a tipping point has been reached and mass tourism has become a local political issue. Mass-tourism carries the same desription  as Overtourism. Goodwin (2017)  views overtourism as a new term that describes “ destinations where hosts or guests, locals or visitors, feel that there are too many visitors and that the quality of life in the area or the quality of the experience has deteriorated unacceptably” .  Masstourism and Overtorism are  the opposite of Responsible Tourism which is about using tourism to make better places to live in and better places to visit.  The  lack of responsible tourism practice often results in visitors and guests experiencing the deterioration concurrently and rebel against it (Goodwin, 2017).

 Literature suggest that, there are several tourist attractions in the country (Mpumalanga Province) that are deserted in a way. These tourist attractions have the potential to realise their sustainable development while practicing responsible tourism in their daily operation activities and generating healthy revenue. Tourism development in South Africa is guided by the key principles of Responsible Tourism stipulated in the 1996 White Paper.  South Africa was the first country to include Responsible Tourism in its national tourism policy, as outlined in the 1996 White Paper on the Development and Promotion of Tourism in South Africa. Goodwin (2007) noted that Responsible Tourism is about "making better places for people to live in and better places for people to visit" Responsible Tourism requires operators, hoteliers, governments, local people, and tourists take responsibility and action to make tourism more sustainable.  

 

The problem to investigate:

The objective of the project is to develop a viable mathematical model to determine the social carrying capacity of a tourist attraction to mitigate the negative impact of over-tourism while providing a high-quality experience for visitors. The model should consider the available infrastructure, activities, natural and cultural resources, and accommodation. The tourist attractions in consideration are Manyeleti Nature Reserve, Mariepskop Nature Reserve, Bushburckrigde Nature Reserve, and Injaka Dam.  

The solution to the problem will be of global importance to the tourism sector.

 The Study Group is asked to develop further a simple model  proposed at MISG 2023 to determine the carrying capacity by simulating the visitor flow of a given tourist attraction. The model is derived from traffic flow problems and principles of advection-reaction-diffusion equations. In the case of visitor flow into the tourist attractions, a traffic light network is analogous to a network of different attractions (activities). The model  considers the following parameters for a particular tourist attraction: the arrival rate, length of stay, activities, and attraction capacity. The aim is to avoid long queues for any activity (attraction) in the network. The Study Group is asked to use simulations to determine the carrying capacity of each tourist attraction, which is the  point at which long queues are unavoidable in a network. The model is straightforward, easy to understand, easy to implement, and immensely useful. It will require data to simulate real-world situations and identify key model parameters. The  Study Group is asked to  validate the  model by fitting it to real-world situations.

 

References

Goodwin, H., 2017. The challenge of overtourism. Responsible tourism partnership, 4, pp.1-19.

Goodwin, H., 2007.Taking responsibility for Tourism. Accessed date: 21 November 2022. Available from: https://haroldgoodwin.info/responsible-tourism/

Morgan, D. and Lok, L., 1999. Social Comfort Within Natural Tourist Attractions: A Case Study of Visitors to Hanging Rock, Victoria.

Supporting Material

MISG2024 Study Group Problem 5 Supporting Material

Presentation

MISG 2024 Study Group Problem 5 Presentation - Tourism

Report-back Presentation

Study Group Problem 6: Acoustic Signal in Wind Turbine

Presenter: Adam Mielke

Problem Statement:

The overall aim is to investigate how to detect faults in wind turbine blades through acoustic emission (AE) measurements. In lab experiments, we move the blade up and down and detect the sound it makes through analogue sensors. The blade in this experiment is 14.3 m long. Four embedded faults are made during the construction of the blade and AE-sensors are placed around them, see Figures 2 and 3. The blade is then excited by an arm that moves it in two different directions at once: Flap-wise (bending along the narrow direction of the blade) and edgewise (bending along the wide direction of the blade). The total movement is perhaps best approximated by a figure eight. The frequency of the flap-wise movement matches an eigenfrequency of 2.305 Hz, which has a quite narrow band and therefore can be regarded as a constant. Note that this is of course a resonance that is avoided in the real world, because the large amount of movement puts a bigger stress on the blade. (If we didn't use this frequency in the lab, it would take forever to test the growth of faults.) As the blade is alternating between C-shaped and figure of 8 movements, the edgewise motion has twice the frequency of the flap-wise motion. The conversion from analogue to digital is a bottleneck, so the full signal is not saved. Instead, a number of waveform parameters are saved, see Table 1 for a list of these and Figure 1 for an illustration of them. The wave packets are strongly dampened, so hardly any acoustic events reach more than the two adjacent sensors. We use a value for the speed of sound along the blade of around 2 km/s, but unfortunately the actual speed of sound is anisotropic (and most likely non-uniform).

Read the full problem statement here: MISG 2024 Study Group Problem 6 - Acoustic signal in wind turbine

References

 

Supporting Material

 

Presentation

MISG 2024 Study Group Problem 6 Presentation - Wind Turbine

Report-back Presentation

MISG 2024 Study Group Problem 6 Report-back presentation - Acoustic Signal

Share