Next Article in Journal
Tourism Adaptation to Coastal Risks: A Socio-Spatial Analysis of the Magdalen Islands in Québec, Canada
Previous Article in Journal
A Gateway to Rapid Prediction of Water Quality: A Case Study in China’s South-to-North Water Diversion Project
Previous Article in Special Issue
Advanced Drinking Groundwater As Phytofiltration by the Hyperaccumulating Fern Pteris vittata
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Proposal of a Laboratory-Scale Anaerobic Biodigester for Introducing the Monitoring and Sensing Techniques, as a Potential Learning Tool in the Fields of Carbon Foot-Print Reduction and Climate Change Mitigation

by
Saulo Brito-Espino
1,
Federico Leon
1,*,
Jenifer Vaswani-Reboso
1,
Alejandro Ramos-Martin
2 and
Carlos Mendieta-Pino
1
1
Institute of Environmental Studies and Resources (IUNAT), Campus Tafira, University of Las Palmas de Gran Canaria, 35017 Las Palmas de Gran Canaria, Spain
2
Departament of Process Engineering, Campus Tafira, University of Las Palmas de Gran Canaria, 35017 Las Palmas de Gran Canaria, Spain
*
Author to whom correspondence should be addressed.
Water 2021, 13(17), 2409; https://doi.org/10.3390/w13172409
Submission received: 4 August 2021 / Revised: 27 August 2021 / Accepted: 29 August 2021 / Published: 1 September 2021
(This article belongs to the Special Issue Water Treatment: Desalination, Treatment, Reuse and Management)

Abstract

:
This article presents a proposal of an anaerobic biodigester on a laboratory scale for introducing the monitoring and sensing techniques of the growth of microorganisms according to different parameters, where the redox potential, pH, pressure, and temperature have been measured in quasi-continuous mode. For this task, a microcontroller system was used (Atmega328—Arduino). Importantly, the design is based on flexible and open-source software, hardware, and firmware (Scilab, Arduino, Processing), facilitating its modification for other related studies. This design was developed to help engineering students to learn and to understand the operation of an anaerobic biodigester, which allows us to know various properties of the system at any time, as well as its evolution over time. In this way, property curves can be drawn and related to each other to obtain a better understanding of the biodigester operation. In this context, the relationship between the oxide-reduction reaction and microbial activity was studied so that the redox potential can be a way of measuring the growth of microorganisms in an anaerobic environment. With all this, through these parameters, it is possible to introduce to engineering students the operation of this technology used normally like a very powerful tool for the control of the carbon footprint, for example in wastewater sector, and consequently for the mitigation of the climate change.

1. Introduction

Anaerobic digestion technologies, applied to organic water treatment, are efficient ways to solve environmental problems also provide energy. They are considered sustainable, safe, and efficient biotechnologies in which carbon footprint reduction, by CH4 capture and fossil fuel replacement, is clearly a factor to take under advisement [1,2,3,4,5,6,7]. The process of anaerobic digestion has been known and applied since ancient times; however, it was understood in terms of its final products and not its processes [3]. The versatility of anaerobic digestion applied as an effective technology in the face of certain fundamental challenges has found its usefulness in biotechnological industries [4,5,6,7]. Unlike aerobic processes where dissolved oxygen can be measured continuously, there is a great challenge for fermentative processes in anaerobic organisms where the technologies referring to control processes are currently insufficient [6,7,8,9,10]. Since pH detection has been commonly used in fermentation processes, where only the activity of the proton is reflected, it is not sensitive to small changes in the intracellular metabolism. The redox potential (ORP) known as oxidation-reduction or oxide-reduction potential, reflects all the electrons transferred and reflects the intracellular metabolism [6,11,12,13,14,15].
Recent advances in analytical technologies allow complex bioconversion processes to be controlled and deciphered. Few parameters in this process are routinely recorded continuously and immediately, such as pH, ORP, gas production rate, and flow rates [16,17,18,19].
Within this context, and due to the great number of applications that are being given to the anaerobic biodigesters, that is why it would be necessary to develop a strategy so that the engineering students could understand and develop its operation as well as the parameters that govern it, all applied to different situations. The learning of these strategies could be achieved using this equipment or through experimental designs carried out by students. This educational proposal is based on theoretical psychological studies published in numerous articles [5,7,8,20,21,22,23,24] which emphasize that students can reinforce their learning through an appropriate teaching environment, as well as through the use, construction, and design of equipment.
After describing the importance of the anaerobic process and the consequent need to control it, at all times, the decision was taken to present the design and manufacture of an anaerobic batch biodigester on a laboratory scale, as well as its implementation through a practical application in which brewer’s yeast (Saccharomyces cerevisiae) was introduced, and it was subjected to a semi-continuous control process, the results of which were subsequently compared with a previously proposed theoretical model [25,26,27,28,29,30].
The main objective of this article is to show an experimental design at laboratory scale of an anaerobic biodigester. In the same way, a series of tools, software, and hardware are proposed, which are easy to use and of low cost, and which will allow engineering students to see that they can develop autonomous elements to control an element, as well as to transform the information received into parameters that can later be interpreted in a computer [31,32,33,34]. The specific objective of this research is to relate the process of anaerobic digestion in a sequentially loaded reactor on a laboratory scale for a known microorganism and with a substrate prepared in the laboratory, to the profiles of oxidation-reduction potential, pH, temperature, and absolute pressure [35,36,37,38,39].

2. Materials and Methods

2.1. Diagram of the Laboratory Reactor

The designed bioreactor (Figure 1) can be grouped into three distinct parts: (1) Digestion system—which includes the digester itself, as well as those elements that are in direct contact with it (sensors, heating cable, loading, and unloading supply, etc.); (2) Control system, circuits, and voltage source—it receives data and sends the orders necessary for the proper functioning of the system; and (3) Computer system, communication interface and software [40,41,42,43,44].

2.1.1. Digestion System

