MATLAB Algorithm Research Articles

Optimal configuration and operation for user-side energy storage considering lithium-ion battery degradation

Battery energy storage systems (BESSs) have been widely employed on the user-side such as buildings, residential communities, and industrial sites due to their scalability, quick response, and design flexibility. However, cell degradation is caused by the charging and discharging of batteries, which reduces the economy of BESSs. For the optimal configuration and operation of BESSs, the battery degradation process should be integrated into the optimal problem by considering the equivalent full cycles over the optimization time horizon. This paper proposes a two-layer optimization frame to estimate and improve the net profit of BESSs in the whole life cycle, the outer layer optimizes the rated capacity and power of BESSs, and the inner layer optimizes the daily operation strategy of BESSs. To this end, the semi-empirical degradation model of lithium-ion batteries and economic models of BESSs are embedded into the optimization frame. Particle swarm optimization (PSO) algorithm and fmincon toolbox of MATLAB are adopted to solve the two-layer frame to maximize the net profit of BESSs. Simulation results of the BESS for a typical industrial user in China demonstrate that the proposed frame can effectively improve the net profit of BESSs.

Use of Automated Kinematic Diadochokinesis Analysis to Identify Potential Indicators of Speech Motor Involvement in Children With Cerebral Palsy.

This study examined multiple variables obtained from an automated measure of lip movement during a diadochokinesis (DDK) task to identify those with potential to detect mild speech motor involvement in school-age children diagnosed with cerebral palsy (CP). Eight children with CP and high speech intelligibility and a matched group of eight children with typical development (TD) completed a DDK task while their lip and jaw movements were recorded. A custom MATLAB algorithm was used to automatically extract 23 kinematic measures of children's lip movements during production of the DDK sequences. Mann-Whitney U tests were used to compare groups on the kinematic measures, and receiver operating characteristic (ROC) analysis was used to evaluate the diagnostic accuracy of measures that significantly differed between groups. Five of the 23 kinematic variables differed significantly between the CP and TD groups. These were two measures of overall DDK performance (i.e., duration of the DDK sequence and number of cycles) and three spatial and temporal measures of lip movement. Duration of the DDK sequence and the mean displacement of the lips across cycles had the highest diagnostic accuracy, differentiating CP and TD groups with 88% sensitivity and 88% specificity. Automatically derived kinematic measures of DDK sequences differentiated children with CP and high intelligibility from typically developing children. Future research is needed to determine the clinical utility of these measures for detecting speech motor impairment.

Parameter Matching and Structural Optimization Design of H-SFCL in MMC Ship MVDC System

As a high-density and high-load system, ship DC system has the characteristics of short power supply line, low system damping. When a short-circuit fault occurs, it will cause serious damage and the direct current circuit breaker (DCCB) is required to cut off the fault current quickly. In our previous paper, a hybrid-type superconducting fault current limiter (H-SFCL) has been proposed, which can help DCCB reliably remove faults and realize fast reclosing. On this basis, the H-SFCL is selected and applied in ship medium voltage direct current (MVDC) system. In this paper, not only the current limiting principle of H-SFCL is introduced, but also its current limiting characteristics in modular multilevel converter (MMC) ship MVDC system are simulated and analyzed. In addition, the paraments selection and coil structural design of H-SFCL are optimized to improve the economy and ensure the safety and stability of the coil. Firstly, the parameter matching method is proposed based on the optimization algorithm in Matlab, and then the electromagnetic force coupling model in COMSOL is built to complete the structural optimization of the superconducting coil. The effectiveness of the optimization results is verified by stress analysis.

Hyper-Spherical Search Algorithm for Maximum Power Point Tracking of Solar Photovoltaic Systems under Partial Shading Conditions

Maximum Power Point Tracking (MPPT) for photovoltaic systems has widely been studied. However, studying the impact of partial shading (due to buildings, trees, clouds, adjacent arrays, and so on) on the MPPT and achieving the global maximum is still a challenging topic. This is because not only there is a global maximum power point (GMPP) but also there are several local maximums due to the non-linear nature of the power-voltage curve after partial shading. Therefore, the conventional MPPT methods fail to track the GMPP and are commonly trapped at one of the local maximum power points. This study utilizes the Hyper-Spherical Search (HSS) algorithm in MATLAB to achieve the GMPP while improving the efficiency, convergence speed, dynamic response, and reducing the losses. To track the GMPP using the HSS algorithm, the output power for the array (PPV) is defined as the objective function, and the duty cycle of the DC-DC converter is selected as the particles’ position (control variable). The performance of the proposed algorithm has been studied in three different partial shading patterns, and the simulation results confirm the capability of the algorithm in the rapid tracking of GMPP points. In addition, if the GMPP position changes over time, it will track the new GMPP with minimal oscillations. The proposed method along with PSO, P&O, ABC, and Dragon algorithms has been applied for various scenarios, and the obtained results using HSS have been compared with the four mentioned algorithms, which confirmed the effectiveness of the proposed method. Briefly, the advantages of the HSS algorithm in finding GMPP can be stated as simple implementation with few parameters, strong exploration, and exploitation during the tracking process, fast-tracking and low fluctuations during tracking, low oscillation at steady-state, high dynamic efficiency, and non-convergence to LMPP points.

Smoke Emission and Distribution Characteristics of Overloaded Wire Insulations under Microgravity

With the development of manned space technology, spacecraft safety and fire prevention have attracted more and more attention. Due to the disappearance of buoyancy in microgravity, the fire early monitoring, detection and alarm technologies designed based on the ground experimental results are not suitable for spacecraft. It is necessary to develop fire early monitoring technology in microgravity. Smoke is an important early monitoring signal for fire prevention both in normal gravity and microgravity. Under microgravity, fire is mostly caused by overload or aging of wire insulations. In order to study the smoke emission characteristics of wire insulations under microgravity, we carried out the overload experiments of wire insulations on board the SJ-10 Chinese recoverable satellite. The smoke generation characteristics captured by laser extinction methods, and a large number of experimental data in the real microgravity environment were obtained for the first time. In this paper, the smoke volume fraction in the early and axisymmetric stages of smoke emission from the wire insulating layer are obtained by using the method of Abel transform and convolution, and the MATLAB algorithm program is compiled. In the later stage of smoke emission, it does not show axisymmetric distribution, but the laser extinction results can be used for obtaining the smoke emission trajectory. According to the results, edge jet smoke emission modes in the early stage of ignition of wire insulation in microgravity are quantitatively analyzed. The effects of insulation thickness, overload current and insulation material on smoke emission are discussed.

Study on ignition process and flame expansion and propagation characteristics in jet-cooled pilot flameholders using image processing techniques

Wall-type pilot flameholder has been used in aero-engines to effectively create a stable ignition source and achieve an excellent ignition performance in high-velocity flows. The ignition delay, flame propagation, and flame boundary are crucial characteristics to evaluate the performance of pilot flameholders. Previous work displayed the significantly changed ignition performance and flow features of the wall flameholder when jet cooling schemes were adopted, indicating that the ignition process and the flame features strongly depend on the cooling jet and the optimal design of jet cooling schemes needs further investigation. Therefore, in this study, two practical image processing algorithms are developed to characterize the flame development of the ignition process, reveal the effect of cooling schemes on the ignition delay time, and quantitatively analyze the flame area and expansion ratio of the wall flameholder during stable combustion. Two image processing methods proposed based on MATLAB algorithms exhibit effectiveness in studying the ignition process and flame expansion characteristics. The ignition process is divided into four phases: kernel generation, kernel propagation, flame growth, and stable combustion. Meanwhile, the ignition delay time is identified by counting the total pixels of the flame projected region during the ignition process. The edge detection and pixel statistics of the flame projected region in the binary luminosity image show the flame boundary and expansion ratio under various working conditions. It was found that the cooling jet angle has a significant impact on the ignition process, ignition delay time, flame propagation, and steady flame features. The pressure-driven gas jet cooling scheme with a smaller jet angle and the external-inhaled air cooling scheme with a larger jet angle present a better capability of ignition and flame expansion and propagation.

Indifference Computer Experiment for Mathematical Identification of Two Variables

In order to understand the two types of nonlinear differential equation problems in engineering dynamics, the author proposes a numerical analysis method for the two types of nonlinear differential equations based on computer simulation. This method establishes the MATLAB algorithm structure of the numerical solution of the fourth-order fixed-step Runge-Kutta and Lorenz models, discusses the error control in the case of variable step size, and plots the numerical solutions of the Lorenz system based on MATLAB in two-dimensional and three-dimensional space graphics. The x -direction displacement and y -direction displacement data are extracted from the Lorenz equation as iterative samples of the model, the regression curve obtained after iteration has a slope of 0.996, and the iterative regression model reflects the basic characteristics of the data well. This method presents the basic idea of numerical solution verification within acceptable error limits. For solving engineering problems with differential equations as mathematical models, an effective numerical solution method is provided, and further discussion on the numerical solutions of partial differential equations is of great significance.