It is made up of an insulated hermetic container, with a feeding and evacuation system arranged so that the mixture is guaranteed in each loading and unloading process. The upper part of the container has a series of sensors that are defined below:
  • pH Sensor: Scientific Grade Silver/Silver Chloride pH), 10 sensor with a response speed of 95% in one second.
  • ORP sensor, (E): High quality sensor from Atlas Scientific [10,45,46,47,48]. The data transmission mode is through an integrated system and with a simple serial communication protocol which gives us an immediate response of the E value.
  • Absolute Pressure Sensor: Phidgets mod. 1141-0—Absolute Sensor of gas pressure from 15 to 115 kPa [49,50]. This is a high-level sensor with analogue input, with input proportional to the of the environment. The pressure measurement for this sensor is 15 kPa. The formula used to translate the sensor value into pressure was the following [2,51,52]:
Pressure (kPa) = (Sensor Value/9.2) + 10.
  • Temperature sensor: Two miniature Vishay NTC thermistors (Figure 2) were used to take external and internal temperature readings. Their main characteristics are described in Table 1.
For the calculation of the temperature from the analogical measurement, considering the resistive values depending on the temperature provided by the manufacturer, and with an algorithm in Scilab, we obtain Equation (2).
f(x) = 2.249 × 10−5x2 + 0.06872x − 16.03

2.1.2. Circuits and Control System

Arduino Uno (Figure 3) microcontroller model ATmega328 (Atmel) was implemented within an embedded system in order to control the measuring processes providing bidirectional communication with the circuit of the electrical conductivity probe, transferring the respective data to the PC via USB for archival purposes.
Figure 4 and Figure 5 show the general circuit diagram and pictures where all the elements necessary for the correct operation of the system are collected, as well as the data collection to be processed later. It is composed of a voltage source of 12 V that feeds: the transistor, a heat source, thermistors and a stirring system, a microcontroller hard plate, a PWM plate through which the bioreactor temperature is controlled by a transistor, a plate for the temperature sensors, auxiliary connection plates and sensor plates, resistors, diodes, and wiring.

2.1.3. Computer System, Communication Interface, and Software

As for the computer system, a data transmission source code was developed for the pH, E, pressure, and temperature sensors with Processing software (source code in Appendix C). The output data was transferred to the Arduino ide serial monitor for checking and control, and then through the Processing software interface for saving into file the sensor samplings.

2.1.4. Auxiliary Equipment and Laboratory Material

The auxiliary equipment were the following: an Atago RX-7000 Alfa3 refractometer used for the effluent measurements, the refraction product, and the Brix degrees (°Bx); and precision weights. As far as materials are concerned, all those related to the sampling and measurement of volumes (typical of a laboratory) were used, such as flasks, pipettes, spoons, etc.

2.2. Preparation and Testing

2.2.1. Preparation of the Substrate

For the preparation of the substrate, the procedure described in bibliography was followed [4]. (1) The substrate was prepared, and its suitability checked. (2) The Brix level was measured and verified to be between 17 and 20 degrees. (3) Once this process was finished, 200 mL of the must was taken, and the yeast was added to it (approximately 2–4 g/L). Although an activation temperature of 37 °C is normally required, it was left at room temperature, as it has been previously proven that the inoculum is active under these conditions for working yeast.

2.2.2. Microbial Inoculum

Brewer’s yeast was used whose species includes Saccharomyces cerevisiae with a yield of 0.25–0.33 kg of dry cell weight per kg of substrate.

2.2.3. Control and Saving Data

Through the Arduino (Appendix B shows the Arduino source code), the temperature was controlled and the signals from the pH, ORP, absolute pressure, and temperature sensors (inside the digester and outside the environment) were read. Data were sent to the computer where they were stored by means of the use of Processing tool. Shown in Appendix C is the Processing source code, and Figure 6 displays the interface. Once all the information was entered, it was saved in a file on the computer’s hard disk.

3. Results

3.1. Digestion Model

Dynamic simulation between reality and model is a very important way for research, as it enables to provide strategies for the digester operation. In the model, the description of all metabolic rates is based on the classical Monod equation,
μ = μ m a x S K s + S
where μ is the growth rate of a considered biomass, μ m a x is the maximum growth rate of this microorganism, S is the concentration of the limiting substrate S for growth, K s is the half-velocity constant for the substrate S when μ μ m a x   =   0.5 .
The matrix differential equation related to biomass and substrate dynamic is as follows:
d X d t = ( μ m a x S K s + S ) X
d S d t = ( β s μ m a x S K s + S ) X
where S and X are the substrate and biomass concentration respectively, β s is the stoichiometric ratio for S.
Figure 7 shows a simulation with discontinuous dynamics (Scilab source code in Appendix A) with period T = 24 h. Both the evolution of the substrate and that of the biomass, as can be seen in the system, start from an initial state and, after the transitory process, reach a stationary state.

3.2. Anaerobic Digester Start-Up and Operation

An average representation of the tests, carried out in the biodigester over 5 weeks, is shown in Table 2. The data was processed using a computerised tool from Scilab. Scilab is a software for numerical analysis, with a high-level programming language for scientific calculation. With the obtained data, a series of graphs were elaborated and the most relevant ones are presented in Figure 8, Figure 9 and Figure 10.
The expression of the oxidation-reduction reactions can be expressed by the Nernst Equation (3).
E = E 0 + 2.303 R T n F l o g 10 ( Product   of   activities   of   oxidized   species   Product   of   reduced   species   activities )
where E0 is the standard ORP, n is the number of exchanged electrons, and F is the Faraday constant (96.42 kJ/g equivalent volts).
The three periods from Table 2 are represented in the graphs as they are the most illustrative. It had taken a time of 1–3000 min for each of them.

3.2.1. First Stage

This stage includes from the start-up of the bioreactor to the first charge/discharge process. Figure 8 shows the ORP, pH, and temperature versus the time, in minutes. The pH profile remains stable with adequate value for fermentative processes, around 4.4, with slight oscillations, while the ORP values indicate, practically, reducing conditions, in a slightly wider range that the previous one. This last achieves its maximum value, 100 (mV), at 1000 min, and minimum value, −250 mV at 1300 min; however, it tends to stabilize with fermentation time. This performance is reflected as well in other works [53,54]. On the other hand, during the whole of this period, temperature moved between 23.6 and 24.6 °C. It is observable that for high value of temperature, the ORP ’s graph tends to drop.
In general, the ORP graph shows a downward slope with a very irregular profile.

3.2.2. Second Stage

In this occasion, the system was fed back by 75 mL of a new mixing (see Table 2). From Figure 9, it is observed, in the first third of the stage, a sharp fall in ORP, −85 mV; it immediately increases until it reaches a maximum value, 40 mV, then it begins to decrease gradually. This may have occurred due to the supply of substrate indicating an increase in activity in the first third of the period recorded. The difference in feedstock could change the microbial community and dominant species in anaerobic digestion process. In relation to temperature, the ORP maintains the same vein as that in the previous case.

3.2.3. Third Stage

For this stage, the alkalinity in the reactor was increased by the adding NaOH mixed with the substrate in the follow feeding. The basic environment of the system reflects negative ORP values, between −300 and −470 mV. The graph, Figure 10, shows a gradual decreasing trend of ORP along this period. It is due to the buffering capacity of anaerobic digestion. Similar results have been achieved by other authors [55,56,57,58].

4. Conclusions

The proposed anaerobic biodigester monitored system works well, allowing for small-scale testing. Its use in research and teaching will allow the development of new research projects in the same way that it will help engineering students in their learning. With this design, it will be easy to determine the factors that can affect the growth of the anaerobic microorganism through the continuous data collection by the ORP, pH, temperature, and pressure sensors. Subsequently, with the subsequent processing of the data, it is possible to make graphs to be contrasted with a previously proposed theoretical model, and at the same time it can be compared with the equations that govern its behavior (Nernst equation).
On the other hand, this design has been supported by flexible and easily accessible free software, this being an important advantage for students since it offers the possibility of adapting this experimental design to each specific case, and the whole software used is free and open source.
The results show that the experimental design is feasible for the control and data collection of magnitudes related to the growth of an anaerobic bacteria in a digester.
Finally, it should be remembered that the design and construction of a laboratory-scale biodigester, due to its economic viability, is a tool available to engineering students for the development of their knowledge and learning.

Author Contributions

Conceptualization, S.B.-E. and A.R.-M.; methodology, S.B.-E. and A.R.-M.; software, S.B.-E. and A.R.-M.; validation, F.L., A.R.-M., S.B.-E., J.V.-R. and C.M.-P.; formal analysis, S.B.-E. and A.R.-M.; investigation, S.B.-E. and A.R.-M.; resources, A.R.-M. and J.V.-R.; data curation, F.L. and A.R.-M.; writing—original draft preparation, F.L. and A.R.-M.; writing—review and editing, F.L., C.M.-P. and J.V.-R.; visualization, S.B.-E., F.L. and A.R.-M.; supervision, F.L., C.M.-P., J.V.-R. and A.R.-M.; project administration, S.B.-E., F.L. and A.R.-M.; funding acquisition, J.V.-R., F.L. and A.R.-M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable for studies not involving humans or animals.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

In this section, please provide details regarding where data supporting reported results can be found, including links to publicly archived datasets analyzed or generated during the study.

Acknowledgments

This research has been co-funded by the INTERREG V-A Cooperation Spain-Portugal MAC (Madeira-Azores-Canarias) programme MITIMAC project MAC2/1.1a/263.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. SCILAB Source Code

clear; clc;
mu_max=.1/3600;
k_s=100; k_d=.00001; k_x=.00000150;
V_1=.01; V_2=3.99; V_T=V_1+V_2;
S_min=k_d*k_s/(mu_max-k_d); S_in=130;
dt=100; time=200*24*3600; Delta_t=.05*24*3600;
q=round(Delta_t/dt);
t=0:dt:time;
X=zeros(1,length(t)); S=zeros(1,length(t));
X(1)=.1; S(1)=0; S_in=S_in; X_in=.1;
e=1; f=0; g=0;
 
for k=2:length(t)
 
  f=f+1;
  if f>=q //Impulse input
     S(k-1)=(V_2*S(k-1)+V_1*S_in)/V_T;
     X(k-1)=(V_2*X(k-1)+V_1*X_in)/V_T;
     f=0;
  end
 
  //Prediction step.
  X(k)=(1+(mu_max*S(k-1)/(k_s+S(k-1))-k_d)*dt)*X(k-1);
  S(k)=(1-k_x*X(k-1)/(k_s+S(k-1))*dt)*S(k-1);
 
  //Initialization and recursive correction steps.
  While (e>=.01)
     x=X(k);
     s=S(k);
     X(k)=X(k-1)+dt/2*((mu_max*S(k-1)/(k_s+S(k-1))-k_d)*X(k-1)+ +(mu_max*S(k)/(k_s+S(k))-k_d)*X(k));
     S(k)=S(k-1)-k_x*dt/2*(S(k-1)/(k_s+S(k-1))*X(k-1)+S(k)/(k_s+S(k))*X(k));
     e=sqrt((x-X(k))^2+(s-S(k))^2);
  end
end

Appendix B. Microcontroller Source Code (Arduino)

// *****Digital inputs for pH y ORP****
#include <SoftwareSerial.h> //add the soft serial libray
#define rxph 4 //set the RX pin to pin 2
#define txph 5 //set the TX pin to pin 3
#define rxorp 2#define txorp 3
// ****pH data***
SoftwareSerial phserial(rxph, txph); //enable the soft serial port
String inputstringph = ““; // string to hold incoming data from the PC
String sensorstringph = ““; //a string to hold the data from the Atlas Scientific product
boolean input_stringcompleteph = false;//have we received all the data from the PC
boolean sensor_stringcompleteph = false; //have we received all the data from the Atlas Scientific product
// ****ORP data***
SoftwareSerial orpserial(rxorp, txorp);
String inputstringorp = ““; //a string to hold incoming data from the PC
String sensorstringorp = ““; //a string to hold the data from the Atlas Scientific product
boolean input_stringcompleteorp = false; //have we received all the data from the PC
boolean sensor_stringcompleteorp = false;//have we received all the data from the Atlas Scientific product
// ****Analog inputs***
String sensorpresion=““;//presion
float ntcbiodig=0; //termistor NTC lectura sistema de control
float ntcext=0; //termistor NTC lectura entorno
int analogPin1 = A1; //definimos los pines de entrada para la temperatura entorno
int analogPin2 = A2;//definimos los pines de entrada para la temperatura sistema
// **** SETUP ****
void setup(){ //set up the hardware
 
Serial.begin(9600);
 
phserial.begin(9600); //set baud rate for software serial port to 38400
 
orpserial.begin(9600);
 
phserial.print(“r\r”);
 
orpserial.print(“r\r”);
 
phserial.print(“r\r”);
 
orpserial.print(“r\r”);
 
inputstringph.reserve(5); //set aside some bytes for receiving data from the PC
 
sensorstringph.reserve(30); //set aside some bytes for receiving data from Atlas Scientific product
 
inputstringorp.reserve(5);//set aside some bytes for receiving data from the PC
 
sensorstringorp.reserve(30);
 
pinMode(analogPin1, INPUT); //def de los pines de entrada
 
pinMode(analogPin2, INPUT);
}
//**** loop ****
 