The Effects of Compressibility on the Performance and Modal Structures of a Sweeping Jet Emitted from Various Scales of a Fluidic Oscillator

Investigations of fluidic oscillators, or sweeping jet actuators, have primarily been conducted within the incompressible flow regime, which limits the accuracy of estimating fluidic oscillator performance for compressible flows. The objective of this study was to evaluate the effects of gas compressibility on the performance of a fluidic oscillator. A commonly used fluidic oscillator geometry (the Bray geometry) was scaled to five different sizes, 3D printed, and tested over a range of air flow rates. High-speed Schlieren images captured the sweeping jet exiting the fluidic oscillators, and custom MATLAB algorithms were used to calculate the oscillation frequencies and angles. A spectral proper orthogonal decomposition (SPOD) method was used to identify and compare the mode structures within the flow fields. All the results were compared using dimensionless parameters to observe performance trends. The results showed that the oscillation frequencies were directly proportional to the flow rate, while the oscillation angles were inversely proportional to the flow rate, regardless of scale size. The angular velocities were not proportional to the flow rate or scale size and exhibited maxima within the evaluated ranges. For all scale sizes, the mode structures were symmetric across the centerlines of the fluidic oscillators and extended further beyond the fluidic oscillators at higher flow rates. These results enable the prediction of fluidic oscillator performance, which can significantly improve the design process for an application where a fluidic oscillator may be used, such as aerospace applications, power generation, heat exchangers, or medical devices.

FEA of friction stir welding of AA1100 and parametric optimization using genetic algorithm

3D uncoupled thermo-mechanical analysis of friction stir welding (FSW) of aluminum alloy AA1100 was carried out using FEM. The heat generated in the FSW process is applied via DFLUX subroutine of ABAQUS. Simulations with different process parameters such as welding speed and rotational speed were performed for the estimation of temperature generated during the FSW process. Thermal analysis’s output is applied as thermal loads during mechanical analysis. Subsequently, residual stress developed during the FSW process was found out. Residual stresses were measured along the selected path, which included node to node stress variations in longitudinal directions. Maximum temperature appears to occur at the center of the weld pool and decreases away from the weld line towards the base metal. The numerical and experimental results were successfully validated with error in the range of 0.5 - 6%. Then, a semi-analytical equation for temperature distribution considering moving point heat source is proposed to carry out single objective unconstrained optimization of process parameters (welding speed and rotational speed). Finally, optimization was carried out to estimate the maximum temperature generated during the FSW process using a genetic algorithm in MATLAB.

A computational approach to determine key anatomic landmarks on pelvis and its application to acetabular orientation assessment and hip computational biomechanics.

The determination of anatomic landmarks on the pelvis is an important procedure in orthopaedics, computational biomechanics and morphometrics. This study aimed to develop a robust method for determining the anatomic landmarks on the pelvis and demonstrate the application of the method in assessing the native acetabular orientation and developing multiscale computational modelling for hip biomechanics. Two Matlab algorithms, named "Pelvisor" and "Acetabulor", were developed to determine four anatomic landmarks on the pelvis and assess the acetabular orientation in the natural hip. The algorithms were then applied to: (1) assess the native acetabular orientations for a cohort of subjects; (2) develop a multiscale computational modelling to investigate the biomechanical behaviour of the biphasic cartilage in the natural hip joint. Results showed that the present method can accurately determine the key anatomic landmarks on the pelvis and assess the acetabular orientation in the natural hip joint. Significant differences in acetabular inclination and anteversion between males and females were identified, with larger orientations in females. The interstitial fluid pressure in the cartilage supported over 90% of loading during gait, and the stresses and strains within the solid matrix of the cartilage were as low as 1.8 MPa and 0.12, respectively.

Location optimization of emergency medical facilities for public health emergencies in megacities based on genetic algorithm