void loop(){
// ******Digital inputs *****
 
delay(60000);
 
input_stringcompleteph = true;
 
input_stringcompleteorp = true;
 
inputstringph=“r\r”;
 
inputstringorp=“r\r”;
 
bool mandar=true;
 
 
// ******pH*****
 
if (input_stringcompleteph){ //if a string from the PC has been recived in its entierty
 
phserial.print(inputstringph); //send that string to the Atlas Scientific product
 
inputstringph = ““; //clear the string:
 
input_stringcompleteph = false;//reset the flage used to tell if we have recived a completed string from the PC
 
}
 
while(mandar){
 
phserial.listen();
 
while (phserial.available()>0) { //while a char is holding in the serial buffer
 
char inchar = (char)phserial.read();//get the new char
 
sensorstringph += inchar; //add it to the sensorstringph
 
if (inchar == ‘\r’) {sensor_stringcompleteph = true;} //if the incoming character is a <CR>, set the flag
 
}
 
if (sensor_stringcompleteph){ //if a string from the Atlas Scientific product has been received in its entirety
 
Serial.print(sensorstringph); //use the hardware serial port to send that data to the PC
 
sensorstringph = ““; //clear the string:
 
sensor_stringcompleteph = false; //reset the flag used to tell if we have received a completed string from the Atlas Scientific product
 
mandar=false;
 
}
}
 
 
// ******ORP*****
 
if (input_stringcompleteorp){ //if a string from the PC has been recived in its entierty
 
orpserial.print(inputstringorp); //send that string to the Atlas Scientific product
 
inputstringorp = ““; //clear the string:
 
input_stringcompleteorp = false;//reset the flage used to tell if we have recived a completed string from the PC
 
}
 
mandar=true;
 
while(mandar){
orpserial.listen();
 
 
while (orpserial.available()>0) { //while a char is holding in the serial buffer
 
char inchar = (char)orpserial.read();//get the new char
 
sensorstringorp += inchar; //add it to the sensorstringorp
 
if (inchar == ‘\r’) {sensor_stringcompleteorp = true;} //if the incoming character is a <CR>, set the flag
 
}
 
 
// ******Analogic inputs*****
// ******THERMISTOR*****
 
int ntcext=analogRead(analogPin2); // leemos
// ******TERMISTOR SISTEMA*****
 
int ntcbiodig=analogRead(analogPin1); // leemos
//****print ****
 
int sensorpresion= analogRead(A0); // lee valor de presion
 
if (sensor_stringcompleteorp){//if a string from the Atlas Scientific product has been received in its entirety
 
Serial.print(sensorstringph);//use the hardware serial port to send that data to the PC
 
Serial.print(“,”); //use the hardware serial port to send that data to the PC
 
Serial.print(sensorstringorp);
 
Serial.print(“,”);
 
Serial.print(sensorpresion); // imprime valor de presion
 
Serial.print(“,”);
 
Serial.print(ntcext);
 
Serial.print(“,”);
 
Serial.println(ntcbiodig);
 
sensorstringorp = ““; //clear the string:
 
sensor_stringcompleteorp = false; //reset the flag used to tell if we have received a completed string from the Atlas Scientific product
 
mandar=false;
 
}
 
}
}
 

Appendix C. Processing Source Code

import processing.serial.*;
import controlP5.*;
 
String nom_archivo=“name.txt”;
String hora_inicio=“Begining: “+hour()+”:”+minute()+”:”+second();
PrintWriter archivo;
ControlP5 cp5;
Serial myPort; // The serial port
boolean serialInited;
int xPos = 1; // horizontal position of the graph
int esquina_x=600;
int esquina_y=100;
int ancho=500;
int alto=300;
int pocision_alto=0;//posicion donde se colocan los datos en la grafica, menor que alto.
int guardar_datos=-1;
String dato_1=“hola”,dato_2=“adios”,dato_3=“hasta”,dato_4=“luego”,dato_5=“vengo”,dato_6=“voy”;
PShape bot;
int ancho_ventana=720,alto_ventana=400;
int x_r1=470,y_r1=100,dx_r1=200,dy_r1=230;
int mx_r1=10,my_r1=30,dy_t=30;
int x_r2=40,y_r2=100,dx_r2=385,dy_r2=130;
int mx_r2=10,my_r2=30;
int ancho_campo_texto=210;
int x_resto=40,y_resto=250,dy_resto=25;
int x_r3=40,y_r3=300,dx_r3=210,dy_r3=35;
int mx_r3=10,my_r3=30;
int ancho_boton=200,alto_boton=20;
int lf = 10; // Salto de linea en ASCII
int BAUD_RATE=9600;
String inString=null;
 
 
 
void setup () {
 
cp5 = new ControlP5(this);
 
 
cp5.addButton(“Save_data”)
 
.setValue(0)
 
.setPosition(x_resto,y_resto)
 
.setSize(ancho_boton,alto_boton)
 
;
 
 
 
cp5.addButton(“Stop_save_data”)
 
.setValue(0)
 
.setPosition(x_resto,y_resto+dy_resto)
 
.setSize(ancho_boton,alto_boton)
 
;
 
 
 
cp5.addTextfield(“nom_archivo”)
 
.setPosition(x_r2+mx_r2+ancho_campo_texto, y_r2+my_r2/3)
 
.setSize(150,25)
 
.setCaptionLabel(““)
 
.setColorBackground(color(50,50,50))
 
.setColorActive(color(100,0,0))
 
.setColorForeground(color(100,0,0))
 
.setFont(createFont(“arial”,15))
 
;
 
//size(ancho_ventana,alto_ventana);
size(720,400);
stroke(127,0,0);
//rect(esquina_x, esquina_y, ancho, alto);
 
println(Serial.list());
 
myPort = new Serial(this, Serial.list()[0], 9600);//Se crea una comunicación serial en el puerto 0, con 9600bd.
myPort.clear();
myPort.buffer(50);
inString = myPort.readStringUntil(lf);
inString = null;
 
 
background(0,255,200);
 
bot = loadShape(“ULPGC.svg”);
shape(bot, 20, 20, 70, 50);
 
String s_1=“University of Las Palmas de Gran Canaria”;
String s_2=“Laboratory-scale biodigester.”;
String s_3=“Saved data”;
String s_4=“T= “;
 
textSize(20);
fill(0, 102, 153);
text(s_1, 95, 45);
text(s_2, 95, 65);
text(s_3, 95, 205);
fill(255, 255, 255);
}
 