PurposeThe purpose of this study is to discuss the principles and factors that influence the site selection of emergency medical facilities for public health emergencies. This paper discusses the selection of the best facilities from the available facilities, proposes the capacity of new facilities, presents a logistic regression model and establishes a site selection model for emergency medical facilities for public health emergencies in megacities.Design/methodology/approachUsing Guangzhou City as the research object, seven alternative facility points and the points' capacities were preset. Nine demand points were determined, and two facility locations were selected using genetic algorithms (GAs) in MATLAB for programing simulation and operational analysis.FindingsComparing the results of the improved GA, the results show that the improved model has fewer evolutionary generations and a faster operation speed, and that the model outperforms the traditional P-center model. The GA provides a theoretical foundation for determining the construction location of emergency medical facilities in megacities in the event of a public health emergency.Research limitations/implicationsFirst, in this case study, there is no scientific assessment of the establishment of the capacity of the facility point, but that is a subjective method based on the assumption of the capacity of the surrounding existing hospitals. Second, because this is a theoretical analysis, the model developed in this study does not consider the actual driving speed and driving distance, but the speed of the unified average driving distance and the driving distance to take the average of multiple distances.Practical implicationsThe results show that the method increases the selection space of decision-makers, provides them with stable technical support, helps them quickly determine the location of emergency medical facilities to respond to disaster relief work and provides better action plans for decision makers.Social implicationsThe results show that the algorithm performs well, which verifies the applicability of this model. When the solution results of the improved GA are compared, the results show that the improved model has fewer evolutionary generations, faster operation speed and better model than the intermediate model GA. This model can more successfully find the optimal location decision scheme, making that more suitable for the location problem of megacities in the case of public health emergencies.Originality/valueThe research findings provide a theoretical and decision-making basis for the location of government emergency medical facilities, as well as guidance for enterprises constructing emergency medical facilities.

The Measurement of Contact Angle, pH, and Conductivity of Artificial Saliva and Mouthwashes on Enamel, Glass-Ionomer, and Composite Dental Materials.

The aim of this study was to tackle the topic of appropriate recommendations for artificial-saliva and mouthwash usage. The contact angle, pH, and conductivity of two artificial saliva solutions, four mouthwashes, and their mixtures on enamel, glass-ionomer, and composite dental materials were measured. The measurements were conducted with a MATLAB algorithm to minimize human error. The obtained values for the contact angle were in the range from 7.98° to 52.6°, and they showed completely nonlinear and nonuniform behavior for all investigated liquids and on all investigated substrates. Results reveal statistically significant differences among all tested liquids on all investigated substrates (p < 0.05). pH values ranged from 1.54 to 7.01. A wide range of conductivity values were observed, from 1205µS/cm in the saliva-stimulating solution to 6679 mS/cm in the artificial saliva. Spearman’s test showed a moderate positive correlation between the pH and conductivity of the tested fluids (R = 0.7108). A comparison of the data obtained using Image J software and the MATLAB algorithm showed consistency, not exceeding 5% error. When an experiment uses human material and bioactive materials THAT are used in biomedicine as substrates, an additional definition of protocols is highly recommended for future research on this topic.

Radiative MHD flow of Rivlin–Ericksen nanofluid of grade three through porous medium with uniform heat source

BackgroundThe paper investigates the analytical and numerical solution of the radiation effect on MHD flow of Rivlin–Ericksen nanofluid of grade three through a porous medium with a uniform heat source between two vertical flat plates. The governing equations are solved analytically using multi-step differential transform method (MDTM) and numerically using finite difference method (FDM) and shooting method by designing MATLAB and Mathematica algorithms. The study discovered that the MDTM, FDM, and shooting methods are effective for solving nonlinear differential equations like this one.ResultsGraphs and tables show the influence of different parameters on velocity and temperature. Figures and tables show the comparisons between current outcomes and previous results that are accessible.ConclusionsThe present results showed that the analytical and numerical solutions agree well with previously published outcomes.

An Adaptive EWMA Control Chart Based on Principal Component Method to Monitor Process Mean Vector

The special causes of variations, which is also known as a shift, can occur in a single or more than one related process characteristics. Statistical process control tools such as control charts are useful to monitor shifts in the process parameters (location and/or dispersion). In real-life situation, the shift is emerging in different sizes, and it is hard to identify it with classical control charts. Moreover, more than one process of characteristics required special attention because they must monitor jointly due to the association among them. This study offers two adaptive control charts to monitor the different sizes of a shift in the process mean vector. The novelty behind this study is to use dimensionally reduction techniques such as principal component analysis (PCA) and an adaptive method such as Huber and Bi-square functions. In brief, the multivariate cumulative sum control chart based on PCA is designed, and its plotting statistic is utilized as an input in the classical exponentially weighted moving average (EWMA) control chart. The run length (RL) properties of the proposed and other control charts are obtained by designing algorithms in MATLAB through a Monte Carlo simulation. For a single shift, the performance of the control charts is assessed through an average of RL, standard deviation of RL, and standard error of RL. Likewise, overall performance measures such as extra quadratic loss, relative ARL, and the performance comparison index are also used. The comparison reveals the superiority over other control charts. Furthermore, to emphasize the application process and benefits of the proposed control charts, a real-life example of the wind turbine process is included.