void draw () {
 
//background(0,255,200);
fill(155, 255, 255);
rect(x_r1,y_r1,dx_r1,dy_r1);
color(250,0,0,205);
fill(0, 102, 153);
text(“Measured variables”, x_r1+mx_r1, y_r1+my_r1);
text(“Tamb = “+dato_1+ “ºC”, x_r1+mx_r1, y_r1+my_r1+dy_t);
text(“Treact = “+dato_2+ “ºC”, x_r1+mx_r1, y_r1+my_r1+2*dy_t);
text(“P = “+dato_3+ “ kPa”, x_r1+mx_r1, y_r1+my_r1+3*dy_t);
text(“ORP = “+dato_4+ “ mV”, x_r1+mx_r1, y_r1+my_r1+4*dy_t);
text(“pH = “+dato_5, x_r1+mx_r1, y_r1+my_r1+5*dy_t);
text(“Aux = “+dato_6, x_r1+mx_r1, y_r1+my_r1+6*dy_t);
fill(155, 255, 255);
rect(x_r2,y_r2,dx_r2,dy_r2);
//rect(40, 130, 385, 130);
color(250,0,0,205);
fill(0, 102, 153);
text(“File name:”,x_r2+mx_r2, y_r2+my_r2);
text(nom_archivo,x_r2+mx_r2, y_r2+my_r2+dy_t);
text(hora_inicio,x_r2+mx_r2, y_r2+my_r2+2*dy_t);
text(“Saving: “+hour()+”:”+minute()+”:”+second(),x_r2+mx_r2, y_r2+my_r2+3*dy_t);
fill(255, 255, 255);
strokeWeight(1);
color(250,0,0,255);
strokeWeight(2);
fill(155, 255, 255);
rect(x_r3,y_r3,dx_r3,dy_r3);
color(250,0,0,205);
fill(0, 102, 153);
 
if(guardar_datos>0){
  color(255,255,0,0);
  fill(255, 0, 0);
  ellipse(x_r3+2*mx_r3, y_r3+my_r3/2, 15, 15);
  text(“Saving”,x_r3+4*mx_r3,y_r3+my_r3/1.25);
}
else{
  color(250,0,0,0);
  fill(0, 0, 0);
  ellipse(x_r3+2*mx_r3, y_r3+my_r3/2, 15, 15);
  text(“Showing”,x_r3+4*mx_r3,y_r3+my_r3/1.25);
}
}
 
void serialEvent (Serial myPort) {
// Toma una cadena de caracteres ASCII:
inString = myPort.readStringUntil(‘\n’);
inString = trim(inString);
 
if (inString != null) {
  println(inString);
  String[] lista = split(inString,”,”);
  dato_1=lista[0];
  dato_2=lista[1];
  dato_3=lista[2];
  dato_4=lista[3];
  dato_5=lista[4];
  dato_6=“---”;//lista[5];
 
 
if(guardar_datos>0){
  archivo.print(month( )+” “);
  archivo.print(day( )+” “);
  archivo.print(year( )+” “);
  archivo.print(hour( )+” “);
  archivo.print(minute( )+” “);
  archivo.print(second( )+” “);
  archivo.print(dato_1+” “);
  archivo.print(dato_2+” “);
  archivo.print(dato_3+” “);
  archivo.print(dato_4+” “);
  archivo.print(dato_5+” “);
  archivo.println(dato_6+” “);
  archivo.flush(); // Writes the remaining data to the file
}
}
}
 
public void Save_data(int theValue) {
 
println(“a button event from colorA: “+theValue);
 
//++dato_prueba;
 
archivo = createWriter(nom_archivo);
 
hora_inicio=“Begining: “+hour()+”:”+minute()+”:”+second();
 
if(guardar_datos>0){int r=0;}
 
else{
 
guardar_datos=−1*guardar_datos;
 
}
}
public void Stop_save_data(int theValue) {
 
println(“a button event from colorA: “+theValue);
 
//++dato_prueba;
 
guardar_datos=-1;
 
archivo.flush(); // Writes the remaining data to the file
 
archivo.close(); // Finishes the file
 
nom_archivo=“name.txt”;
}
 