Identification of precipitation trend and landslide susceptibility analysis in Miandoab County using MATLAB.

Iranian plateau is a seismically active region. Within this region, northwestern Iran is of high importance. Before proper planning for mitigating the earthquake-induced hazards can be achieved, it is necessary to identify high-risk areas in terms of susceptibility to earthquakes. In this study, landslide susceptibility in Miandoab Country was modeled using the so-called random forest algorithm (RFA) in MATLAB based on records acquired at 67 earthquake hotspots considering 9 factors affecting the earthquake occurrence (i.e., height, slope, direction, distance from fault, distance from river, distance from road, land use, geology, and precipitation). Predictive power of the model and validity of its results were evaluated using relative operating characteristic (ROC) curve and area under the curve (AUC). The assessment results showed very good accuracy of the model (0.97). It was further found that digital height layer, geology, and distance from fault impose the largest contributions into earthquake potential. The results also showed that 53%, 8.3%, and 38.4% of the study area were classified as being at low risk, moderate risk, and high risk of earthquake. Among other climatic elements, the precipitation exhibits the largest fluctuations; we proceeded to evaluate precipitation trends in the study area for a statistical period of 30years. This was practiced by implementing Mann-Kendall nonparametric test in MATLAB. This subject-matter is especially important in Iran where the annual precipitation level is as low as 250mm. The results showed that the precipitation follows an increasing trend in the region.

Analysis and Research on the Automatic Control Systems of Oil–Water Baffles in Horizontal Three-Phase Separators

The three-phase separator is one of the most important pieces of equipment in the combined station of the oilfield. The control level of the oil–water interface directly affects the energy consumption of the subsequent production of the combined station and the effect of oil, gas and water separation. In order to avoid these situations, the Siemens PLC control system, configuration software WinCC and MATLAB were used. The OPC technology is used to connect communication between WinCC and MATLAB, and the genetic algorithm in MATLAB is used to obtain the optimal separation height of the oil–water interface under the produced liquid in different periods. Subsequently, through the Siemens PLC system and WinCC configuration software, the automatic control of the three-phase separator is achieved, and finally the water content of crude oil is significantly reduced. The system provides a visual interface function. In the future, it will also provide an effective simulation platform for the theoretical research and design of an automatic control system of an oilfield combined station.

DIRECT COMPARISON OF INTERSECTING-TANGENTS VERSUS SECOND DERIVATIVE FOR DETERMINATION OF PULSE TRANSIT TIME AND PULSE WAVE VELOCITY

Objective: Arterial stiffness is a non-traditional risk factor for cardiovascular disease. Aortic stiffness is assessed by determination of pulse wave velocity using pulse transit time and the distance between carotid and femoral arteries. Transit time is obtained by using the foot-to-foot method to define the transit time through intersecting tangents algorithm or the point of maximal upstroke during systole (2nd derivative). Millasseau et al have proposed a formula for converting transit time between methods using SphygmoCor (intersecting tangents) and the Complior Analyse (2nd derivative). Based on a mathematical modeling of the proposed formula, there is discrepancy between values of pulse wave velocities, especially in subjects with higher aortic stiffness. The objective of this study is to directly compare the two methods using the same pressure waveforms obtained by the newer generation of Complior Analyse and using Millasseau’s formula. Design and method: In a cross-sectional study of heterogeneous subjects, aortic stiffness was assessed by the Complior Analyse device which uses 2nd derivative. The pulse waveforms were extracted and used for the analysis by custom MATLAB algorithm for intersecting tangents, and the results were compared to the formula proposed by Millasseau. Results: The preliminary results of the first 24 patients (men: 71%; mean age: 61 ± 18 years) show that Millasseau’s formula underestimates the transit times values by about 19% in comparison with the transit times obtained by the intersecting tangents method using MATLAB software (49,8 ± 18,8 ms vs 61,6 ± 18,6 ms; P &lt; 0,001). This results in an overestimation of the pulse wave velocities values by about 30% (13,8 ± 3,8 m/s vs 10,6 ± 2,7 m/s; P &lt; 0,001). Conclusions: Our preliminary results allow us to conclude that the values of pulse wave velocities obtained with Millasseau’s formula are overestimated values when comparing with values obtained by using the intersecting tangents method. Increasing the number of subjects will allow us to examine the possibility of a more reliable formula for converting transit times from one method to the other.