References

  1. Brito-Espino, S.; Ramos-Martín, A.; Pérez-Báez, S.O.; Mendieta-Pino, C.; Leon-Zerpa, F. A Framework Based on Finite Element Method (FEM) for Modelling and Assessing the Affection of the Local Thermal Weather Factors on the Performance of Anaerobic Lagoons for the Natural Treatment of Swine Wastewater. Water 2021, 13, 882. [Google Scholar] [CrossRef]
  2. Cano, E.; Ruiz, J.G.; Garcia, I.A. Integrating a learning constructionist environment and the instructional design approach into the definition of a basic course for embedded systems design. Comput. Appl. Eng. Educ. 2015, 23, 36–53. [Google Scholar] [CrossRef]
  3. Garcia, I.A.; Cano, E.M. Designing and implementing a constructionist approach for improving the teaching-learning process in the embedded systems and wireless communications areas. Comput. Appl. Eng. Educ. 2014, 22, 481–493. [Google Scholar] [CrossRef]
  4. Available online: http://www.atago.net/Spanish/download.htmlRX-7000i (accessed on 11 June 2021).
  5. Refractómetro Digital Automático Atago. Available online: http://www.atago.net/Spanish/download.html (accessed on 14 June 2021).
  6. Jagnow, G.; Dawind, W. Biotecnología: Introducción Con Experimentos Modelo; Acribia S.A.: Zaragoza, Spain, 1991. [Google Scholar]
  7. Liu, C.-G.; Xue, C.; Lin, Y.-H.; Bai, F.-W. Redox potential control and applications in microaerobic and anaerobic fermentations. Biotechnol. Adv. 2013, 31, 257–265. [Google Scholar] [CrossRef]
  8. Madsen, M.; Holm-Nielsen, J.B.; Esbensen, K.H. Monitoring of anaerobic digestion processes: A review perspective. Renew. Sustain. Energy Rev. 2011, 15, 3141–3155. [Google Scholar] [CrossRef] [Green Version]
  9. Mekic, E.; Djokic, I.; Zejnelagic, S.; Matovic, A. Constructive approach in teaching of voip in line with good laboratory and manufacturing practice. Comput. Appl. Eng. Educ. 2016, 24, 277–287. [Google Scholar] [CrossRef]
  10. Pantaleo, A.; De Gennaro, B.; Shah, N. Assessment of optimal size of anaerobic co- digestion plants: An application to cattle farms in the province of bari (Italy). Renew. Sustain. Energy Rev. 2013, 20, 57–70. [Google Scholar] [CrossRef]
  11. Atlas Scientific. Atlas Scientific. Orpatlas. Available online: https://atlas-scientific.com/?gclid=EAIaIQobChMIvbPk8o3Z8gIVied3Ch03AgqaEAAYASAAEgKLzPD_BwE (accessed on 19 July 2021).
  12. Sorathia, K.; Servidio, R. Learning and experience: Teaching tangible interaction & edutainment. Procedia—Soc. Behav. Sci. 2012, 64, 265–274. [Google Scholar]
  13. Taylhardat Arjona, L.A. El biogas. Fundamentos e Infraestructura Rural; Instituto de Ingenieria Agri- cola; Facultad de Agronomia U.C.V: Maracay, Venezuela, 1986. [Google Scholar]
  14. Leon, F.; Ramos, A.; Vaswani, J.; Mendieta, C.; Brito, S. Climate Change Mitigation Strategy through Membranes Replacement and Determination Methodology of Carbon Footprint in Reverse Osmosis RO Desalination Plants for Islands and Isolated Territories. Water 2021, 13, 293. [Google Scholar] [CrossRef]
  15. Products for USB Sensing and Control. Products for Usb Sensing and Control. Available online: www.phidgets.com (accessed on 1 July 2021).
  16. Parralejo, A.; Royano, L.; González, J.; González, J. Small scale biogas production with animal excrement and agricultural residues. Ind. Crops Prod. 2019, 131, 307–314. [Google Scholar] [CrossRef]
  17. Holm-Nielsen, J.; Seadi, T.A.; Oleskowicz-Popiel, P. The future of anaerobic digestion and biogas utilization. Bioresour. Technol. 2009, 100, 5478–5484. [Google Scholar] [CrossRef]
  18. Park, J.H.; Park, J.H.; Lee, S.H.; Jung, S.P.; Kim, S.H. Enhancing anaerobic digestion for rural wastewater treatment with granular activated carbon (GAC) supplementation. Bioresour. Technol. 2020, 315, 123890. [Google Scholar] [CrossRef] [PubMed]
  19. Jiang, Y.; Bebee, B.; Mendoza, A.; Robinson, A.K.; Zhang, X.; Rosso, D. Energy footprint and carbon emission reduction using off-the-grid solar-powered mixing for lagoon treatment. J. Environ. Manag. 2018, 205, 125–133. [Google Scholar] [CrossRef] [PubMed]
  20. Duan, N.; Zhang, D.; Khoshnevisan, B.; Kougias, P.G.; Treu, L.; Liu, Z.; Lin, C.; Liu, H.; Zhang, Y.; Angelidaki, I. Human waste anaerobic digestion as a promising low-carbon strategy: Operating performance, microbial dynamics and environmental footprint. J. Clean. Prod. 2020, 256, 120414. [Google Scholar] [CrossRef]
  21. Mendieta-Pino, C.A.; Ramos-Martin, A.; Perez-Baez, S.O.; Brito-Espino, S. Management of slurry in Gran Canaria Island with full-scale natural treatment systems for wastewater (NTSW). One year experience in livestock farms. J. Environ. Manag. 2019, 232, 666–678. [Google Scholar]
  22. Muga, H.; Mihelcic, J. Sustainability of wastewater treatment technologies. J. Environ. Manag. 2008, 88, 437–447. [Google Scholar] [CrossRef]
  23. Wu, B.; Chen, Z. An integrated physical and biological model for anaerobic lagoons. Bioresour. Technol. 2011, 102, 5032–5038. [Google Scholar] [CrossRef] [PubMed]
  24. Wu, B. Advances in the use of CFD to characterize, design and optimize bioenergy systems. Comput. Electron. Agric. 2013, 93, 195–208. [Google Scholar] [CrossRef]
  25. Donoso-Bravo, A.; Sadino-Riquelme, C.; Gómez, D.; Segura, C.; Valdebenito, E.; Hansen, F. Modelling of an anaerobic plug-flow reactor. Process analysis and evaluation approaches with non-ideal mixing considerations. Bioresour. Technol. 2018, 260, 95–104. [Google Scholar]
  26. Rajeshwari, K.; Balakrishnan, M.; Kansal, A.; Lata, K.; Kishore, V. State-of-the-art of anaerobic digestion technology for industrial wastewater treatment. Renew. Sustain. Energy Rev. 2000, 4, 135–156. [Google Scholar] [CrossRef]
  27. Lauwers, J.; Appels, L.; Thompson, I.P.; Degrève, J.; Impe, J.F.V.; Dewil, R. Mathematical modelling of anaerobic digestion of biomass and waste: Power and limitations. Prog. Energy Combust. 2013, 39, 383–402. [Google Scholar] [CrossRef] [Green Version]
  28. Wade, M.; Harmand, J.; Benyahia, B.; Bouchez, T.; Chaillou, S.; Cloez, B. Perspectives in mathematical modelling for microbial ecology. Ecol. Model. 2016, 321, 64–74. [Google Scholar] [CrossRef]
  29. Batstone, D.; Keller, J.; Angelidaki, I.; Kalyuzhnyi, S.; Pavlostathis, S.; Rozzi, A.; Sanders, W.; Siegrist, H.; Vavilin, V. The IWAAnaerobic Digestion Model No 1 (ADM1). Water Sci. Technol. 2002, 45, 65–73. Available online: https://library.lanl.gov/cgi-bin/getfile?00285556.pdf (accessed on 1 September 2020). [CrossRef]
  30. Kleerebezem, R.; van Loosdrecht, M.C.M. Critical analysis of some concepts proposed in ADM1. Water Sci. Technol. 2006, 54, 51–57. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Li, D.; Lee, I.; Kim, H. Application of the linearized ADM1 (LADM) to lab-scale anaerobic digestion system. J. Environ. Chem. Eng. 2021, 9, 105193. [Google Scholar] [CrossRef]
  32. Fleming, J.G. Novel Simulation of Anaerobic Digestion Using Computational Fluid Dynamics. Ph.D. Thesis, North Carolina State University, Raleigh, NC, USA, 2002. [Google Scholar]
  33. Goodarzi, D.; Sookhak Lari, K.; Mossaiby, F. Thermal effects on the hydraulic performance of sedimentation ponds. J. Water Process. Eng. 2020, 33, 101100. [Google Scholar] [CrossRef]
  34. Brito-Espino, S.; Ramos-Martín, A.; Pérez-Báez, S.; Mendieta-Pino, C. Application of a mathematical model to predict simultaneous reactions in anaerobic plug-flow reactors as a primary treatment for constructed wetlands. Sci. Total Environ. 2020, 713, 136244. [Google Scholar] [CrossRef]
  35. Mahmudul, H.; Rasul, M.; Akbar, D.; Narayanan, R.; Mofijur, M. A comprehensive review of the recent development and challenges of a solar-assisted biodigester system. Sci. Total Environ. 2021, 753, 141920. [Google Scholar] [CrossRef] [PubMed]
  36. Atelge, M.; Atabani, A.; Banu, J.R.; Krisa, D.; Kaya, M.; Eskicioglu, C.; Kumar, G.; Lee, C.; Yildiz, Y.; Unalan, S.; et al. A critical review of pretreatment technologies to enhance anaerobic digestion and energy recovery. Fuel 2020, 270, 117494. [Google Scholar] [CrossRef]
  37. Tumilar, A.S.; Milani, D.; Cohn, Z.; Florin, N.; Abbas, A. A Modelling Framework for the Conceptual Design of Low-Emission Eco-Industrial Parks in the Circular Economy: A Case for Algae-Centered Business Consortia. Water 2021, 13, 69. [Google Scholar] [CrossRef]
  38. Haßler, S.; Ranno, A.M.; Behr, M. Finite-element formulation for advection–reaction equations with change of variable and discontinuity capturing. Comput. Methods Appl. Mech. Eng. 2020, 369, 113171. [Google Scholar] [CrossRef]
  39. Mirza, I.A.; Akram, M.S.; Shah, N.A.; Imtiaz, W.; Chung, J.D. Analytical solutions to the advection-diffusion equation with Atangana-Baleanu time-fractional derivative and a concentrated loading. Alex. Eng. J. 2021, 60, 1199–1208. [Google Scholar] [CrossRef]
  40. Singh, S.; Bansal, D.; Kaur, G.; Sircar, S. Implicit-explicit-compact methods for advection diffusion reaction equations. Comput. Fluids 2020, 212, 104709. [Google Scholar] [CrossRef]
  41. Zeng, L.; Chen, G. Ecological degradation and hydraulic dispersion of contaminant in wetland. Ecol. Model. 2011, 222, 293–300. [Google Scholar] [CrossRef]
  42. Bozkurt, S.; Moreno, L.; Neretnieks, I. Long-term processes in waste deposits. Sci. Total Environ. 2000, 250, 101–121. [Google Scholar] [CrossRef]
  43. Song, L.; Li, P.W.; Gu, Y.; Fan, C.M. Generalized finite difference method for solving stationary 2D and 3D Stokes equations with a mixed boundary condition. Comput. Math. Appl. 2020, 80, 1726–1743. [Google Scholar] [CrossRef]
  44. Ukai, S. A solution formula for the Stokes equation in Rn+. Commun. Pure Appl. Math. 1987, 40, 611–621. [Google Scholar] [CrossRef] [Green Version]
  45. Reddy, J.; Gartling, D. The Finite Element Method in Heat Transfer and Fluid Dynamics, 3rd ed.; CRC Press: Boca Raton, FL, USA, 2010; pp. 1–483. [Google Scholar]
  46. Alvarez-Hostos, J.C.; Bencomo, A.D.; Puchi-Cabrera, E.S.; Fachinotti, V.D.; Tourn, B.; Salazar-Bove, J.C. Implementation of a standard stream-upwind stabilization scheme in the element-free Galerkin based solution of advection-dominated heat transfer problems during solidification in direct chill casting processes. Eng. Anal. Bound. Elem. 2019, 106, 170–181. [Google Scholar] [CrossRef]
  47. Guldentops, G.; Van Dessel, S. A numerical and experimental study of a cellular passive solar façade system for building thermal control. Sol. Energy 2017, 149, 102–113. [Google Scholar] [CrossRef]
  48. Lawrence, M.G. The Relationship between Relative Humidity and the Dewpoint Temperature in Moist Air: A Simple Conversion and Applications. Bull. Am. Meteorol. Soc. 2005, 86, 225–234. [Google Scholar] [CrossRef]
  49. Çengel, Y. Heat Transfer: A Practical Approach. In McGraw-Hill Series in Mechanical Engineering; McGraw Hill Books: London, UK, 2003. [Google Scholar]
  50. Walton, G.N. Thermal Analysis Research Program Reference Manual; NBSIR, Department of Energy, Office of Building Energy Research and Development: Washington, DC, USA, 1983. [Google Scholar]
  51. Monod, J. The Growth of Bacterial Cultures. Annu. Rev. Microbiol. 1949, 3, 371–394. [Google Scholar] [CrossRef] [Green Version]
  52. Rosso, L.; Lobry, J.; Flandrois, J. An Unexpected Correlation between Cardinal Temperatures of Microbial Growth Highlighted by a New Model. J. Theor. Biol. 1993, 162, 447–463. [Google Scholar] [CrossRef] [PubMed]
  53. Herus, V.A.; Ivanchuk, N.V.; Martyniuk, P.M. A System Approach to Mathematical and Computer Modeling of Geomigration Processes Using Freefem++ and Parallelization of Computations. Cybern Syst. Anal. 2018, 54, 284–292. [Google Scholar] [CrossRef]
  54. Donoso-Bravo, A.; Bandara, W.; Satoh, H.; Ruiz-Filippi, G. Explicit temperature-based model for anaerobic digestion: Application in domestic wastewater treatment in a UASB reactor. Bioresour. Technol. 2013, 133, 437–442. [Google Scholar] [CrossRef] [Green Version]
  55. Donoso-Bravo, A.; Retamal, C.; Carballa, M.; Ruiz-Filippi, G.; Chamy, R. Influence of temperature on the hydrolysis, acidogenesis and methanogenesis in mesophilic anaerobic digestion: Parameter identification and modeling application. Water Sci. Technol. 2009, 60, 9–17. [Google Scholar] [CrossRef]
  56. Wang, R.; Lv, N.; Li, C.; Cai, G.; Pan, X.; Li, Y.; Zhu, G. Novel strategy for enhancing acetic and formic acids generation in acidogenesis of anaerobic digestion via targeted adjusting environmental niches. Water Res. 2021, 193, 116896. Available online: https://www.sciencedirect.com/science/article/pii/S0043135421000944 (accessed on 23 June 2021). [CrossRef]
  57. Weißbach, M.; Drewes, J.E.; Koch, K. Application of the oxidation reduction potential (ORP) for process control and monitoring nitrite in a Coupled Aerobic-anoxic Nitrous Decomposition Operation (CANDO). Chem. Eng. J. 2018, 343, 484–491. Available online: https://www.sciencedirect.com/science/article/pii/S1385894718303929 (accessed on 30 June 2021). [CrossRef]
  58. Ao, T.; Chen, L.; Zhou, P.; Liu, X.; Li, D. The role of oxidation-reduction potential as an early warning indicator, and a microbial instability mechanism in a pilot-scale anaerobic mesophilic digestion of chicken manure. Renew. Energy 2021, 179, 223–232. Available online: https://www.sciencedirect.com/science/article/pii/S0960148121010521 (accessed on 1 July 2021). [CrossRef]