Inverse algorithm for extraction of multiple parameters using analytical model of eddy current response

Eddy current (EC) technology is commonly used for detecting flaws, measuring geometric parameters, or determining properties of conducting materials. However, the measurement of a particular parameter can become more challenging if multiple influential parameters vary simultaneously. In particular, eddy current-based measurement of separation (gap) between a pressure tube (PT) and a calandria tube (CT) in the fuel channels of CANDU® reactors is made more difficult by variations in PT wall thickness, resistivity, and probe lift-off. An analytical model of the EC response to changes in PT–CT gap has been developed by approximating the geometry of the PT within the larger diameter CT as a pair of concentric tubes, where gap is varied by changing the CT radius. In this article, this model is used in combination with an error minimization algorithm to construct an inverse algorithm for the extraction of PT–CT gap, PT resistivity (ρ), and PT wall thickness (WT) from measured multi-frequency eddy current signals. Application of a linear regression tool in MATLAB, with fourth-order polynomial fitting of modeled data with varying ρ and WT as a function of PT–CT gap, is used to obtain coefficients that depend on ρ and WT. Output of multidimensional fitting of these coefficients is scaled and rotated to calibration data. Finally, implementation of an error minimization algorithm in MATLAB is used to produce estimates of multiple target parameters from experimental data. Simultaneous extraction of either one, two, or three parameters is examined, using experimental data obtained at frequencies used for in-reactor inspection of 4.2, 8, and 16 kHz, or just two frequencies of 4.2 and 8 kHz. Under full gap variation conditions, the inverse algorithm predicts gap to within 0.1 mm at gaps between 0 and 9 mm and to within 0.4 mm at gaps between 9 and 18 mm. PT resistivity is predicted to within 1 μΩ cm (2% relative error) and PT wall thickness within 0.03 mm (1% relative error) when each is the only extracted parameter. An excellent agreement between actual and predicted values of gap demonstrates the potential of the inverse algorithm for application to in-reactor gap measurement and simultaneous extraction of either PT wall thickness or resistivity when the other parameter is known. The extraction of PT resistivity may be particularly useful, as this parameter cannot otherwise currently be measured in-reactor.

Environmental and social sustainability automated evaluation of plazas based on 3D visibility measurements

The purpose of this paper is to present an innovative method for measuring and evaluating 3D visibility indicators to design sustainable urban plazas. Sustainable urban plazas are studied in five physical, social, semantic, spatial, and ecological dimensions. Based on the literature review, automated 3D geometrical evaluation calculated by five normalized properties taken from a 3D space of Plaza included: Isovist, Smallness, Compactness, Enclosure, and Regularity. Human subjective evaluations are tested in 20 worldwide plazas. The research was conducted through survey methods and MATLAB software analysis, image processing, algorithm finding, machine learning, and inferential statics. Accordingly, the question was answered: how does the citizens’ mental evaluation change based on the squares’ three-dimensional characteristics? Findings show a significant relationship between people’s evaluation of the plazas and the three-dimensional geometric evaluations. The research results showed that three-dimensional geometric evaluations are a more accurate and complete method for recognizing high-quality plazas. Also, five normalized properties could explain 0.67% of the place quality and sustainability. Comparing this study’s findings with two-dimensional research shows a significant difference in evaluating the quality of two-dimensional and three-dimensional fields.

The Characterization of Liquefied Petroleum Gas (LPG) Using a Modified Bunsen Burner

The equivalence ratio ranges were found between 22.77 and 42.93 for the Saudi LPG/air mixture using a traditional Bunsen burner. An operation problem was found with a traditional Bunsen burner for the Saudi LPG/air mixture, especially in a lean mixture. Therefore, a Bunsen burner was successfully modified to overcome the limits of operation with different mixtures of Saudi LPG/air and a stable flame was obtained. The equivalence ratio ranges were found between 0.68 and 1.30 using the modified Bunsen burner. A premixed flame was used for the modified Bunsen burner. A MATLAB algorithm was successfully applied to flame image processing and measurement of laminar burning velocity. The laminar burning velocity was determined to be approximately 35 ± 0.91 cm/s under stoichiometric conditions using the modified Bunsen burner for the Saudi LPG/air mixture. The half-cone angle of the flame was found to be 16.20 ± 0.76°. The minimum flame height was observed to be 21.50 ± 0.22 mm above the Bunsen burner exit.