Figure 1. Basic diagram of experimental design.
Figure 1. Basic diagram of experimental design.
Water 13 02409 g001
Figure 2. Thermistor NTC Vishay, dimensions in millimeters.
Figure 2. Thermistor NTC Vishay, dimensions in millimeters.
Water 13 02409 g002
Figure 3. Arduino Uno.
Figure 3. Arduino Uno.
Water 13 02409 g003
Figure 4. Photos of the reactor and sensors (a) and the circuits with the controls systems (b).
Figure 4. Photos of the reactor and sensors (a) and the circuits with the controls systems (b).
Water 13 02409 g004
Figure 5. Circuits and control system design.
Figure 5. Circuits and control system design.
Water 13 02409 g005
Figure 6. Processing PC interface data logger.
Figure 6. Processing PC interface data logger.
Water 13 02409 g006
Figure 7. Concentration of biomass and substrate where μ m a x = 0.75   ( h 1 ) ,   K s = 10 ,   β s = 1.8 , X(t=0) is 0.1 (g/m3) and S(t = 0) is 4 (g/m3).
Figure 7. Concentration of biomass and substrate where μ m a x = 0.75   ( h 1 ) ,   K s = 10 ,   β s = 1.8 , X(t=0) is 0.1 (g/m3) and S(t = 0) is 4 (g/m3).
Water 13 02409 g007
Figure 8. ORP, pH, and temperature in the first period.
Figure 8. ORP, pH, and temperature in the first period.
Water 13 02409 g008
Figure 9. ORP, pH, and Temperature in the second period.
Figure 9. ORP, pH, and Temperature in the second period.
Water 13 02409 g009
Figure 10. ORP, pH, and Temperature in the third period.
Figure 10. ORP, pH, and Temperature in the third period.
Water 13 02409 g010
Table 1. Characteristics of the thermistor NTC, Vishay.
Table 1. Characteristics of the thermistor NTC, Vishay.
ParameterValueUnits
Resistence value at 25 °C10 K
Tolerance of R25±3%
B25/85 (Beta)3984K
Temperatura range of operation−25 to 105°C
Table 2. Laboratory-scale biodigester; performance periods, operating volume, and Brix measures.
Table 2. Laboratory-scale biodigester; performance periods, operating volume, and Brix measures.
StageDateFeeding/Evacuation Brix Grades Remarks
(mL)1st Lecture2nd LectureAverage
113 June30020.30 20.30
218 June7520.3120.2920.30
323 June5020.3320.2420.26Addition NaOH (↑alkalinity)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Brito-Espino, S.; Leon, F.; Vaswani-Reboso, J.; Ramos-Martin, A.; Mendieta-Pino, C. Proposal of a Laboratory-Scale Anaerobic Biodigester for Introducing the Monitoring and Sensing Techniques, as a Potential Learning Tool in the Fields of Carbon Foot-Print Reduction and Climate Change Mitigation. Water 2021, 13, 2409. https://doi.org/10.3390/w13172409

AMA Style

Brito-Espino S, Leon F, Vaswani-Reboso J, Ramos-Martin A, Mendieta-Pino C. Proposal of a Laboratory-Scale Anaerobic Biodigester for Introducing the Monitoring and Sensing Techniques, as a Potential Learning Tool in the Fields of Carbon Foot-Print Reduction and Climate Change Mitigation. Water. 2021; 13(17):2409. https://doi.org/10.3390/w13172409

Chicago/Turabian Style

Brito-Espino, Saulo, Federico Leon, Jenifer Vaswani-Reboso, Alejandro Ramos-Martin, and Carlos Mendieta-Pino. 2021. "Proposal of a Laboratory-Scale Anaerobic Biodigester for Introducing the Monitoring and Sensing Techniques, as a Potential Learning Tool in the Fields of Carbon Foot-Print Reduction and Climate Change Mitigation" Water 13, no. 17: 2409. https://doi.org/10.3390/w13172409

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